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Abstract 



The present review gives a survey of recent developments and applications of the Nambu- 
Jona-Lasinio model with Nf = 2 and Nf = 3 quark flavors for the structure of baryons. The 
model is an effective chiral quark theory which incorporates the SU(Nj)£(g>SU(Nj)^(g>U(l)y 
approximate symmetry of Quantum chromo dynamics. The approach describes the spontaneous 
chiral symmetry breaking and dynamical quark mass generation. Mesons appear as quark- 
antiquark excitations and baryons arise as non-topological solitons with three valence quarks and 
a polarized Dirac sea. For the evaluation of the baryon properties the present review concentrates 
on the non-linear Nambu-Jona-Lasinio model with quark and Goldstone degrees of freedom 
which is identical to the Chiral quark soliton model obtained from the instanton liquid model of 
the QCD vacuum. In this non-linear model, a wide variety of observables of baryons of the octet 
and decuplet is considered. These include, in particular, electromagnetic, axial, pseudoscalar 
and pion nucleon form factors and the related static properties like magnetic moments, radii 
and coupling constants of the nucleon as well as the mass splittings and electromagnetic form 
factors of hyperons. Predictions are given for the strange form factors, the scalar form factor 
and the tensor charge of the nucleon. 



Though this be madness, yet there is method in it. 

Shakespeare, Hamlet, II-2 
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1 Introduction 



The structure and the dynamics of hadrons are generally believed to be described by Quantum 
chromodynamics (QCD). A prominent feature of QCD is the asymptotic freedom which allows the 
processes at high energies to be described perturbatively in terms of these degrees of freedom. How- 
ever, at energies comparable to the low-lying hadron masses QCD shows non-perturbative phenomena 
such as confinement of quarks and gluons and the spontaneous breakdown of chiral symmetry. It 
makes the description in terms of quarks and gluons enormously complicated and the only direct 
source of information from QCD is the lattice Monte Carlo simulations. In this situation it is natural 
to invoke effective models of the strong interaction which incorporate the relevant degrees of freedom 
in the low energy regime. Those are all related to chiral symmetry and its dynamical breaking. In 
fact, the spontaneous breakdown of the chiral symmetry has been phenomeno logically recognized to 
be an essential feature of the strong interaction long before the advent of QCD. One has by now 
accumulated a large number of empirical facts and, following from this, somewhat generally accepted 
features of low energy hadronic dynamics, which have to be respected by effective models: 

• The existence of light pseudoscalar mesons (pions) and the fact that there is no degenerate part- 
ner of the nucleon with negative parity suggests that the pions are Goldstone bosons associated 
with the dynamical breaking of SU(2)^®SU(2)/j chiral flavor symmetry down to SU(2)y. This 
picture can be extended to SU(3), although obviously the larger masses of kaons and rj indicate 
that the explicit chiral symmetry breaking is larger in this case. 

• The relatively large mass of rj' suggests that the flavor [/^(l)-current is not conserved even in 
chiral limit. 

• The hadron mass differences in the SU(3)-flavor multiplets appear as being governed by a linear 
breaking of the SU(3)y-symmetry in a rather simple fashion and follow phenomenologically the 
Gell-Mann-Okubo mass formula. 

• The masses of the low lying mesons and baryons, apart from the Goldstone bosons, are quite 
well described in terms of the constituent quark model whose basic ingredients are massive 
constituent up and down (M Ujd ps 350 MeV) and strange (M s ~ 500 MeV) quarks. 

Actually the above empirical facts can be understood from the viewpoint of QCD. Its lagrangian is 
given by 

£qcd = + g-fAp - m)V - ^Sp V" . (1) 

Here \I> comprises a quark field with N c = 3 colors and Nf = 6 flavors, and and F^ v are the 
corresponding SU(3)-color gauge field and field strength tensor, respectively. 

The current quark masses are represented in (|l]) by the matrix rh = diag(m u , m^, m s , m c , rrib, m t ). 
Since up-, down- and strange current quark masses are much smaller than those of other quarks, 
there is a natural separation between "light quark physics", relevant for energies ~ 1 GeV, and 
"heavy quark physics" for higher energies. Thus for systems composed of up-, down- and strange 
quarks one can ignore the heavy quarks and work with Nf = 2 or 3. Accordingly it is physically 
reasonable to consider the limit of vanishing (light) current quark masses and if needed the current 
quark masses can be treated as perturbations. This limit corresponds to exact chiral symmetry. The 
latter is broken spontaneously which results in a non-zero chiral condensate (0|\I/\I/|0). The breakdown 
of chiral symmetry leads to the appearance of Goldstone bosons and to a dynamical generation of 
a constituent quark mass of a few hundred MeV. Although for many observables one can probably 



neglect m u and rrid, the finite value of m s m 175 MeV is too large to be ignored. In particular, the 
relatively large current strange mass m s is responsible for the baryon mass splitting summarized in 
the Gell-Mann-Okubo mass formula. 

The number of colors iV c plays a special role in QCD. The inverse number of colors 1/N C can be 
considered ('t Hooft, 1974) as a small parameter of QCD. Although in the real world N c = 3, the 
1/N C expansion is a very useful tool to analyze non-perturbative phenomena in the QCD. We will use 
this technique also in the effective theory considered below. In fact, concepts like that of constituent 
quarks and the idea of describing baryons as chiral solitons with a polarized Dirac sea emerge naturally 
from the large N c limit of QCD and have proven to be phenomeno logically quit useful. 

In this report we shall concentrate on theory and applications of an effective model, which on one hand 
is formulated in terms of quarks interacting via a simple and tractable quark- ant iquark interaction, 
and on the other hand fulfills basically the above mentioned phenomenological demands, in particular 
the dynamical chiral symmetry breaking and the explicit symmetry breaking due to current quark 
masses. The model has first been proposed by Nambu and Jona-Lasinio (1961a,b) in a different 
context. It has become quite popular and successful in a form where the fermion fields are re- 
interpreted as quark fields. In the meson sector it has already been the subject of several review 
articles (Vogl and Weise, 1991, Klevansky, 1992, Hatsuda and Kunihiro, 1994). In the baryon sector 
the most interesting aspect is the applications to physical observables of nucleon and hyperons and 
this is the subject of the present report. Previous reviews cover only partially these topics (Meissner 
et al., 1993, Alkofer et al., 1995). The Nambu - Jona-Lasinio (NJL) model provides an excellent 
account of the today's experimental data if both the collective quantization and the perturbative 
expansion in m s are performed properly. The model describes the physics in a way quite different 
from the other well known baryonic effective model, i.e. the Skyrme model (Skyrme 1961, 1962). 
In the Skyrme-approach (Adkins et al., 1993, Holzwarth and Schwesinger, 1986, Zahed and Brown, 
1986) one does not have quarks but meson fields and the baryons are described as topological solitons. 
In fact, the NJL model contains this physical picture as a limiting case, though being a genuine non- 
topological model. The relationship between these two approaches will also be a subject of discussion 
of our review. 

The NJL model in its simplest SU(iV/ = 2)-form is given by 

C = ^ 7 M <9 M - m )^ + ^-[(M) 2 + (^ 7 5T#) 2 ] . (2) 

Here \l/ are quarks fields, which are flavor doublets with N c = 3 colors, and itlq is the average up and 
down quark current mass. The corresponding semibosonized lagrangian is 

£ = - a - «r • r 75 )^ - ^(a 2 + vr 2 ) + ^-a , (3) 

where a and 7? are auxiliary meson fields. The model is not renormalizable and lacks confinement. 
Although these features were for many years a reason to ignore the NJL model, nowadays the view is 
different. For an effective model in the low energy regime, the problem of non-renormalizability can 
be avoided if one has a reasonable means to fix a suitable ultraviolet cut-off. The lack of confinement 
provides in NJL model unphysical ^^/-thresholds. However, there are strong arguments that the 
basic properties of light hadrons can be reproduced without invoking confinement and one of them is 
the recent lattice QCD results for the correlation functions of hadronic currents in the cooled vacuum 
(Chu et al, 1994). This point of view will be also supported by the results reviewed in the present 
paper. 

In the NJL model, baryons appear as large N c solutions consisting of valence quarks bound to the 
Dirac sea via classical (mean-field) meson fields in accordance with the arguments of Witten (1979). 



This physical picture was proposed by Kahana and Ripka (1984), and by Diakonov et al. (1986). 
However, the mean-field solutions exist only if the classical meson fields are subjected to the condition 
of the chiral circle 

a 2 + v? 2 = M 2 , (4) 

where M is the dynamically generated constituent quark mass. Under this constraint the NJL 
lagrangian (|^) reduces to an effective lagrangian 

C! = ^{i^d^ -mo- MC/ 75 )# , (5) 

in which the quarks interact only with Goldstone fields and the SU(2) matrix U corresponds to a non- 
linear representation of the pion field MU 15 = a + I^ttt. For SU(3) the chiral field configuration U 
also encompasses kaon and eta fields. In the actual calculations U will be determined selfconsistently 
by minimizing the effective action which corresponds to C @. The model with constraint (§]) is 
usually referred to as the non-linear NJL model or, when applied to baryons, the Chiral quark soliton 
model. 

The assumption of the chiral circle (non-linear pion field configuration) is not an ad hoc device to allow 
for solitonic solutions of NJL. Physically, it is a means of separating low-energy mesonic degrees of 
freedom (Goldstone bosons) from heavier ones (sigma meson) and to "freeze" the latter. The physics, 
which apparently goes beyond the scope of the NJL model, can be understood in the framework of 
the underlying theory, the low-energy QCD. 

There are various attempts (McKay and Munczek, 1985, Cahill and Roberts, 1985, Diakonov and 
Petrov, 1986, Ball, 1989, Schaden et al., 1990) to relate the NJL model to some low-energy limit of 
QCD. Further, arguments to illuminate the chiral structure of the NJL model are given by Ebert 
and Volkov (1983), Volkov (1984), Dhar and Wadia (1984) and Ebert and Reinhardt (1986). The 
relationship to the chiral perturbation theory (Gasser and Leutwyler, 1982) is studied by Hansson et 
al. (1990), Bernard and Meissner (1991), Ruiz Arriola (1991), Schiiren et al. (1992) and Bijnens et 
al. (1993). In fact, the physics behind the present model can be well understood in the framework 
of the instanton liquid model of QCD (Shuryak 1982, 1983, Diakonov and Petrov, 1984, 1986) which 
is also supported by recent lattice QCD calculations (Chu et al, 1994). The assumption that the 
non-perturbative ground state of QCD is governed by a dilute liquid of interacting instantons and anti- 
instantons allows to obtain in large N c limit an effective low-energy theory (Diakonov and Petrov, 
1986) which can be used to justify the NJL lagrangian on the chiral circle, i.e. eq. ([5]). In fact, 
by assuming the instanton model of the QCD vacuum, Diakonov and Petrov (1984, 1986) derived 
an effective low-energy QCD lagrangian in terms of a 't Hooft-like 2iVj-fermion interaction. The 
crucial quantity in the instanton vacuum is the ratio p/R of the average size p of instantons to the 
average distance R between them. Both the phenomenological analysis (Shuryak 1982, 1983) and 
the variational estimates (Diakonov and Petrov, 1984) suggested a relatively small value p/R ~ 1/3 
which has been confirmed by the lattice calculations (Chu et al, 1994). Thus, besides the inverse 
number of colors 1/N C this small ratio plays also an important role of a small parameter in the 
theory. If one places light quarks in the instanton vacuum, one observes that chiral symmetry is 
dynamically broken due to the derealization of the would-be fermion zero modes related to the 
individual instantons. Spontaneous chiral symmetry breaking manifests itself in a non-zero chiral 
condensate (0|\i/\I/|0), massless Goldstone bosons, and in a momentum dependent dynamical quark 
mass M(k) with M(0) « 350 MeV and M(k) vanishing for momenta k 3> 1/p. Due to the instanton- 
induced 't Hooft-like interaction the rf is not a Goldstone boson, so that the ?7a(1) problem is 
properly solved. The momentum dependence of M(k) provides a natural cut-off 1/p which makes 
the resulting theory finite. It is important that the dynamical quark mass is parametrically small 
Mp = 0(p 2 /R 2 ) <C 1 compared to the cutoff 1/p. As a final result, due to the diluteness parameter 
p/R in this instanton model of QCD vacuum, at low momenta k < 1/p and large N c the QCD is 
reduced to an effective quark meson theory with a lagrangian which includes only constituent quarks 



of mass M(k) interacting with Goldstone pion fields. If both quark and pion momenta are smaller 
than 1/p the dynamical mass M(k) can be approximated by M = M(0). In this approximation one 
reproduces the semibosonized NJL lagrangian (|5|). Since the inverse size of nucleons corresponds to 
~ 300 MeV this can be considered as "much" less than the cutoff 1/p and the effective lagrangian (|5|) 
can be also applied to describe baryons. Hence a systematic treatment of the nucleons and hyperons 
involves only quarks and Goldstone bosons, and the implementation of the chiral circle in (|3[) is 
justified. It should be noted that the only relevant degrees of freedom are pions and quarks and 
excitations connected with other degrees of freedom like vector mesons have k ~ 1/p and do not 
contribute essentially to the dynamics at momenta k <C 1/p. Following this logic one realizes that 
for a proper description of the baryon ground-state properties the only relevant degrees of freedom 
are those of quarks and pions, and heavier meson degrees of freedom like scalar and vector mesons 
are not important. 

The present paper concentrates on the results for the baryon properties in the non-linear NJL model. 
The paper reviews applications in flavor SU(2) and SU(3) versions to electromagnetic, axial, scalar 
and strange vector form factors and the related static properties like magnetic moments, radii and 
coupling constants of the nucleon, and to the mass splittings and electromagnetic properties of the 
hyperons of the octet and decuplet. Calculations of the spin properties are discussed in connection 
with the "spin crisis". Predictions are reviewed also for nucleon strange form factors and tensor 
charges. For all observables, the calculations reported are performed with one and the same set of 
parameter values. These values are fixed in the meson sector by fitting the physical values of the 
pion mass and the pion decay constant, and additionally to the kaon mass in the case of SU(3). 
The only parameter, which is free, is the constituent mass M. However, as we will show for the 
constituent mass of about M = 420 MeV the model is rather successful in describing the baryon 
properties. Actually, the comparison of the model results with the experimental data of the nucleon 
and also of the hyperons in the octet and decuplet shows that the NJL model on the chiral circle 
indeed is able to reproduce most of the data within 15 %. This includes critical observables like the 
axial vector coupling constant, g^, and the magnetic moments of proton and neutron. Only a few 
observables are less successfully described (within 30 %) as e.g. the neutron electric form factor and 
the corresponding charge radius. 

Since the present theory yields a rather good agreement with the experimental data one can use the 
model picture to obtain some insight into the low-energy structure of the nucleon and the other octet 
and decuplet baryons. The basic mechanism is obviously the interaction of quarks with pseudoscalar 
Goldstone bosons, the latter being themselves quark- ant iquark excitations of the spontaneously bro- 
ken chiral vacuum. Heavier mesons like scalar mesons or vector mesons are apparently not needed 
in this scenario. In addition, the confining forces seem also not to be much relevant in describing 
the ground-state properties of low mass mesons and baryons. Hence, the only important mechanism 
seems to be chiral symmetry breaking. The calculations support the physical picture of baryons as 
non-topological solitons of N c = 3 valence quarks bound to a polarized Dirac sea via classical (mean- 
field) meson fields and being well separated from the negative Dirac continuum. The polarized Dirac 
sea phenomenologically corresponds to a meson cloud. Most of the observables receive about 30 % 
contributions from the Dirac sea quarks. This valence picture of the baryon in the present model 
differs from the one of the topological solitons favored in the Skyrme model. In principle, the NJL 
model can be reduced to a topological approach similar to the Skyrme model by an expansion of the 
effective action in terms of gradients of the meson fields (gradient expansion). This expansion does 
only converge if the valence level is close to the negative Dirac continuum. This, however, requires in 
the NJL model unphysically large values of the constituent mass for which the baryon properties are 
not reproduced. Hence, within the NJL model there is no way to recover Skyrme type models. 



2 SU(2) NJL model: vacuum sector and meson properties 



In this chapter we deal with the Nambu - Jona Lasinio model with two lightest quark flavors. To 
begin with, we study the NJL lagrangian and its symmetries. Using the functional integral approach 
we bosonize the theory and solve the effective bosonized theory in the mean-field approximation which 
corresponds to the leading order of the expansion in inverse number of quark colors. In the functional 
integral approach the mean-field solution naturally appears as a saddle point of the bosonized action. 
We analyze the structure of the ultraviolet divergences of the model and introduce the regularization 
scheme. We show that in the mean-field treatment of the model chiral symmetry is spontaneously 
broken, which leads to a generation of the constituent quark mass. An important consequence of 
the spontaneous breakdown of chiral symmetry is the appearance of light Goldstone mesons (pions). 
We discuss shortly the description of mesons in the model and sketch the derivation of the meson 
propagators, masses and of the meson decay constant needed to fix the model parameters. 



2.1 Bosonization of the NJL model 



In the simplest SU(2) version it includes only scalar and pseudoscalar couplings: 



NJL lagrangian and its symmetries. The NJL lagrangian contains a local four-fermion interaction. 

(6) 



C = V(x) 



i0 — m 



G 



X 



2i 



where one recognizes the free and the interaction parts 

iff(x) 



£free 

G 
~2 



%0 — m 



C 



int 



^(x)^(x)) 2 + Nf(x)ij 5 T^\ 



X 



(7) 



respectively. Here ty(x) denotes a quark field with u and d flavors and N c colors. The matrices t° 
(a = 1, 2, 3) are usual Pauli matrices. The rh represents the current quark mass matrix given by 



m 



m u 
m d 



m l + r 3 m 3 



It is useful to divide it in the average mass itlq and in the isospin mass difference m 3 defined as 



m = -(m u + m d ), m 3 



m u - m d 



(8) 



(9) 



In @ G is the coupling constant of the four fermion interaction. 



Lagrangian ([]) possesses a number of symmetries which we also have in QCD. In particular, it is 
symmetric under the global U(l) transformation 



-> ty'(x) = exp(iS)^(x) . 



(10) 



If one neglects the quark isospin mass difference m 3 and considers m u = m d = mo, which corresponds 
to exact isospin symmetry, lagrangian (Q) is also invariant under the transformation 

-> V(x) = exp(zf- p)^(x) . (11) 



If we additionally neglect the current mass mo, it becomes invariant under global chiral transforma- 
tions: 

ty(x) -> V(x) = exp(ry 5 r • k)^(x) . (12) 

Thus, in chiral limit m u = m d = 0, the NJL lagrangian has an SUl(2) ® SUr{2) ® Uy{l) symmetry. 
Therefore, many well-known results obtained within the current algebra can be reproduced in the 
NJL model. Later we will show that chiral symmetry is spontaneously broken which leads to such 
fundamental consequences as the appearance of Goldstone bosons (pions) and to the generation of 
the dynamical quark mass. 



The above transformations correspond to the following Noether currents 

B^(x) = ^f(x) r y^(x) (baryon current) , (13) 



T a 



ft 

whose divergences are given by 



V£(x) = ty(x)'y fl —'fy(x) (vector current) , (14) 

A^(x) = ^(x)jfj,j 5 —^(x) (axial current) , (15) 



0"B„ = O, (16) 



T 



= i^[rh,— (17) 



2 



9^ = ztt 75 {m,y}*. 



The baryon current and the third component of the vector isotopic current VJf are conserved for any 
values of the quark masses. For the conservation of V^, we need equal quark masses m u = m d , 
whereas the axial current is conserved only in chiral limit m u = = 0, so that the last equation 
(p~8|) expresses the partial conservation of the axial current (PCAC). 



Below for simplicity we work with the NJL model in the approximation m u = m d = ttlq. How- 
ever, most of our formalism can be easily generalized to the case m u ^ at least in the linear 
approximation in small quark masses. 

The NJL lagrangian is assumed to mimic the low energy effective theory of QCD with gluon degrees of 
freedom integrated out. However the quark fields carry a color index that runs from 1 to N c (N c = 3 
in the real world). Instead of the SU(N C ) local gauge invariance of QCD, in the NJL lagrangian, 
we have only the global SU(N C ) invariance. Although the color degrees of freedom enter the NJL 
lagrangian in a rather trivial way, their role will be very important for the justification of the mean- 
field approximation to the NJL model. Later we shall show that the mean-field approach becomes 
asymptotically exact in the limit of large number of quark colors N c . 



Auxiliary boson fields. The NJL lagrangian ffl) contains only quark degrees of freedom. On the other 
hand, the light mesons which we are going to describe by using this effective lagrangian are pions. 
Therefore, it is desirable to find an equivalent formulation of the theory that involves meson fields 
and not only quarks. Another reason why we need the bosonization is that we are going to solve 
the model by using the mean field approximation. In terms of the functional integral approach the 
mean field solution corresponds to the saddle-point approximation. The saddle-point solution can be 
obtained only in terms of boson fields so that the bosonization of the NJL lagrangian is a natural 
step towards the mean-field solution of this model. 



Historically, the bosonization has first been formulated in solid state physics where it is known as 
Hubbard-Stratononovich transformation. In the present theory it has been first applied by Eguchi 
(1976), Kikkawa (1976) and Kleinert (1976). 

For simplicity we describe the bosonization of NJL on the level of the partition function 

V^fVWf d4xC{x) (19) 
keeping in mind that the generalization to other quantities is straightforward. 
Under proper normalization of the integration measure one can write the following identity 

1 = JvaVnexp^-t J d A x^ (a + G - + (t? + G^if 7s ^) 2 j (20) 

with a gaussian path integral over auxiliary boson fields a and tt on the rhs. Inserting this identity 
into the partition function Z we obtain: 

Z= [vVVVVoVTre'S 6 **^ , (21) 

where the path integral runs over both the quark and auxiliary meson fields. The corresponding 
semibosonized lagrangian is 

d = -a -it- f 75)* - ^ 2 + n 2 ) + ^-a . (22) 

One should keep in mind that on this stage the introduced auxiliary fields a and n are non-dynamical 
fields and no kinetic term \ [d^ad^a + d^nd^rr] appears in lagrangian (|22|) . However, after integrating 
out the quarks in (|2T|), i.e. after including the quark loop effects, the fields will be "dressed" and can 
describe the physical mesons. It should be stressed that since the bosonization procedure is exact, 
the semibosonized lagrangian CJ has the same symmetries as the initial lagrangian (|6]). Note that the 
meson fields a, 7? transform under the chiral rotation (12) as follows 

(0 -f- in ■ t) — > (a + in ■ r) = exp(— ir ■ k) (a + itf ■ f) exp(— if ■ k) . (23) 



Diagrammatically, the described bosonization procedure corresponds to rearrangement and resumma- 
tion of the graphs of the 4-fermion point interaction into quark-meson vertices of the Yukawa type by 
introducing collective scalar-isoscalar (a) and pseudoscalar-isovector (jf) auxiliary boson fields, both 
carrying the quantum numbers of the composite operators fifty) and fififj^ty) but not color. The 
idea has a long history and traces back to the sixties (see for instance (Coleman, 1985) and references 
therein). 



Effective action . Now instead of the four-fermion interaction we have a theory of quarks interacting 
with boson fields. Next we notice that the integral over quark fields in eq. (|2T|) is gaussian and can 
be done exactly, which leads to a determinant of the Dirac operator 



V^Vty exp 



: j d 4 xty(i$ - a - m ■ f 75)^ 



Det c {i(j) — a — in ■ 7^75) 

exp [iV c Tr log(i# -a -in- f 75)] . (24) 



Here the power of N c appeared due to iV c colors of the quark fields. 



The Dirac operator can be rewritten in the form 

ift — a — in ■ r 75 = -y [id t — h(a, tt)] . (25) 
where h(a, 7?) is the one-particle Dirac hamiltonian 

h(a, tT) = -i7°7 fe V fe + 7V + wr ■ r 75 ) . (26) 

Below we prefer to work with the euclidean time r and perform the Wick rotation 

r = it . (27) 

In the euclidean case the following change should be done in (|24]) 

Det(i$ -a -in- f 75) = Det[id t - h(a, tt)] -> Det[<9 r + /i(cr, 7?)] . (28) 

Integrating out the quarks as it is done in (p4[) we arrive at the following completely bosonized 
euclidean effective action 

S(a, 7?) = -iV c Tr log £>(<r, 7?) + ^ J & 4 x(a 2 + vf 2 ) - ^ | rf 4 a;a . (29) 

Here we use the notation 

D(a,if) = d T + h(a,n) . (30) 



Large iV c and mean-field approximation. In terms of the effective action ( p9|) the partition function 
of the model is given by equation 

Z = J VaVne- sia '^ . (31) 
All our formulas, which have been derived so far, are exact. However, in order to evaluate the 



functional integral ( pT|) we need some approximation. Note that the fermion determinant contribution 
to the effective action ([29]) has a prefactor of iV c . For a consistent large N c treatment of the effective 
action (]29f) we take the coupling constant of the NJL model to be of order 

G = OiN' 1 ) . (32) 

Then all terms in ( p9|) will be of order 0(N C ) and in large N c limit the functional integral (|3~TD can be 
evaluated in the saddle-point approximation. The stationary meson configuration a c , tt c minimizing 
the effective action S(a, 7?) can be found by solving the saddle point equations 



SS(a, 7?) 5S(a,jr 

07T 



5a 



TT — TYc W " 7T = lTc 



. (33) 



Thus, in the leading order of the 1/N C expansion the functional integral ( |3~TD is simply given by the 
value of the integrand at the saddle point. One can show that this approximation is nothing else 
but the Hartree mean-field approximation in the original 4-fermion version, which has been used 
repeatedly (Bernard et al., 1984, Bernard, 1986, Ferstl et at, 1986, Bernard et al, 1987, Providencia 
et al, 1987, Bernard et al, 1988, Klimt et al, 1990) and has been reviewed by Vogl and Weise (1991) 
and Klevansky (1992). 



2.2 Divergences and regularization 



Like the Fermi theory of weak interaction the Nambu-Jona-Lasinio model is non-renormalizable, 
because the coupling constant of the 4-fermion interaction G has the dimension [G] = mass~ 2 . This 
means that in each increasing order in G new graphs with a higher degree of ultraviolet divergence 
appear. In order to get a well defined theory it is therefore necessary to specify how the infinities of 
the model have to be treated. 

One faces the ultraviolet divergences in the effective action (|29|) already in the leading order in N c 
which are due to the Dirac determinant Det D(a, tt). It is convenient to consider separately the 
divergences of the real part 

Re Tr log D = - Tr log {D ] D) , (34) 

as well as of the imaginary part 

Im Tr log D = — Tr log (D / D ] ) . (35) 



UV divergences. Let us start from the real part of the effective action fl34l) . The functional determi- 
nant Det(DW) is a complicated non-local functional of fields a, tt and the simplest way to study its 
divergences is to consider the case of slowly changing fields a, tt and to expand this functional deter- 
minant in a series of integrals of local polynomials in derivatives of fields a, n (Aitchison and Fraser, 
1984, 1985a, 1985b). It should be noted that in the case of the proper-time regularized effective action 
(see next subsection eq. (pEIp) for a systematic expansion it is more convenient to use the heat kernel 
method (Kleinert, 1976, Nepomechie, 1985, Ebert and Reinhardt, 1986). However, the structure of 
divergent terms can be determined by a simple dimensional analysis without explicit calculations. 
Since the effective action is dimensionless and the divergent terms should contain a positive power 
or the logarithm of the cutoff, the local integrals of the fields and their derivatives in the divergent 
terms should have a dimension of mass m where m > 0. Combining this dimensional argumentation 
with the chiral symmetry it is easy to see that the divergent terms are of the form 

Re Tr log D = c x J <$x{o 2 + vr 2 ) + c 2 J d^x[a 2 + tt 2 ) 2 

+ c 3 J d A x[(d»a){d,a) + {d^){d^)\ + 5>, vr) . (36) 

Here ci, C2, C3 are some divergent coefficients and S(a, 7?) stands for the ultraviolet finite part. Ap- 
parently, in order to examine the particular divergent structure of coefficients q one should perform 
explicitly a gradient expansion of the real part of the effective action (0). Note that among the 
divergent terms there appeared the kinetic term for meson fields a, n. Later we will introduce a regu- 
larization so that this kinetic term will be finite but non-vanishing. The presence of this kinetic term 
means that the mesonic fields, which initially have been introduced as auxiliary unphysical fields, 
become dynamical fields after taking into account the quark loop effects. 



Imaginary part of the effective action. In the case of the 577(2) flavor group the imaginary part of 
the effective action vanishes. It is connected to the fact that operators D(a,jr) and D*(a,7r) (the 
asterisk stands for usual complex, not hermitian, conjugation) are unitary equivalent in the SU(2) 
case 

D*(a, 7?) = Vr 2 D(cr, vr)(\/r 2 ) t . (37) 



Here V is some unitary spin matrix defined by 

lt = Vi,V- 1 Gu= 1,2, 3, 4, 5), (38) 
where the asterisk stands for complex conjugation and 7 M are euclidean hermitian Dirac matrices. 

Using equality ([37]) one immediately concludes that 

Det {D/D ] ) = 1 , (39) 
which means that the imaginary part of the effective action vanishes exactly in the SU(2) case. 

In the case of higher flavor groups the imaginary part of the effective action is not necessarily vanishing. 
However, it still remains finite. The simplest way to check it is again to consider the gradient 
expansion. 



Regularization scheme. The NJL model is a non-renormalizable effective theory which should be 
considered as a low-energy approximation to QCD and hence is physically relevant only up to some 
momentum scale. Actually, in the instanton liquid model (Diakonov and Petrov, 1986) the inverse 
size of the instantons 1/p plays the role of a natural cutoff in the theory free from the UV divergences. 
Therefore, we have to regularize the present theory introducing a cutoff at this scale. There are many 
ways to do it and there is no reason to prefer some scheme in favor of the other provided certain 
constraints such as Lorentz invariance, gauge invariance, current conservation, integer baryon number, 
are obeyed (Ball, 1989). We have already shown that only the real part of the effective action contains 
ultraviolet divergences and hence needs regularization. In general, one is free to regularize the finite 
imaginary part as well. However, it would lead to additional complications related for instance to 
the definitions of the baryon number and the charge. Some insights into the problem can be gained 
from the underlying non-local theory (Diakonov and Petrov, 1986) which because of the momentum 
dependent dynamical mass M(k) is finite. In fact, Diakonov and Petrov (1986), and Ball and Ripka 
(1994) have shown that the results for the amplitude of the anomalous processes 7r° — > 27 in the 
non-local theory are very similar to those of the NJL model with non-regularized imaginary part of 
the effective action. 



In this review we regularize only the divergent real part and keep the finite imaginary part unaffected. 
We will discuss the results obtained in the proper-time scheme based on the identity 
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(40) 



In the case of our interest K = D^D and the ultraviolet divergence appears as a singularity in the 
integration over small u. We make this integral finite by simply introducing a cutoff A" 2 in the lower 
limit 
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Here we subtract the fermion determinant of the "free" Dirac operator Dq = D(ao, v?o) corresponding 
to some fixed fields <jq, ttq. One can generalize the proper-time regularization in the form 



^Tr log (rfD) reg - ^Tr log (dJd ) p 
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e -uDlD 



(42) 



where <p(u, A) is a function properly chosen to make the integral finite. In particular, for the specific 
choice 4>(u) = 0(-^ — u) one recovers (0). Different cutoff schemes have been considered in detail by 



Meissner Th. et al. (1990) for the vacuum, and by Blotz et al. (1990) and Doring et al. (1992a) in the 
soliton sectors. In both regimes the results of using other regularization schemes, Pauli-Villars and 
both sharp and soft 0(3)- and 0(4)-cuttofs, were rather similar to those of the proper-time scheme 
if the real part was only regularized. 



2.3 Spontaneously broken chiral symmetry and constituent quark mass 



The main observation of Nambu and Jona-Lasinio (1961a,b) in their original work is that if the 
coupling constant G is larger than a certain critical value, the four fermion interaction leads to a non- 
trivial vacuum solution in which chiral symmetry is spontaneously broken and in this way fermions 
acquire a dynamical mass. This was done in the canonical formalism using a Bogoliubov-Valatin 
transformation from bare massless fermions to constituent massive fermions as quasiparticles. This 
phenomenon is called dynamical mass generation. In the present section we will see explicitly how 
the spontaneous chiral symmetry breaking is realized within the path-integral approach to the NJL 
model in large N c limit. 



We look for the vacuum saddle-point solution of eqs.(33) in the form 

o c 7^ 0, 7? c = . 



where the expectation value a c is a constant field. Using non-regularized effective action 
saddle-point equation (|33| ) reads 



5S(<j, 7?) 
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(43) 
) the 

(44) 



In the case of the proper-time regularization of the fermion determinant ([IT]) the divergent integral 
I\{M) in the saddle-point eq. (|44T) should be replaced by: 
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(45) 



Eq. (f44|) is a non-linear equation for a c which is known as gap equation. In fact, it corresponds to the 
Schwinger-Dyson equation in the mean-field Hartree approximation. 



The non-zero value a c enters the Dirac operator fl26[ ) in the same way as the quark current mass and 
as such plays a role of a quark mass which is usually called constituent quark mass 

M = a c (46) 



The gap equation is illustrated in fig. |TJ. Because of the quark-loop contribution the quark line gets 
"dressed" and the quarks become massive even in chiral limit. As a result in the Dirac quark spectrum 
there appears a mass gap of 2M between the positive and negative energy continua. 

Note that the non-vanishing saddle-point expectation value a c ^ is not invariant under the chiral 
transformation ( |23"D even in chiral limit. This means that in the saddle-point solution chiral symmetry 
is spontaneously broken and the non-zero value for the constituent quark mass (|46D is the consequence 
of this spontaneous breakdown of chiral symmetry. 

Analogous to the magnetization in ferromagnets an order parameter for the non-trivial vacuum can 
be introduced. In fact, it is the quark condensate (^ r \l / ) = (uu) + (dd), which characterizes the 



+ 




Figure 1: Diagram representation of the gap equation. The thin line corresponds to "bare" quarks 
whereas the thick line represents the constituent "dressed" quarks. 



strength of the spontaneous breakdown of chiral symmetry, and as such plays the role of the chiral 
order parameter. In the leading order of the large N c approximation we obtain 
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—S(a c , 7r c ) - — 



(47) 



Here V4 is the four-dimensional euclidean space-time volume. Note that the additional term m\j (2G) 
appearing here should have been written already in the semibosonized lagrangian (^) but we have 
dropped it there since it is a field independent constant. Differentiating the bosonized action S(a c , n c ) 
with respect to m at the saddle point a c = M, n c = we get: 



(48) 
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< > = - -(a c -m ) = - — (M - mo) • 

Using explicitly the gap equation (Ell) we can rewrite this result in the form 



< ^ > = -8iV c MJi(M) -> —8N C M 



00 du 



-uM z 



(47r) 2 A- 2 u 



(49) 



One realizes that the condensate < > is a quadratically divergent quantity and depends strongly 
on details of the regularization scheme. 



2.4 Mesonic properties 



The evaluation of the mesonic two-point functions within the NJL model has been undertaken by 
several authors, both in a pure fermionic theory using a Bethe-Salpeter formalism in the ladder 
approximation (Blin et al, 1988, Bernard et at, 1988, Bernard and Meissner U.-G., 1988, Klimt 
et al, 1990, Takizawa et al, 1990) and in a bosonized version (Jaminon et al, 1989, 1992) in the 
path integral formalism. In this section we sketch the calculation of mesonic spectra in the bosonized 
version of the model as it has been done by Jaminon et al. The advantage of their calculations is 
that using path integrals one is able to start directly from the regularized theory. 

The mesonic two-point function K a b can be obtained from the generating functional 

Z(j) = J V<p a e~ s+ ^ 

with explicit meson source j a included 

6HnZ 
Kab{x - y) = — 

Here we use an abbreviation j ■ (f> — f d 4 x j a (x)(f) a (x) where <p a = (a, 7?) stands for both meson fields. 



(50) 



(51) 



In the model the mesons appear as low-lying collective ^f^f excitations. In order to describe these 
modes one should consider meson fluctuations around the stationary meson field configuration: 



+ 



(52) 



In the leading order in N c , the saddle-point value 0°, which minimizes the effective action 
coincides with the "classical" value < 4> a > 

< (f) a (x) > 



5\nZ 



Sj a (x) 



(53) 
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Since the meson fluctuations are small (lowest excitations) we can expand the effective action (|2 
up to the second order in the fluctuating meson fields (small amplitude approximation): 



Z(j) 

where we use the short notation 
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In the exponent we have a bilinear form and the integral can be easily evaluated 
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(56) 



For the last term in (|56|) we use the same short notation as ( |55"1) and Tr is the functional trace over 
the meson degrees of freedom. The first term in rhs is the effective action in leading order in N c fl29|) . 
The third one is the one-meson loop contribution. Compared to the leading term 5(0°) ~ iV c it is 
suppressed by 1/N C . It means that it is parametrically small (large N c ) and we do not include it in 
our considerations (zero-meson-loop approximation). From eq.(|56|) one also realizes that the inverse 
meson propagators in leading order in iV c are given by the second variation of the effective action S 
with respect to the meson fields at the stationary point 0o : 



Kl L (x-y) 
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(57) 



The particular calculations of the r.h.s. of the above equation are straightforward but quite lengthy 
and we skip it here. We refer the reader to (Jaminon et al., 1989, 1992) for details. After Fourier 
transform to the momentum space the final result for the meson propagator in leading order in iV c is 
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(58) 



GMZ p (q 2 ) 



The function Z p (q 2 ) corresponds to a quark-loop with two pseudoscalar-isovector insertions (275 r a ). 
Here we present only the proper-time regularized expression for it: 
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(59) 



The Pauli-Villars regularized expression can be found in (Schuren et al., 1993). The on-shell meson 
mass corresponds to the pole of the meson propagator 
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(60) 



From (HD one can see that in chiral limit mo — > 0, the pions become massless Goldstone bosons 
whereas m a = 2M. The meson coupling constants are given by the residue of the propagator at the 
pole 

gl= lim 2 (g 2 + mj)^(g 2 ), (61) 



n O 



and as usual the physical meson fields are defined as 

T a 



7T a 
Tip h = — and cr p h = — • (62) 

9n 9<r 



Another important quantity, which we need to fix our parameters, is the pion decay constant. It is 
defined as the matrix element of the axial current between the physical vacuum and a physical pion 
state of a four momenta q: 

< | A%(z) | rc\q) >= iq»Ue-^ z 5 ab . (63) 

Using the Lehmann-Symanzik-Zimmermann reduction formula one can write the pion decay constant 
as 

U = f(q 2 = -ml) . (64) 

which is given by 

^ f 4\^4^ e ~ iqZ =< I T K(z)7T b ph (0)) | > , (65) 
J (2-7r) 4 q z + m~ h ' 

where T stands for time ordering. In order to calculate the rhs of eq.(p5D we couple an external axial 
source J" 

D (<7, tt) — > D(a, vr) + J^zh^ . (66) 

in the effective action ( p9| ) and in eq.(^). Using the generating functional Z which now explicitly 
depends on we can write 

<0|71W*(0)]|0>=£ (67) 



From ([56]) one can calculate directly 



SJ;{z)Sj b (0) J 5J°(z)SMx) 
and combining eqs. (|65|) and (|68|) one gets 

q »fW = + / dS^l^e- . (69) 

The calculation of S jl$ is almost identical to one of the inverse pion propagator (|57D . The final 
result for the physical pion decay constant reads 

U = M 97T Z p (q 2 = -ml) , (70) 

where the regularized function Z p (q 2 ) is defined in (|59D . In chiral limit the pion coupling constant 
(|6~lD is given by 

^m- (71) 

Inserting this result in (|70f) one recovers the Goldberger-Treiman relation on the quark level: 

Ug* = M. (72) 

In chiral limit using the above relation and the explicit expression for Z p (0) ([59]) we get for the pion 
decay constant a simple expression 
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Figure 2: Cutoff A on the left and the quark condensate (uu) (solid line) and current quark mass mo 
(dashed line) on the right in the case of the proper-time regularization as a function of the constituent 
quark mass M. 

2.5 Fixing of the parameters in the vacuum sector 



The SU(2) NJL model, treated in the leading order in N c (no meson loops), contains three free 
parameters. Two of them, the coupling constant G and the quark current mass mo, are presented 
in the NJL lagrangian (|6|). The third one, the cutoff A is needed to make the theory finite. The 
parameters mo and A are fixed to reproduce the physical pion mass m n = 139 MeV and the physical 
pion decay constant f n = 93 MeV. As usual, using the gap equation ((PJ), the last model parameter, 
the coupling constant G, is eliminated in favor of the constituent quark mass M. 

Thus, starting with a particular value for M one fixes the cutoff A reproducing the physical value of 
the pion decay constant (f70|). For the pion mass we use the equation obtained from ( |58|) and (|60| ) 



and also the gap equation ([44]): 



m m < > 



GMZ p (q 2 = -mi) f4 



where the quark condensate is given by (f49|). In fact, we recover the Gell-Mann - Oakes - Renner 
(GMOR) relation: 

mlfl = -m < M > . (75) 

Now, we are left with only one free parameter, the constituent quark mass M. In principle, the latter 
can be related to the phenomenological value of the quark condensate (^^I/). We remind that we work 
in an effective non-renormalizable theory with a finite cutoff, so that the quantities which in QCD 
depend on the normalization point should be compared to the particular values of these quantities in 
our model at the scale of order of the model cutoff. Analysis of (Gasser and Leutwyler, 1982) shows 
(M) = (-280 ± 30MeV) 3 at MS renormalization scale /i = 1 GeV. (Gasser and Leutwyler, 1982). 

The dependence of A on M can be seen on fig.|2]. The numerical value of the cutoff A depends 
noticeably on the regularization scheme chosen. The results for the quark condensate (mm) and the 
current quark mass mo can be also seen on the right side of fig.|2] as a function of M. For both 
observables there is a plateau for M larger than 350 MeV. 



2.6 Restriction to Goldstone degrees of freedom: The chiral circle 



Up to now we have considered the meson sector of the model bosonized in terms of independent fields 
a and tt (linear NJL model). However, our main aim is to apply the model to baryons. Although 
it is successful in describing the meson sector (for review see Klevansky 1992, Vogl and Weise, 1991, 
Hatsuda and Kunihiro, 1994), the version with independent a and 7? does not provide a soliton 
solution with non-zero baryon number and therefore cannot be applied to baryons. In fact, it has 
been shown (Sieber et al, 1992, Watabe and Toki, 1992) that this solution collapses to an infinitely 
bound point-like solution in the linear NJL-model . Apparently, the success of the NJL model in the 
meson sector shows that the model incorporates some essential physics and should be not discarded 
but rather modified. In order to understand which ingredient of the NJL model should be revised let 
us recall that we are looking for a model describing the low-energy baryon properties. The particle 
spectrum of the NJL model with independent fields a, n includes quarks, pions and sigma mesons, 
and the sigma meson appears to be much heavier than the others. Furthermore, it is precisely the 
sigma field which drives the collapse of the soliton in the linear NJL model. Thus, we argue that 
the NJL model in its linear formulation does not correctly describe the degrees of freedom related 
to the sigma meson. Since we are interested in the low-energy properties of pions and baryons this 
degree of freedom seems to be not relevant. Therefore, a natural solution to this problem is simply 
to freeze the sigma-degree of freedom. This will be done by invoking the chiral circle condition. The 
corresponding model is called non-linear NJL model which applied to the baryon sector is equivalent 
to the Chiral quark soliton model of Diakonov et al, (1988). 

In order to separate the degrees of freedom corresponding to pions from the sigma meson degrees of 
freedom, let us note that in chiral limit we have a whole family of vacuum solutions connected by the 
chiral transformations ( [231 ) 

a + m ■ f = MU , (76) 
where U is a constant £77(2) flavor matrix, and the meson fields fulfill the chiral circle constraint (01) 

a 2 + tt 2 = M 2 . (77) 

The pions are massless Goldstone particles and in order to construct the effective theory of soft pions 
it is sufficient to consider only slowly changing SU(2) fields U(x) of the form (|76|). With meson fields 



of the form (76) the corresponding semibosonized lagrangian (p2|) can be written in the form 

L = tf[d T + h(U)]V, (78) 

where the current quark mass tuq is included in the one-particle hamiltonian 

h(U) = -nVVfc + M 7 °t/ 75 + m 7° . (79) 

In fact, this treatment of tuq corresponds to a redefinition of auxiliary scalar field a used in the 
identity (ISO). Here we use notation 

^75 = l+Jl u+ [/t. ( 80 ) 

Integrating out the quarks we arrive at the bosonized non-local effective action 

S(U) =-AT c Tr log £>([/), (81) 



where 

D(U) = d T + h{U) . (82) 



From now on we will work with this effective action, i.e. with the non-linear version of NJL model. 



Although, as we will see, the solutions of the non-linear and linear version of the NJL model are 
qualitatively different in the soliton sector, they are equivalent in the meson sector as far as soft pions 
are concerned. This has the important consequence that we do not have to change the way in which 
the parameters of the model are fixed. The pion decay constant is given by the same expression fl73| ) 



and the GMOR relation (|75|) also remains. In fact, the GMOR formula is valid with accuracy 0(mfy 



in any model where chiral symmetry is spontaneously broken. 



As explained above, the lagrangian fl78|) is suitable for the description of soft pions. In the case of 
baryons, however, the pions are no longer "soft" , although their momenta are low (of the order of the 
inverse baryon size) compared to the baryon mass. This leads us to ask: can one use the lagrangian 
(|78| ) for the description of baryons? Apparently the answer to this question goes beyond the model 
itself. It depends on the dynamics of the underlying theory, which in our case is QCD, related to 
the low-energy baryon properties. The physical aspect of the problem is to what extent the chiral 
lagrangian can mimic the "exact" effective low-energy lagrangian of QCD. Indeed, we understand 
that the NJL model is only a model and that the "true" effective lagrangian of QCD has a much 
more complicated structure including non-local multi-fermion interactions. Various attempts (McKay 
and Munczek, 1985, Cahill and Roberts, 1985, Diakonov and Petrov, 1986, Ball, 1989, Schaden et 
al, 1990), have been made to "derive" the lagrangian (|78|) from QCD. In the instanton liquid model 
(Diakonov and Petrov, 1984, 1986) it has been shown that indeed the lagrangian (|78|) can be used to 
describe the pion dynamics relevant for the low-energy baryon structure. 

The instanton liquid model assumes that the QCD vacuum is described by a gluon field configuration 
corresponding to interacting instantons and anti-instantons. The model successfully describes the 
spontaneous breakdown of chiral symmetry which can be understood already on the level of quarks 
interacting with instantons before the reformulation of this interaction in terms of the effective quark- 
pion interaction. Due to the spontaneous breakdown of chiral symmetry quarks acquire a dynamical 
mass M(k). In contrast to the NJL model, however, this dynamical quark mass depends on the quark 
momentum k and the effective lagrangian describing the quark interactions induced by instantons 
is non-local. Therefore, again in contrast to the NJL model, which due to the local four fermion 
interaction is not renormalizable and has to be regularized, the instanton liquid model is ultraviolet 
finite. The typical momentum scale above which the non-locality of the effective quark action becomes 
important is the inverse average size of the instanton 1/p. In particular, the dynamical quark mass 
M(k) can be treated as a constant at momenta < 1/^ and vanishes at momenta k ^> 1/p. Thus 
the scale 1 /p effectively plays a role of a cutoff. 

Besides the average instanton size p there is another important parameter of the instanton model, 
namely the average distance R between instantons. The instanton medium is rather dilute: p/R ~ 
1/3, so that one can develop a systematic expansion in this small parameter. Since the dynamical 
quark mass M(k) appears as a collective effect of the whole instanton medium, it is parametrically 
small for a dilute medium. Indeed, a detailed analysis (Diakonov and Petrov, 1984) shows that for 
p < R 

M(0)p = o(j^) <1 (83) 

which means that the dynamical quark mass M(k) is much less than the cutoff 1/p. On the other 
hand, as mentioned above, M(k) remains constant as long as k is less than this cutoff. Thus the 
presence of the small parameter p/R in the instanton model allows to approximate this non-local 
effective theory by a theory with a local (non-renormalizable) lagrangian and with a cutoff A of the 
order of 1/p. 



A detailed analysis of the non-local four quark interaction of the instanton model leads to the con- 
clusion that in the case of the small parameter p/ R, the pion and quark degrees of freedom can be 
consistently treated within the local lagrangian (|78|) up to momenta k ~ M -C 1/p. The degrees 
of freedom, which lead to the violation of the chiral circle condition, correspond to excitations with 
momentum k of order of Indeed, the specific characteristics of the pion and quark degrees of 
freedom are that the pions are light Goldstone bosons and the quarks have a parametrically small 
dynamical mass M. 

It is important to note that in the lagrangian (|78|) and the effective action (0), the only interactions 
are those of quarks with Goldstone degrees of freedom. These are pions in SU(2) and in addition, 
kaons and //-mesons in SU(3). Furthermore, there is no need to include heavier mesons like e.g. vector 
mesons p, u and A\ in the lagrangian ( fTSf) on the chiral circle. 

Altogether, the instanton model not only naturally leads to the lagrangian (|78|) with meson fields 
restricted to the chiral circle, it also determines the region of the applicability of this lagrangian. 
In particular, it can be used to describe the baryons since the characteristic pion momenta are low 
compared to the cutoff. 



2.7 Generalization to higher flavor groups and relation to the topological approach 

The NJL lagrangian (§) has a structure which suggests a trivial generalization to three and more 
flavors. In such a case, however, the resulting lagrangian is invariant under U(3)l ® U(3)# and some 
terms have to be added in order to break the £7a(1) symmetry (Vogl and Weise, 1991, Klevansky, 
1992, Hatsuda and Kunihiro, 1994). On the other hand, if one works in the instanton liquid model 
(Diakonov and Petrov, 1986) one arrives at a SU(3) effective quark lagrangian of 't Hooft-like form 
which breaks the Ua(1) symmetry explicitly. Further, after bosonization of this 't Hooft-like terms one 
ends up after a couple of approximations with a lagrangian describing the interaction of quarks with 
SU(3) Goldstone bosons. This lagrangian has then exactly the structure (f78[) after simply replacing 
the SU(2) matrix U by a SU(3) one in it. In this way, the lagrangian ([78"D can be immediately 
generalized to any number of flavors. In fact, in section 5 we will eventually apply it to hyperons. 

Although the present model with the meson fields restricted to the SU(Nf) has the same meson 
degrees of freedom as the Skyrme model, in contrast to it our bosonized action (|8l|) is a complicated 
non-local functional. For slowly changing fields U(x) we can perform a gradient expansion for the 
action (|8~ip . From the regularized real part one obtains (in chiral limit mo = 0) 

f 2 r N r 

ReS(U) = ^J d 4 xS P (d^U)(d^) - j d 4 x[2Sp(«9^) 2 + Sp^L^L,)] + . . . , (84) 

where L M = ill^d^U and the ellipsis stands for the terms with higher derivatives. One recognizes the 
kinetic term and the four derivative terms similar to those of the Skyrme lagrangian although the 
latter differ in signs and coefficients. 

The analogy with the Skyrme model gets closer if one considers the imaginary part of the effective 
action. As we already stated before, the imaginary part of the effective action exactly vanishes in 
the SU(2) case. However, for higher flavor groups SU(Nf) the effective action S(U) is generally 
speaking complex. In contrast to the real part, the imaginary part of the effective action is free of 
the ultraviolet divergences. A straightforward way to show this is to perform a gradient expansion 
for this imaginary part. In the leading order of the gradient expansion one finds (Dhar and Wadia, 



1984, Dhar et ai, 1985) the famous Wess-Zumino term (Wess and Zumino, 1971, Wess, 1972) 
lmS(U)\ grad = J d 5 xe mifl4flr Sp [(tfd^u) (u^d, 2 U) (rfd„u) (u*d„u) (tfd» B U 



(85) 



Note that for time-independent fields U(x) the imaginary part ImlogDet.D(£/) vanishes, because the 
one-particle hamiltonian h(U) is hermitian. The same is also valid for the Wess-Zumino term. 

In the Skyrme model the Wess-Zumino term generates the baryon current 

B»(x) = -■^e^Sy(rfd v Urfd a UrfdpU) . (86) 
The corresponding baryon charge 

B = J d 3 xB°(x) = Q T (87) 
is exactly the winding number Qr of the chiral field U(x): 

Qt = ~^2 eijk J ^xSp^diUU^djUU^dkU) . (88) 

Thus, we see that using the gradient expansion one can trace a rather interesting relation between 
the model with the non-local action fl8"T| ) and a local Skyrme-type lagrangian with the Wess-Zumino 
term. However, although this relation is rather attractive from the mathematical point of view, we 
have the problems of the higher order terms as well as with the validity of the gradient expansion 
itself. We remind that the above gradient expansions of both the real and the imaginary parts of 
the effective action are valid only for slowly changing fields U(x). However, as we shall see in the 
next part, in the case of the nucleon the corresponding saddle-point field U(x) varies rather fast, so 
that the gradient expansion generally cannot be used. Therefore, in the present model, the relation 
between the baryon and topological charges used in the Skyrme model is not valid. 



3 Nucleon as a non-topological soliton in SU(2) NJL model 



This chapter is devoted to the description of the nucleon in the NJL model. Using the functional 
integral approach we show that in the large N c approximation the nucleon naturally appears as a 
many-body bound state of quarks in the selfconsistent mean chiral field U. The spectrum of the 
Dirac one-particle hamiltonian with this field U contains a discrete valence level. The nucleon state 
corresponds to the Hartree picture where all negative Dirac sea levels and the valence level are 
occupied by quarks. In fact, a similar physical picture of the nucleon has been first considered by 
(Kahana et al, 1984) in the context of the sigma model. 

Although the mean chiral field U resembles to some extent the soliton of the Skyrme model, the 
important difference is that, whereas in the Skyrme model the baryon number of the nucleon is 
identified with the winding number of the chiral mean field U, in the NJL model the nucleon baryon 
number B — 1 is explicitly given by the presence of the valence level occupied by iV c quarks. In 
order to distinguish the NJL mean-field solution from the chiral soliton of the Skyrme model, it is 
commonly called a non-topological soliton. 

Depending on constituent mass M in the NJL model we have two different and to some extent 
controversial physical pictures of the nucleon: 

• At M ~ 400 MeV a clear valence picture prevails. The valence quarks occupy a bound state in 
the gap between the positive and negative continuum and its single-particle energy is around 
200 MeV. The number of quarks in this valence level determines the baryon number and the 
polarization of the Dirac sea plays the role of the pion cloud. The mean field is noticeably 
varying over the size of the soliton and the gradient expansion of the action does not converge. 

• At M > 1000 MeV the valence level has come asymptotically close to the negative continuum 
and its single-particle energy is slightly above —M. Still the number of quarks in this valence 
level determines the baryon number and the polarization of the Dirac sea plays the role of the 
pion cloud. However, the gradient expansion in this situation works well, and the non-local 
effective action of the NJL model can be well approximated by a local one. Hence, a Skyrme- 
type picture emerges in which the winding number of the meson field equals the baryon number 
of the system. 

• Obviously, since M is a free model parameter, only the agreement with the experiment can give 
a preference to one of these physically different pictures. In fact, all calculations within the 
present model, which reproduce the experimental data, support the picture of valence quarks 
with polarized Dirac sea and not the Skyrme picture. 

Actually, there has been a debate in the literature whether the valence quark picture is still valid if 
vector mesons are explicitly taken into account in the NJL lagrangian. In fact, several theoretical and 
numerical investigations have been performed to clarify the role of p, uo and A\ mesons in the solitonic 
sector (Alkofer and Reinhardt, 1990, Alkofer et al, 1992, Doring et al, 1992b, 1993, 1995, Schiiren 
et al, 1992, 1993, Zuckert et al, 1994, Weigel et al, 1995). Unfortunately, no clear conclusion on the 
formalism and on the effects of the vector mesons on the nucleon observables has been obtained and 
none of the above approaches is beyond criticism. In the present paper, we do not consider vector 
quark couplings or vector meson fields at all. 

The picture of bound valence quarks interacting with a polarized Dirac sea bears some similarity to 
the idea of chiral quark models where the Dirac sea is replaced by the dynamical meson clouds and 



the quarks are confined by a bag (Theberge et al, 1980, Thomas, 1993). In the chiral sigma model 
one does not have confining walls but the system is bound by sigma and pion fields respecting the 
spontaneous breaking of chiral symmetry (Gell-Mann and Levy, 1960, Birse and Banerjee, 1984, 1985, 
Birse 1990). This model allows to use the Peierls-Yoccoz projection method for the quantization of 
collective degrees of freedom (Birse, 1986, Fiolhais et al, 1988, Neuber et al, 1993). There are also 
hybrid models which combine valence quarks with meson clouds. This is e.g. the chiral bag model 
of Brown and Rho (1979). In the model, the pion field outside the bag carries a fractional baryon 
number, and it is connected to the interior by the axial current. The linear chiral sigma model of 
Kahana and Ripka (1984) involves both sea quarks and dynamical pion and sigma fields, and is 
renormalizable. In fact, these authors were the first to treat the polarization of the Dirac sea in the 
context of the valence picture for the baryon (Kahana and Ripka, 1984, Soni, 1987, Li et al., 1988). 
The model shows, however, a vacuum instability and hence has never been applied in a selfconsistent 
way to baryons. 

One should also mention some papers in which the nucleon is not described as a soliton but as a 
three-body bound state as a solution of Faddeev-type equations (Johii et al., 1993, Huang and Tjon, 
1994, Buck et al, 1993, Weiss et al, 1993b, Hellstern and Weiss, 1995). 

In order to describe the nucleon in the NJL model, we consider the nucleon correlation function at 
large euclidean time separation. First, we concentrate on the solution in strict large N c limit - the 
"classical" soliton solution. We discuss the classical part of the soliton mass and the baryon number. 
The selfconsistent hedgehog soliton is presented and the results are discussed. As a second step, 
the quantization of the hedgehog soliton is done. The rotational zero modes are considered and the 
related moment of inertia and the energy of the quantized soliton are given. The translational zero 
modes are also shortly discussed. In the end, we present the evaluation of the nucleon matrix element 
of a quark current in this quantization scheme. 



3.1 Nucleon correlation function 



The correlation functions are a general tool to study the structure of any kind of matter and in 
particular, they provide a powerful method to study the structure of the vacuum and the hadronic 
spectrum in the QCD. The main idea behind this is that the correlation function of a given current 
6 

ll(t)=<0\d(t)d*(0)\0> (89) 
can be written in terms of intermediate physical states 

II(t) =J2\ < 0\6\n> \ 2 e- tEnt , (90) 

n 

which in euclidean space-time (r = it) becomes a sum of decreasing exponents. Hence, from the 
behavior of the correlator at very large euclidean separation r, 

n(r) ~ exp(— mr) , 

one is able to extract the energy (mass) m of the lowest state with quantum numbers of the current 
6. 

In the case of the nucleon, the correlation function is given by 

U N (T) =< 01^(0^/2)4(0, -T/2)\0 > , (91) 



where the iV c -quark current Jn carries the nucleoli quantum numbers JJ3, TT3 and no color (Ioffe, 
1981): 

M% *) = j^j e A '" /,WB T^.TTa <W & *) • • • *Av c / Wc t) ■ (92) 

f f 

Here /3j are color indices and r^j ^ is a matrix with standing for both flavor and spin indices. J 
and T denote the nucleon spin and isospin, respectively. 

The asymptotics of the correlation functions at large euclidean time separations determines the nu- 
cleon mass Mjv 

lim U N (T) ~ e ~ MNT . (93) 

The nucleon correlation function can be written as a path integral with the action (|31|) 
U N (T) =< 0\J N (0,T/2) 4(0, -T/2) 1 >= i J T>&VWDUJ N (0, T/2)4(0, -T/2)e-f d4 * L , (94) 

The effective lagrangian (|78|) is quadratic in the quark fields and a gaussian integration can be 
performed. In this we apply the Wick theorem contracting \fr and iff* from the currents which leads 
to a product of N c quark propagator (diagonal in color indices) in the background of meson field U : 

1 r Nc 1 
n N (T)=T { J } rf*- J VU]J(0,T/2\—-\0,-T/2) fi9i e- s ^. (95) 



3.2 "Classical" soliton 



In order to integrate over the chiral field U, as in chapter 2, we use the large N c argumentation 
and perform the integration in the saddle-point approximation which is exact in large N c limit. In 
this we make use of the fact that in the limit of large time T the saddle-point field configuration is 
time independent which allows us to consider only static fields U below. Here we concentrate on the 
properties of the "classical" soliton which is given by this large N c saddle-point solution. 



Saddle-point solution and the nucleon mass in the leading order of iV c . A specific feature of the nu- 
cleon case is that the functional integral contains now not only the action exponential e~ s ^ with 
S(U) = 0(N C ) but also a product of iV c quark propagators which also depend on the chiral field U. 
Thus, in the leading order of N c these quark propagators also contribute to the saddle-point solution. 



We start with the pre-exponential factor in eq.(|95|). For the quark propagator in the case of the static 
field U we use its spectral representation 

' " ' \X,X 4 ) = e(x' A -X 4 )J2 e ~ en(< ~ X4) ®n(x')®l l (x) 

e„>0 

-0(z 4 - x' A ) e~ eni <- X4) $ n (f) (2) (96) 

en<0 



D(U) 



where $ n and e n are the eigenfunctions and eigenvalues of the one-particle hamiltonian h(U) ( [79]) in 
the static background field U 

h{U)<$> n = e n <f> n . (97) 



In the case x' A — 24 = T > 0, we have 



(0,T/2|— Ly |0,-T/2> = £ e-- T $ n (0)$t(0) (98) 



where the sum goes over all levels of positive energy. At T — ► 00 only the level with the lowest 
positive energy will survive in the sum and hence, for the product of N c quark propagators we have 

N c 1 

n< o . T / 2 i^i°.- r / 2 > T ^ e " m6 ^ (to - (") 

Here, we refer to this lowest energy level as valence level e va i- 

Now, we turn to the fermion determinant. Although this determinant contains ultraviolet divergences, 
we start with the non-regularized version, because it will clarify the physics behind the saddle-point 
approximation. For a static field U following (|81|) we can write 



+00 

da; 



S(U) - S(U ) = -N c Tr\og[D(U)/D(U )} 

= -iV c Tr[log(<9 T + h)- \og(d T + ho)] = -N C T J ^- J d 3 xSp(x | log(ia; + h) - \og(iu + h )\x) 

—00 

= N C T [ +C ° p- [ d 3 x Sp(x \-^— - -^—\ x) = TN C ^ (e a - e°J = TE sea {U) . (100) 
j-00 in J too + n too + h E < 

Here D{Uq) is the Dirac operator corresponding to the vacuum saddle point solution U = 1. The 
vacuum part iV c Tr log.D(£/o) appears due to the normalization l/Z factor in (p5|). h is the Dirac 
hamiltonian corresponding to the vacuum solution. In ( |100| ) we used that Sp(/i — ho) = and 



integrated by parts over u. The energies e a (e°) are the eigenvalues of h (h ) 

We see that S(U) is proportional to the sum of all negative energies, i.e. to the energy of the Dirac 
sea from which the analogous vacuum contribution is subtracted. 

Combining (^9]) and (|100|) we see that 
1 N ° 1 

L rj( , r/2|^y|o, -T/2) e - s ^ ^ j-[^m^ am T (101) 



Now the functional integral over U in (Bo]) can be performed in the saddle point approximation 
justified in the large iV c limit. The saddle-point meson configuration U c can be found from the 
stationary condition 

Su[N c e val (U) + E sea (U)} u = (102) 
and the large T asymptotics of Un{T) in the saddle-point approximation is 

U N {T) ~ e -[Nce va i(U c )+E sea (U c )]T _ ^ 103 ^ 
T — >oo 

In contrast to the vacuum case, the stationary meson field configuration U c is a static localized 
solution which is commonly called a non-topological soliton. Comparing (|TU3) with fl9"H| ) we see that 
in leading order in iV c the nucleon mass is 

M d = min[N c e ml (U) + E sea (U)} = N c e val (U c ) + E sea (U c ) (104) 



This result has a simple physical meaning. In our model the nucleon appears as a bound state of 
quarks in a mean field U c and the energy (mass) of the nucleon state is a sum of one-particle energies 



of all occupied states. These states are eigenstates of the Dirac hamiltonian h(U c ). In order to get a 
colorless state, each level is occupied by N c quarks. Therefore, the contribution of each level to the 
total energy of the nucleon ( |104j ) contains a common factor N c . 

Thus, we see that in the path integral approach, the saddle-point approximation for the correlation 
function of nucleon currents agrees with the Hartree picture of the nucleon as a bound state of quarks 
in the self consistent mean field U c . 

In the above euclidean treatment, the lowest positive level in the spectral representation fl98|) plays a 
privileged role and we called it a valence level. Indeed, as we shall see later, in the case of physically 
relevant parameter values in the NJL model the valence energy is positive, so that the euclidean 
calculation with E sea including all negative levels and with e va [ understood as the lowest positive level 
gives correct result. However, let us assume that under some change of the parameters of the model 
the valence level crosses the zero energy and becomes negative. Formally, in the euclidean treatment 
we face a problem, since in this case the large time asymptotics of the euclidean quark propagator 
(|99|) would be dominated by another lowest positive level. On the other hand, from the physical point 
of view, there is nothing special about crossing zero and the nucleon still remains the bound state 
of the Dirac sea and valence quarks although the valence level is negative. Apparently not the sign 
of the valence level but the way in which the levels are occupied is important. In order to overcome 
this problem, we can generalize the euclidean treatment identifying the valence level with the bound 
occupied level which does not belong to the polarized Dirac sea, i.e. which originates from the upper 
positive continuum. In principle, the valence level can have any sign of the energy. Similarly, the 
Dirac sea energy E sea should be understood as the total energy of all levels which belong to the 
polarized Dirac sea. Actually, the case in which the valence level crosses zero and becomes negative is 
important for understanding the relation between the current approach and the Skyrme-like models. 



3.3 Baryon number of the soliton 

Now we turn to the baryon number of the soliton. We remind that starting from the correlation 
function of two N c quark currents Jjy (|92|), simply by construction, the state with baryon number 
B = 1 is created. The saddle-point treatment described above does not destroy this feature. It is 
easy to see that the physical picture associated with the saddle-point (large N c ) approximation really 
agrees with B = 1. The nucleon in the model appears as a bound state of the Dirac sea and valence 
quarks. Since the sea quarks originate from the lower Dirac continuum, they carry the same baryon 
number B = as the vacuum Dirac sea. Hence, the baryon number B — 1 is given by the N c valence 
quarks, each of them carrying the baryon charge 1/N C , i.e. due to the presence of the valence quarks 
the nucleon has a baryon number B = 1. Explicit proof will be given later when we will evaluate the 
nucleon isoscalar electric form factor. 

We conclude that the quark configuration determines the baryon number of the NJL soliton and 
it is carried by the valence quarks. In particular, in order to consider solutions with a larger, for 
instance B = 2, baryon number one should consider a correlation function of two currents made of 
2N C quarks. It corresponds to a quark configuration in which besides the filled Dirac sea the first two 
levels (coming from upper continuum) are additionally occupied. 

Apparently, this physical picture differs from the one of the Skyrme model where the baryon charge 
is identified with the winding number of the meson field and the field configurations with different 
winding numbers are isolated from each other from the very beginning. The effective action of the NJL 
model is non-local and the transitions between different topological sectors are allowed. Therefore, in 



the NJL model there is no a priori connection between the winding number and the baryon charge. 
On the other hand, in section 2 we stated that in the leading order of the gradient expansion the 
imaginary part of the effective action of the NJL model coincides with the Wess-Zumino term which 
generates the topological current in the Skyrme model. As we will see later, indeed the baryon 
charge, which we discussed above in terms of occupied levels, originates from the imaginary part of 
the effective action. Hence, one expect that if the gradient expansion is valid the baryon charge B 
could be related to the winding number Qt- In order to examine this connection let us consider the 
following quantity 

B(U) = ±-J d 3 x{< tfttf >u-< > Uo } , (105) 

where the brackets mean the euclidean functional integral over the fermion field in a classical back- 
ground chiral field U : 

/ exp {- / d 4 z^[d t + h{Um\ 

<W^f> TT — i . (106) 

U fDVDWexp{-Jd 4 zW[d t + h(U)]y} y } 

Formally B(U) looks like the baryon charge. However, the baryon charge of the soliton depends 
not on the mean field U c itself but on the way in which the quark levels of the hamiltonian h(U c ) 
are occupied, whereas the quantity B(U) depends only on the background field U. Indeed, for the 
time-independent field U 

r 11 

y ' J 1 d T + h(U) d T + h(U ) 

/+00 Am r I 1 
— / d 3 xSp(x I — — ——I x) = Tr W[-h(U)] - 6\-h(U )]](W) 
-00 2tt J PX l iu + h(U) iuj + h{U Q y ' L L v n [ v ;JJV ; 

we see that B(U) counts the number of negative states of the hamiltonian h(U) from which the 
number of negative vacuum quark states is subtracted. We rewrite (|105[) in the form 

^Ctl^^r^k)^' (108) 

where 

h 2 {U) = -d\ + M 2 + iM lk d k U 15 and h 2 (U ) = -d\ + M 2 , (109) 
and perform a local expansion in derivatives dkU/M. The leading term of the gradient expansion 

B{U) = ~^e ihl J d 3 xSp(f/t^f/)(f/t 9fcf/ )( [/ t^ [/) + . . . (110) 

gives exactly the winding number Qt{U) (|88|) of the field U. The higher order corrections in the 
gradient expansion vanish. Indeed, from (|107|) follows that B(U) is an integer number and any higher 
order correction would destroy this feature. However, the equality (|110|) exists only if the gradient 
expansion is valid, i.e. for chiral fields U(x) with characteristic size RM 3> 1. We remind that for 
such fields the non-local action can be also approximated by the first terms of the gradient expansion 
which are similar to those of the Skyrme lagrangian. Thus, we see that the NJL model exhibits a 
certain similarity with local Skyrme-like lagrangians in the case of solitons with MR ^> 1. In this 
limit: 



• the baryon charge coincides with the topological charge (if all and only the negative-energy 
levels are occupied); 

• the non-local action can be approximated by the first terms of the gradient expansion which 
are similar to those of the Skyrme lagrangian. 



Regularization. As we already know from the vacuum sector the fermion determinant is divergent 



and must be regularized. We start from the proper time regularization formula (|4T|). In the case of 
static fields U(x) it gives 

ilr log (D^D) reg - ilr log {D\D,) reg = ~T j°° ^ f°° ^Tr \e^ 2 ^ - e^^] , (111) 

where Dq = D(Uq) corresponds to the vacuum. The regularized analog of the Dirac sea energy E sea 
(|100|) is given by 

EL(U) = ^Tr log (rfD) reg - ^Tr log (dJA,)™, . (112) 



We perform the integral over lo in (111) and express the trace through the sum over the eigenstates 
of the Dirac operators h(ho). Finally, we get 

EL = N c J: Rf(e n ) - A c £ Jtf (£) , (113) 

n n 

where the regularization function R±(e n ) is given by 

In the limit of the large cutoff A — ► oo the regularized sea energy E^ ea reduces to the non-regularized 
sum over negative energies (|100| ). If the valence level is positive its energy should be added to this 
sum but if the valence level energy is negative its contribution is already contained in E^ ea , so that 
the total regularized energy of the soliton is 

M d = N c e val 9{t val ) + Ei a . (115) 



3.4 Stationary hedgehog meson field configuration. 



Even in the simplest SU(2) meson field configuration U (ffq) , in which the meson fields are constrained 
on the chiral circle (01), one needs additional assumptions in order to solve eq. ( |102| ) numerically. 
Usually for the stationary meson field configuration a hedgehog structure is used 

U c (r) = e^M^) . (116) 

The hedgehog is the most symmetric meson field configuration in which 

vr a (f) = n a n(r) (117) 

and the classical meson fields are parametrized in terms of the profile function of the soliton @(r) as 
follows 

cr(r) = M cos 0(r) and n(r) = Msin0(r) . (118) 
The chiral profile function satisfies the boundary condition 

9(r) ~ 0, (119) 



which ensures that the physical vacuum U = 1 is recovered at r — > oo. In the Skyrme model, 
additionally to condition flH9p , it is imposed that 

9(r = 0) = -vrQ T , (120) 



where the number Qt (winding number) is identified with the baryon number and as such it should 
be integer. In the present NJL model the Qt is not necessarily an integer number, since due to 
the non-locality of the effective action we could have in principle singular saddle-point solutions. In 
addition, even an integer Qt in the NJL model is not necessary identical to the baryon number. 



At first glance the assumption of the hedgehog structure ( |116| ) appears arbitrarily in the NJL model 
and, indeed, there is no general proof that it corresponds to a stationary solution of ( |102| ). However, 
such a proof exists in the Gell-Mann and Levy sigma model with valence quarks (Ruiz Arriola et al, 
1989). Additionally, there is a general argument in favor of the hedgehog choice. In particular, one 
expects that the classical solution should have the most symmetric structure which is the hedgehog 
one. A deviation from the hedgehog configuration would lead to a rather complicated spectrum of 
rotational excitations of the nucleon which is difficult to be classified and in fact it is not observed in 
nature. 



For the hedgehog meson fields the stationary condition ( |102| ) can be written in terms of the profile 
function Q(r) 



0. 



;m) 



e=e c 



5e[N c e val (Q) + E sea (Q)\ 

In particular, it leads to the following classical equations of motion written in terms of @(r) 

sin e{r)S{r) -cos 6(r)P(r) = 0, (122) 
where the scalar S(r) and pseudoscalar P(r) densities read 



S{r) = N C M 



]Ti^(e n )$ n (r> n (r) + 9(e val )$ val (r)$ 



vai(r 



P(r) = N C M 



^^(e n )$ n (r)^ 75 (r a n a )$ n (f) + 6(e val )^ val (r)i^(r a h a )^ 



vai(r 



(123) 



(124) 



respectively. The regularization function R2 (e n ) is given by (Reinhardt and Wiinsch, 1988, Meissner 
Th. et al, 1988) 

1 r°° du „.,2 



R2 ( e n) 



e„.e 
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(125) 



One-particle hamiltonian h and Dirac spectrum. Due to the hedgehog ansatz ( |116|) , the one-particle 
hamiltonian h ([79]) is invariant only with respect to simultaneous space-isospin rotations. It means 
that it commutes with the grand spin G = j + t, but not separately with the total momentum 
j = 1 + s and the isospin t of the one-particle eigenstates of h. Since it also commutes with the parity 
operator IT, the hamiltonian h, the grand spin G, its ^-component G z and the parity operator II form 
a complete set of commuting operators. Therefore, the eigenstates $ n of h: 

h$ n {r) = e n $ n (r) (126) 

can be characterized by the energies e n and three quantum numbers, which correspond to G, G z and 

n. 



Following Kahana and Ripka (1984) the eigenvalue problem (|126|) can be solved numerically in a 
finite quasi-discrete basis. The basis is made discrete by putting system in a large box of radius D 
and imposing boundary conditions at r = D. Also, it is made finite by restricting momenta of the 




Figure 3: Quark spectrum of the one-particle hamiltonian h for an exponential profile form 0(r) = 
— 7rexp(— r/R) as a function of the profile size R. Both the energies e n and R are given in scaled 
units: (ff and MR). 



basis states to be smaller than the numerical cutoff K max . Both quantities have no physical meaning 
and the results should not depend on them. In particular, K max has nothing to do with the model 



cutoff A (|111|) . The typical values used are D ~ 20/M and K max ~ 7M. 



As an illustration a typical spectrum of h for a simple meson profile function B(r) = — n exp(— r/R) 
are shown in Fig. |3|. This G is not a selfconsistent solution and it is purely chosen for didactic 
purposes. The parameter R is the size of the meson profile. Together with the constituent quark 
mass M, which describes the strength coupling between quarks and mesons (or equally the magnitude 
of the chiral field), the product RM characterizes the deviation of the actual field configuration from 
the vacuum (MR = 0) one. One can vary MR either changing M at a fixed profile size R or vice 
verse. In most of our discussions, we prefer to use a language in which M is kept fixed. 

At some R the lowest quark level G n = + starts to deviate from the positive continuum and at large 
enough (R > 0.5/M) one finds a bound level which we call the valence level | n) = | val). Its single- 
particle energy e va i decreases with increasing R, changes sign and at very large R (RM 3> 1) this 
level approaches the negative continuum (Dirac sea). Thus, for large R the valence level approaches 
the negative continuum. In this cases the chiral field is slow varying (RM ^> 1) and the equality 
( |110p is valid, i.e. the baryon number can be related to the topological winding number. Similar 



considerations hold for higher winding numbers (Kahana et al, 1984). At sufficiently large RM ^> 1 
one therefore gets close to the picture of the topological soliton models like e.g. the Skyrme model 
which contains no quark degree of freedom and the winding number of the chiral field is related to 
the baryon number of the soliton. At moderate values of R (RM « 1), however, the NJL model 
supports a different physical (i.e. valence) picture: there is a positive valence level which originates 
from the positive continuum, and it gets bound due to the interaction of the valence quarks with the 
polarized negative continuum (Dirac sea). This mechanism creates the non-topological soliton as a 
bound many-quark state. In order to illustrate it, we plot the soliton energy as well as the separated 
valence and Dirac sea parts for the same meson profile as a function of the profile size R in fig.f|. One 
clearly recognizes that M c i has a local minimum at R M « 1 which corresponds to a solitonic solution 
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Figure 4: Classical soliton energy M c i split in Dirac sea and valence parts in the case of an exponential 
profile form. The energies are given in unit of the constituent quark mass M and the size of the soliton 
R is given in terms of the scaled unit MR. 



of B = 1. This minimum is due to the competing behavior of the valence energy E va i and of the sea 
contribution E sea . At very small RM < 1 we have actually the vacuum quark spectrum with both 
unpolarized continua which means that E sea = (no polarization). With increasing R due to the 
interaction between the valence quarks and the Dirac sea, the latter gets more and more polarized 
and the polarization energy E sea increases. The valence level becomes more bound and the valence 
energy decreases. Thus, at small R the valence part dominates over the sea one and it changes to 
the opposite at large R. The B = 1 soliton solution which minimizes the soliton energy occurs at 
intermediate RM « 1 and the soliton consists of a positive valence level (well separated from the 
Dirac sea) and a polarized Dirac sea, and both parts contribute almost equally to the soliton energy. 



3.5 Selfconsistent soliton solution 



After some pioneering calculations with fixed meson profile (Diakonov et al., 1988, Meissner Th. 
et al., 1988) the selfconsistent soliton solution with baryon number one has been first obtained by 
Reinhardt and Wiinsch (1988), and by Meissner Th. et al. (1989). Using a selfconsistent procedure, 
similar to ones known from Hartree and Hartree-Fock mean-field calculations in atomic and nuclear 
physics, ( p.22| ) and ( |126| ) have been solved numerically. One starts with a reasonably chosen profile 
0(r), diagonalizes the h in the quasi-discrete finite basis (Kahana and Ripka, 1984) and using the 
eigenfunctions an d the eigenvalues e n obtains a new profile function G'(r) from the equation 

of motion (|122|). The procedure is repeated until a desired degree of selfconsistency is reached. 



Typical results for the energy M c \ of the selfconsistent B = 1 soliton and the separated valence 
and sea contributions are presented in fig.[5| as a function of the constituent quark mass M. The 
solitonic solutions exist if the M exceeds a critical value M > M cr « 350 MeV. This cusp behavior 
is typical for localized (soliton) solutions of a system of coupled non-linear equations (see e.g. Lee, 
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Figure 5: Classical energy M d of the selfconsistent soliton split in sea and vacuum contributions as 
a function of the constituent quark mass M. 

1981 and refs. therein) and has been observed also in other chiral quark meson models (Birse, 1990 
and refs. therein). The soliton energy shows a rather weak dependence on M as it slightly decreases 
with M increasing. However, both the valence and the sea parts depend strongly on the mass M: 
with increasing M the valence energy decreases and correspondingly the sea one increases. It is a 
behavior, similar to the one obtained with the fixed profile and variation of R, and apparently, this 
is an indication for a change in the soliton structure. Below M cr there is no soliton solution and 
the N c free quarks is the most preferable configuration. At M > M cr there is a soliton which is 
characterized by a bound valence orbit and polarized Dirac sea. With increasing M the valence level 
goes down and at very large values M ~ 1 GeV approaches the negative continuum. In the latter 
case, since the mass M is large enough (MR ^> 1) and all negative-energy states are occupied, the 
baryon number can be identified via gradient expansion with the winding number of the chiral field 
as already discussed in the context of the fixed profile. Obviously, B(U) has a non-zero value and it 
is equal to the baryon number B — 1. Thus, the model allows for two different and to some extent 
opposite nucleon scenarios. At M < 600 MeV the soliton consists of a well pronounced bound valence 
level and a polarized Dirac sea which corresponds to the widely accepted phenomenological valence 
picture of the nucleon as valence quarks surrounded by the pion cloud. At much larger values M > 1 
GeV the valence level can be considered as a part of the Dirac sea and we are close to the Skyrme 
picture of the nucleon as a topological soliton. In order to make a definite decision which picture is 
favored one has to calculate the baryonic observables using M as a free model parameter to reproduce 
as many as possible experiment data. As we will see later, fixing the parameters of the model in the 
vacuum to reproduce the physical values of the pion decay constant and the pion mass, the model is 
able to provide an overall good description of the baryon properties for M between 400 and 450 MeV. 
Higher M values are excluded. It means that in the case of phenomenologically acceptable values of 
the parameters of the NJL lagrangian the saddle point field U c does not fulfill the condition MR 3> 1. 
Hence, the gradient expansion cannot be used, and the approximation of the non-local NJL effective 
action by a local one similar to those of the topological Skyrme-like models is not justified. 



On fig.^l one can see the profile function of the selfconsistent B = 1 soliton solutions for three values 
of M as well as the corresponding mean-field densities (saddle-point meson fields) S(r) and P(r). As 
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Figure 6: Selfconsistent profile function (left) and the selfconsistent scalar 5* and pseudoscalar P 
mean fields (right) for three different values of the constituent quark mass M and = 140 MeV. 



one can see the actual form of the profile function and the meson fields as function of the distance 
r are rather independent of M. The typical size of the soliton is about 0.6 fm which means that 
RM fa 1. The profile function fulfills the boundary conditions (|119| ) and (|120|) . There is, however, a 
priori no reason in this model to fix Qt to one because the baryon number B — 1 is carried by 
the valence quarks. This is different from the topological soliton models (e.g. the Skyrme model), 
where the baryon number is identified with the winding number from the very beginning. There is 
even no a priori reason why Qt should be an integer number, if one leaves mathematical arguments 
like continuity of the tt field in the origin. In contrast to the Skyrme model, where an integer Qt is 
required in order to obtain a finite energy, the total regularized mass M c \ is finite for any value of 
Qt in the present approach. Apparently, we cannot expect that the saddle-point field U c describing 
the baryon with the baryon charge B = 1 should have the same topological charge Qt{U c ) = B. 
Looking for the solution of the stationary equation ( |102| ) we have to minimize the functional ( |104| ) 
over all fields U which could have in principle different topological charges. However, the explicit 
calculations show that the saddle-point field for the B = 1 soliton indeed has the winding number 
Qt = 1. Diakonov et al. (1988,1989) and Berg et al. (1992) have investigated the behavior of the 
energy M c i(Qt) of the B = 1 soliton as a function of Qt for a fixed profile form. It turned out, that 
the minimum of the soliton energy M c i definitely lies at the point Qt = 1, which means that this 
feature at least for the case B = 1 seems to be a dynamical consequence of the equations of motion. 
Unfortunately the whole analysis is based purely on numerical arguments and up to now there is no 
general proof of this fact. 



In addition Berg et al. (1992) also studied selfconsistent solutions with higher baryon numbers B > 1. 
They found that a local minimum which corresponds to a solitonic solution only exists if Qt < B and 
in any case these local minima appear for integer winding number Qt- For example, for the B = 2 
system they found a soliton in the Qt = 1 sector with 2 valence orbitals (0 + ,0~) occupied if M > 370 
MeV. In the Qt = 2 sector, solitonic solutions exist if M > 560 MeV. The main difference between 
the two cases is that for Qt = 2 and large mass M both the + and the 0~ orbital cross the zero line 
and approach the Dirac sea, whereas in the case of Qt = 1 only the + orbital comes down and the 
0~ orbital moves back to the positive continuum. 



Sieber et al. (1992) performed a detailed analysis of the NJL model without the chiral circle condition 
(U) and showed that in contrast to claims of Meissner Th. et al. (1989) and Reinhardt and Wunsch, 
(1989), at least in case of the proper time regularization scheme, no localized solitonic solution with a 



finite energy exists. Using a modified version of the NJL model, which in addition to the spontaneously 
broken chiral symmetry also simulates the anomalous breaking of scale invariance in QCD, Meissner 
Th. et al. (1993b) and Weiss et al. (1993a) have shown that in this case stable solitonic solutions 
exist without the non-linear constraint which could be considered also as a physical motivation for 
the chiral circle condition (|j) in the present model. 



3.6 Semiclassical quantization of the soliton 

In the previous section we have shown that in the leading order of the large N c approximation the 
calculation of the nucleon mass reduces to solving the saddle point equation (|104|) . Apparently, if we 
want to compute the 1/N C corrections to it we have to take into account the quantum fluctuations of 
the chiral field U(x) around this classical saddle point solution U c (x). Among all possible fluctuations 
a special role is played by those which do not change the classical soliton energy. Since the hedgehog 
soliton is not invariant under translations, space rotations and isotopic transformations (although it is 
invariant under the simultaneous space and isotopic rotations), applying these symmetry transforma- 
tions to the saddle point solution we can construct a whole family of new solutions which correspond 
to one and the same classical energy M c i. We remind that we had an analogous situation in the 
vacuum sector where in the chiral limit there also was a family of vacuum solutions ([76]) with a spon- 
taneous breakdown of the chiral symmetry. However, since all vacuum solutions reflect the fact that 
the chiral symmetry is spontaneously broken, in the vacuum sector we can choose one of the vacuum 
solutions and forget about the others. In the soliton sector the classical hedgehog solution cannot 
be associated directly with the nucleon physical state. Indeed, the physical nucleon state should 
have proper quantum numbers like spin, isospin momentum. Such quantum states can be obtained 
as quantum superposition of classical states with different space and isotopic orientations and with 
different space locations. The corresponding techniques are well known in the non-relativistic many- 
body physics (see Ring and Shuck, 1980) and they are also used in the field theory. In particular, 
the cranking (semiclassical quantization) scheme is successfully applied for the quantization of the 
Skyrmion soliton (Adkins et al, 1983) as well as in the chiral sigma model with valence quarks (Cohen 
and Broniowski 1986). However, in the present model we have a more complicated case because of 
both the non-locality of the effective action (fermion determinant) and the ultraviolet divergences 
which require regularization. 

The symmetry of the hedgehog solution leads to the equivalence of the space-rotational and of the 
isotopic zero modes, so that we are left with three rotational and three translational zero modes. In 
order to restore the rotational, isotopic and translational symmetries, we have to treat the degrees 
of freedom related to these symmetries on the quantum level. Even if we treat the meson fields 
in the leading order of the 1/N C approximation (no meson loops), in order to obtain states with 
good quantum numbers, we have to include the zero modes in our considerations. In contrast to 
other normal modes the zero-mode fluctuations are not small. Since they cannot be treated in a 
perturbative way, one should integrate over these fluctuations exactly in terms of path integrals. 
In fact, the semiclassical quantization scheme for the soliton of the NJL model, including polarized 
Dirac sea, is elaborated in (Diakonov et al, 1988). The problems arising in this scheme due to the 
regularization are considered at first in (Reinhardt, 1989). 



Rotational zero modes and the quantum rotational corrections to the soliton energy. In this part we 



restrict ourselves to consider only the rotational zero modes needed to construct the nucleon states 
with certain spin and isospin eigenvalues. To this end we restrict the path integral in (|9"4j) to the 



meson fields of the form: 

U(r,x) = R(r)U c (x)R^r), (127) 

where R(t) is a unitary time-dependent SU(2) orientation matrix of the soliton and U c is the sta- 
tionary meson field configuration. Now we are left with a path integral over time dependent matrices 
R 

r Nc 1 
U N (T) = T { J } Tf* J VRl[(0,T /2\ \0, -T/2) hgi e~^ RU ^ . 



For ansatz 



~(RU C W) 
can be rewritten as 



:i28i 



the operator D(U) 

D(U)=R [D(U C ) + R j R] R ] = R [D(U C ) + ifi] R ] , 



(129) 



where the dot stands for the derivative with respect to the euclidean time r. In (|129|) we introduced 
a hermitian angular velocity matrix 



Q = -iR ] R = -^Sp(i? 1 Rr a )T a = ^Q a r a . 



With Dirac operator (|129|) the effective action reads (before regularization) 

S{RU C R ] ) = -N c Trlog[D(U c ) + ifi] . 



(130) 



(131) 



In fact, the D(U C ) + iQ corresponds to the body-fixed frame of the soliton in which the quark fields 
are transformed as 

— ► R(t)V and ¥ — ► &tf(t). (132) 
Similarly, the quark propagator in the background meson field U can be presented in the form 

, , 1 , 1 
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D(U) 



x') f9 = (x\R(x 4 ) 



D(U C ) + in 



R\x' A )\x') 



fg ■ 



(133) 



In order to perform the integration over R in ( 128|) we treat the angular velocity R^R of the soliton 
as a small quantity and neglect its derivatives. In fact, as we will see later this is justified by 1/N C 
arguments. Therefore, in the local expansion in powers of the angular velocity and its derivatives we 
can restrict to the terms 



N c Trlog[D(U c )+iQ] = iy c Trlog[,D(t/ c )] - ^ J drQ 2 a 



(134) 



The coefficient I sea is the Dirac sea contribution to the moment of inertia of the soliton. If we forget 
about the regularization, after some straightforward calculations we get 
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(135) 
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Here, similarly to the baryon number the trace is calculated explicitly in terms of quark matrix 
elements. Since the sea contribution to the moment of inertia originates from the effective action, it 



is divergent and should be regularized. It is done (Reinhardt, 1989) using the proper-time regularized 
expression for the real part (|4"1~D and the expansion 



1 

e A+ *= e A + [dae aA Be^ A 
o 

1 1-/3 

+ J dp J d«e ai Be^5e (1 - t, -" )i + ... (136) 



o o 

The final regularized expression for the soliton moment of inertia reads 



IseaS ab = l -N c £ Kf{e m , e n ){n \r a \ m) (m \r b \ n) , (137) 



where the proper-time regulator is given by (Reinhardt, 1989) 



1 °f du ( 1 e ue " — e ue ™ e n e + e m e l " 7 



^/(^,e n ) = _ / ^ _ - -» . (138) 

For the quark propagator using its spectral representation and treating Q as a perturbation we obtain 

^'^l''-*^ 11 ' (-N^T-'-fJtml) . (139) 
Analogously, J m / is the valence quark contribution to the moment of inertia 

-r r 1 »r \— » (val\r a \n)(n\T b \val) 

Ivarfab = ~N C Y, ~ _ " • (140) 

1 n^val 6 ™ €val 

It originates from the quark propagators and is not affected by the regularization. The total moment 
of inertia of the soliton is now 

I = I val + I sea ~ N c . (141) 
Combining ( |140| ) with the non-regularized I sea ( |135| ) 

•£ < " |Tl '"i ( '" |Tt| " > . (142) 

we recover the Inglis formula for the moment of inertia in the cranking approximation known from 
the many-body theory (Ring and Schuck, 1980). 

Now, using expansions ( |134| ) and ( |139| ) we arrive at a functional integral over time-dependent orien- 
tation matrices R(t) with an action quadratic in Q. The quadratic action is in fact the action of a 
symmetric rotator and hence, the path integral reduces to the euclidean evolution matrix element for 
the corresponding quantum rotator: 



R(T/2)=R 2 

T>Rexp 



( r T l 2 



' 1 drill I = (R* \e- HrotT \ Ri) . (143) 



, 2 .. 

R{-T/2)=R 1 \ -T/2 

Here we fixed the boundary conditions in the path integral at r = ±T/2 postponing the integral over 
_R(±T/2) for later. The quantum hamiltonian H rot acts in the Hilbert space of functions ip(R) where 
R is the SU(2) orientation matrix of the soliton and in the rhs of (|143 ) we have the evolution matrix 



element in this R representation. From quantum mechanics we know that the hamiltonian of the 
symmetric rotator can be expressed through the squared angular momentum: 

j 2 T 2 

H ro t — — > (144) 
21 21 ' v ; 

where the action of generators J a and T a on the wave function ip(R) is as follows 

exjp(iu a T a )ip(R) = 4) (exp(-iu a r a /2)R) , (145) 

exp(iuj a J a )ip(R) = ip {Rexp(iu; a T a /2)) . (146) 

Comparing these formulas with the form of the rotating soliton RU C R^ we see that the left generator 
T a is nothing else but the isospin operator, whereas J a is the spin one. We stress that the identity 
(|143|) corresponds to the canonical quantization rule (in the euclidean time) 



- -ij (147) 

The eigenfunctions of the operators J 2 = T 2 , J 3 , T 3 can be expressed through the Wigner D-functions 
as follows 

ijjj 3 TT 3 (R) = (-1) t+T3 VW+TD% j J3 (R) . (148) 
Obviously these functions are also eigenstates of the hamiltonian H rot ( |144| ) 



Hrot^JJsTT 3 (R) = J(J 2 j 1 V jJ 3 TT 3 (R) • (149) 



The path integral over the orientation matrices R with the local action quadratic in f2 is computed 
exactly in (|143|) although this local quadratic form of the action is an approximation to the non- 



local rotational action ( 134 ) justified in the large N c limit. Indeed, we are interested in nucleon 
rotational excitations with J = T of order of 0(N®). For such excitations the angular velocity of the 
semiclassical trajectories, which saturate our path integral in the large N c limit, is Q a ~ 1/1 ~ 1/N C 

The above argumentation that the left and right generators T a and J a can be interpreted as spin and 
isospin operators is based rather on physical intuition than on a consistent derivation. In fact, the 
spin and isospin of the nucleon state created by current J N are determined by the tensor on the 



rhs of fl9~2|) . Let us show how the quantum numbers T 2 = J 2 , T 3 , J 3 of the current J^ N penetrate into 
the rotational wave functions of the nucleon. We can rewrite the quark propagator in the background 
field of the rotating soliton RU C R) in the form 

{ °^D(MW) l0 -- T/2) = W2)<0,T/2| 5j ^ Tli j|0,-r/2>flt(-r/2). (150) 
Inserting this into ( |128| ) we arrive at the following product 



E Nc IlR^iT/2) . .. Rt , N , tNo (T/2) ^rl^'^^'^P^)] • (151) 

Here we have taken into account that tensor T]^ contains Clebsh-Gordan coefficients which relate 
the isospins of separate quarks to the isospin T of the nucleon, so that the isospin rotation of 
N c quarks is equivalent to the corresponding isospin rotation of the nucleon isospin index T 3 . Next, 

r ft 

tensor T]^ also contains spin indices of separate quarks that are contracted with spin indices of the 
quark propagators ( |150| ). At large euclidean time separation T, the quark propagator is saturated by 



the valence quark contribution (|9"9"|), so that effectively the spin and isospin indices of rffi become 
contracted with the indices of the N c valence wave functions. Using the fact that the valence level 
has the grand spin G = one can see that after some manipulations the valence wave functions 
contracted with the function D^ T , (i?) in (|151|) result in 

(-l) J+J3 ^ 3 ,-^( T /2)] = (-l) T+T3 ^ 3 ,.3^( T /2)] (152) 

which up to a normalization factor is identical with the rotational wave function of the final state 
-0}j 3TT3 [i?(T/2)] ( |148| ). In a similar way, one obtains the wave function ipjj 3TT3 [R(T/2)] of the initial 
nucleon. Note that in eq. (|143|) we fixed the boundary conditions for R(—T/2) = Rx and R(T/2) = 



i?2- Now we can turn back to ( |128| ) and perform the integral over Ri, R2. We get a rather simple 
result 

J VRr Jj3 TTM T / 2 )] ^jj 3 ttM- T / 2 )] ex P {-[ I d ^ 2 } 
= JdRr jj3 TT 3 (R)^jj 3 TT 3 (R)e-^ T =(JJ3Tn\e- H ^ (153) 

From the asymptotics of the nucleon correlation function 



ITv(T) ~ exp 

1 — >oo 



(154) 



we obtain the quantized soliton energy 



Ej = M d + ^pl, (155) 

which contains the familiar rotational corrections. In this formula the classical contribution is M c \ = 
0(N C ) and the rotational correction is of order 1/1 = OlN^ 1 ). On the other hand, there exist 0(N®) 
corrections coming from the non-zero modes which were neglected in our derivation. These non-zero 
mode corrections are parametrically larger than the 0(A^ C _1 ) rotational corrections. However the 
non-zero mode corrections do not depend on the spin and isospin of the nucleon state so that the 
non-zero mode corrections cancel in mass splittings. In particular, in the semiclassical quantization 
scheme the nucleon-delta mass difference is simply given by the rotational corrections 

3 1 

E NA = E 3 / 2 — Ex/2 = — ~ — • (156) 

We get the expected results that in large N c limit the nucleon and Delta are degenerate in mass. 

The moment of inertia has been calculated numerically by Wakamatsu and Yoshiki, (1991) and 
Goeke et al. (1991). The results are shown in fig. |7] as a function of the constituent quark mass M. 
One notices a clear decrease of the moment of inertia with M increasing. At M w 420 MeV the 



experimental nucleon-Delta mass difference of 294 MeV can be reproduced (see eq. ( |156| ). In this 
region, the valence quarks dominate and the sea contribution to / is about 20 %. The N — A mass 
difference increase with M and as can be seen the values for M w 420 MeV are favored. 



Translational zero modes. So far we have considered only the rotational zero modes. However, the 
hedgehog soliton U c breaks not only the rotational and but also the translational symmetries. In 
order to take into account the translational zero modes, we consider the "moving soliton" extending 
the path integral over all time-dependent meson fields of the form: 

U(x) = U c (x - X(r)) = T x{t) U c {x)T ] x[t) . (157) 
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Figure 7: Moment of inertia and its valence (dashed) and sea (dash-dotted) parts (left), and the 
nucleon-Delta mass splitting (right) for the selfconsistent solution and m n = 140 MeV in dependence 
of the constituent quark mass M. The experimental nucleon-Delta mass difference is also shown. 



Here T x is a unitary translation operator which corresponds to a translation X. Analogously to ( |129| ) 
it enters the Dirac operator 

D(U) = T x \D{U C ) + T x f x ] T x = T x [D{U C ) - it ■ V] T x , (158) 

where dot stands for the derivative with respect to the euclidean time r. Further, the effective action 
can be transformed to 

S(T X U C T X ) = -Tr log [D(U C ) - X k d k ] . (159) 

and expanded in X k 

N c Ti log [D(U C ) - X k d k ] = iV c Tr log [D(U C )} - ^ J drX 2 + ..., (160) 

For the quark propagator the expansion leads to 

(0,T/2\T X —. T X \0,-T/2) N < ~ ^n^t-M^l J drX ^ (m) 

D(U C ) - X k d k t^oc 

It can be shown (Pobylitsa et ai, 1992) that 

M sol = M val + M sea ~ N c . (162) 

coincides with the soliton mass M c [. The integral over X with an action quadratic X can be easily 
computed and now we arrive at the following expression for the nucleon mass 

Ej = M d + J(J+1) + . (163) 

As it is expected the large N c limit leads to a non-relativistic approximation. 

One realizes that although the model is able to reproduce the N — A mass difference, the masses itself 
are by 25% too large. The reason is that the hedgehog soliton as a mean-field solution breaks the 
translational and rotational symmetry and because of that the classical mass M c \ contains spurious 
contributions of the center-of-mass motion which should be subtracted. The problem is well-known 



in the many-body theory (Ring and Schuck, 1980) and analogously to so called RPA scheme one 
expects two main subtraction terms, both of order of 0(N®): 



AE = -^-±£±, (164) 
2M sol 21 ' K J 

where p and j are the linear and angular momentum operators, and M so i and I are the mass and the 
moment of inertia of the soliton . A rather crude estimate of the first term can be done in an simple 
oscillator potential model. An estimate of the spurious contributions using directly the formula ( [L64j ) 
and calculating < p 2 > and < j 2 > for the soliton reduces the nucleon mass by more than 30 % 
(Pobylitsa et al, 1992). In the Skyrme model this question has been investigated quite in detail by 
including pion loop effects (Moussallam, 1993, Holzwarth, 1994, Holzwarth and Walliser, 1995). The 
final value for the nucleon mass comes out close to experiment. 

Despite the problem of the spurious contributions to the soliton energy and the related overestimation 
of the masses has attracted a lot of attention (Pobylitsa et al, 1992, Weigel et al, 1995b, Broniowski 
et al, 1995) a strict derivation and a reliable quantitative estimate of those contributions in the 
present model is still missing. The reason is that in order to solve this problem rigorously one has to 
include the meson fluctuations. In the path integral formalism it means that we have to go beyond 
the present large N c saddle-point treatment of the path integral (classical meson fields). Note that 
although using the semiclassical quantization scheme we are able to assign proper spin and isospin 



number to the soliton, now the ansatz ( 127 ) is not enough. The quantum meson contributions to the 
effective action, which we neglect as parametrically small (of order 0(N C ) and hence suppressed by 
l/iV c ) in the present considerations, have to be taken into account. They are given by the third term 
in (|56"D |Tr log Jr|r . However, including this term one encounters the following problems: 



• The functional trace over the meson fields includes new divergences and one needs a new cutoff 
in order to make the theory finite. It should be done in a consistent way for both the vacuum 
and the soliton. 

• The meson-loop term should be evaluated directly in euclidean space-time, since the used 
proper-time regularization makes the analytical continuation to Minkowski space-time ill de- 
fined (Broniowski et al, 1995). The usual dispersion relations for the Green's function are not 
valid and hence, one cannot use the convenient particle-hole description as it is nevertheless 
done in (Weigel et al, 1995b). 



3.7 Nucleon matrix elements of quark currents 



Most of the nucleon observables are related to the evaluation of the nucleon matrix element (N | (p')j \ N(p)) 
of various quark currents j. In this subsection, we will demonstrate how this matrix element can be 
written in terms of path integrals and can be evaluated by using the 1/N C expansion in the semiclas- 
sical quantization scheme. 

We start with a matrix element of an arbitrary quark current ^O^/, where O is some matrix with 
spin and isospin indices. Analogously to the nucleon correlation function (p4]), the nucleon matrix 
element can be expressed as an euclidean functional integral: 

(N'(j?)\&(0)dV(0)\N(p)) = 4 / d 3 xdVe-^'+^ 

jvU Jvm J J N ,tf ,T/2)^ t (0)6^(0)4(x, -T/2)e~ J '*** D V)* . (165) 
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Figure 8: Diagrams corresponding to the expansion in Q of the current matrix element: a) the valence 
contribution and b) the Dirac sea contribution. 



Here the nucleon state \N) is created by the nucleon current j\q (j92|) . The constant Z' in eq. ( |165| ) 
is chosen, so that the nucleon states obey the non-relativistic normalization condition: 



(N(p)\N(p)) = (2ti) 3 5^{0 -p). 



(166) 



Integrating out the quarks in (|165|) the result is naturally split in valence and Dirac sea parts (see 
Fig# 

(N'(^)\¥d^\N(p)) = (N'(fi)\¥6y\N(p)} sea + (N'(^)\¥d^\N(p)) val , (167) 



where 



and 



(N\fi)\¥6^\N(p)) sea = N c -^- f d 3 xd 3 i'e~'^' +i ^ 

T— >+oo Z J 



r Nc 1 —1 

J vu r^rf* n W> T / 2 I n(u) ^ - T / 2 )t*&>{°*> ^ 1 DfUjM**} 

(N'(fi)\^6y(0)\N(p)) val t = N c — J d 3 xd 3 x'e- i[fs ' +i ^ 

r Nc 1 

x J Vl/rtfT^* jl(*,T/2\^\x-T/2) fkgh 

x(a;\T/2|^|0,0> /ld d^ 



1 JV c Trlog[£)(C/) ( 



(168) 



(169) 



As in the case of the nucleon correlation function, the integral over the meson fields can be done in 
the saddle-point approximation. Since the hedgehog stationary solution U c breaks the rotational and 
translational symmetries and we need nucleon states of both good spin and isospin numbers and a 



well defined linear momentum, we have to consider both the rotational and the translational zero 
modes. To this end we extend the path integral over all time-dependent meson fields of the form: 

U(x) = R(t)U c (x - X(t))R?(t) , (170) 

and follow the same steps as in the case of the nucleon correlation function. In the leading order 
in 1/N C , we obtain (neglecting regularization) for the sum of the level and sea contributions to the 
matrix element 



(N'(jf) \&OV{0)\ N(p)) =N C J d z xe i{ ^-^ 

x J dRrj,j>T>T>(R) E & n {Z)tfdm n (Z)if>jj a r m {R) (in) 



en<f„ 



a simple expression with a clear physical meaning: in the mean-field approach the nucleon consists 
of quarks occupying all the levels with the energy e n < e va i- In deriving ( |171| ) we took into account 
the zero modes of the soliton. The rotational motion of the soliton results in the rotational wave 
functions of the soliton ip j j 3 tt 3 {R) and the isospin rotation of the matrix O — > R^OR. Similarly, 
the integral over x and exp(i(f? — p)x) appear due to the translational zero modes. In the large N c 
approach the mass of the nucleon is 0(N C ), whereas the momentum transfer in the nucleon matrix 
element is q = ff — p = 0(N®), so that we work in the non-relativistic limit \q\ -C M c \. 



Expression (171) corresponds to the leading order of 1/N C expansion. In the subleading order in 
l/iV c , there appear various corrections to this leading-order result. Below we concentrate on the 1/N C 
corrections connected to the rotational motion of the soliton. One of the reasons why we do it is that 
for some quantities the leading result vanishes identically and the non-vanishing contributions come 
from the rotational corrections. 

Note that the rotational matrix element in ( |171| ) is trivial if the quark current \l/^0\I/ is isosinglet, 
since in this case R^OR = O and the integral over R in ( |171| ) reduces to the normalization matrix 
element. In the case of the isovector current, we replace in the above formulas: 

6 -> O a = Or a , (172) 

where we separate explicitly the isospin part of O b . Hence, we have 

R ] O a R = D ab {R)Or b , (173) 

where function D a b(R) is defined by 

D ab (R) = h V {R^r a Rr b ). (174) 

We see that the rotational part of the matrix element ( |171| ) reduces to the calculation of the integral 

J VRr Nf [R(T/2)]D ab [R(0)]MR(-T/2)]e- L ^ dTn2 = J dR^R) D ab (R) ^(i?)e"T^. 

(175) 

In the left hand side, we have written the corresponding path integral over the time dependent 
orientation matrices R(t) which appears in the systematic derivation of eq. ( |171| ) from the path 
integral formulas ( |168|) , (|169|) . 



Now we turn to the rotational corrections to formula ( |1 71| ) . As it was explained earlier the angular 
velocity of the soliton is parametrically small (Q = 0(1/N C )) and we can expand in it 

1 _ _ 1 

< x ' t \d(u c ) + isv x ' t >_< x ' r| D(tg |x ' r > 

- J dVdr' < x, T \^—\x! y > z'O(r') < x 1 ,t'\^-\x,t > + . . . . (176) 



Inserting this expansion into the general formulas (|168|) , (|169|) , we arrive at the path integral con- 



taining both angular velocity fi[i2(r)] and the D-function D ab [R(0)]. As long as we are working 
with the path integral approach these two quantities are c-numbers. However, if we compute this 
path integral and express it in terms of the operator approach, we arrive at non-commuting collec- 
tive operators and the result generally depends on the ordering of these operators. This point has 
provoked a lot of discussion (Schechter and Weigel, 1995a,b), (Christov et ai, 1995a), (Wakamatsu, 
1995) and it is often misunderstood. The reason is that the canonical quantization rule is considered 
independently of the path integral. Indeed, using some "naive" ordering of the collective variables Q 
and D a b (e.g. Wakamatsu and Watabe (1993) kept the order in which these variables appear after 
the expansion ( |176|) ) one faces the problem of violation of the Pauli principle and of the G-parity 
symmetry (Schechter and Weigel, 1995a,b). Such a problem does not exist using the 1/N C expansion 
in the present semiclassical quantization scheme (Christov et ai, 1994, 1995a). In this scheme, as it 
was already shown, the canonical quantization rule is a consequence of the exactly calculated path 
integral over R with the action quadratic in Q and hence, it must be considered in the context of this 
path integral. It is well known that the general formula connecting the path integral and the operator 
approaches leads to time-ordering of the operators which correspond to the time-dependent c-number 
quantities of the path integral. Therefore, using the quantization rule (147) we arrive at the following 



result for the rotational path integral appearing in the formula for the rotational corrections to the 
nucleon matrix element of a quark current: 

J VRr N [R(T/2)] D ab [R(0)] n c [R(r)]MR(-T/2)} e~U^ 

= - l -j dRr N (R)[$(-T)D ab (R)J c + $(T) J c D ab (R)]i; N (R)e~ ii i^ lT , (177) 

sandwiched between the nucleon rotational wave functions. Using the spectral representation of the 
quark propagator flDTp we can integrate over r and finally get for the quark valence contribution 



(N'tf)\¥O a *\N(£)) vcd = N C J tfxeW-wJ dRif>* N ,{R)l [& val (x) O b D ab (R) 
^ E {[®UZ)O b $ val (x)](val\T c \n)[e(e n ) J c D ab (R) + 6(-e n ) D ab (R)J c ] 



+[$L(*) O h <b n (x)} (n\r c \val)[6(e n )D ab J c + 6{-e n ) J c D ab {R)\) U N (R) (178) 



In the case of the sea part of the matrix element ( |168|) we have an additional problem with the 
divergences of the model. The integrand in the r.h.s. of eq. ( |168| ) contains the fermionic determinant 
Det[D(Z7)] and the related generally divergent matrix element 

Sp[d<0|^|0>]=^log D 4 o . (179) 

where 

D e = D{U c )+m-iR ] OR. (180) 

As we already stated only the real part of the fermion determinant is divergent and needs to be 
regularized whereas the imaginary part is finite and it is not regularized. 

Let us first consider the sea part of the matrix element of a particular current which contributes only 
to the imaginary part of Trlog-Dg(£7). In this case the matrix element given by ( |179|) is finite and 
does not contain any regulator. Following the same steps as for the valence contribution, namely the 
expansion ( |1 76| ) keeping the terms up to linear order in Q and making use of the identities ( |130| ) and 



(|173|) , we get for the sea part 



(N'tf)\&0 a *\Ntf)) 8ea = N C [ d 3 xe^~^ / dRr N >(R)\ E & n (?) b $ n (x)D ab (R) 



i ]T -^—{j c D ab (R) [$i(x)O b <S> m (x)] (m\T c \n) + [& m (x)O b $ n (x)] J c D ab (R)}(n\r c \m) 



xtp N (R) 



;isi) 



Adding the valence contribution ( |178| ) to the sea part (|181| ), we obtain the final results for the matrix 
element of a given current whose sea part originates from the imaginary part of the effective action: 

(N'(p)\&O a V\N(p))=N c J d 3 xe^-^J dR^ Nt (R)i E & n (x) O b <f> n (x) D ab (R) 

ye n <e va i 

E —^—{i&n(x)O b $ m m (m\r c \n)J c D ab (R) + [& m (x) O b (n\r c \m) D ab (R) J c }\ 



21 e n >e val e n e m 



xip N (R) 



:i82) 



It is easy to see that expression (|182|) has the same structure as the sea contribution but with the 
valence quark included into the occupied states. 



In the case that the matrix element (|179 ) comes from the divergent real part of Trlog D^(U) it 
must be regularized. Actually, in order to make it finite it is sufficient to start with the proper-time 
regularized real part 



Re[Irlog^(C/)] reg = ~Trlog[D\(U)Ds(U)} 
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in the evaluation of the matrix element 



Sp 0(x\ 



D(U) 



\x) 
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6t{x 



-Re(Trlog D e ) 



reg 



--SpR^x^ORix^ixle-^^D^lx) - - Sp J R t (x 4 )O t i?(x 4 )(x|( J D t )- 1 e- DtD / A2 |x) (184) 
2 2. 



Here D stands for 



D = D i=0 = d T + h(U c ) + ifi 



;i85) 



Note that we assume £ (x) to be real but the matrix O is not necessarily hermitian. After quantization 
( 147|) of the time-ordered collective operators and using the spectral representation for the quark 
propagator (|96D we get the regularized sea contribution in the form 

(N'(jf)\¥O a V\N(£)) s r i a g = N C J tfxe'V-Wj dRr N '(R){Y, ^^i(^) O b <$> n {x) D ab (R) 

+77 E [ ° b ^n(x)](n\r a \m) - V [& n (x) O b $ m (f)] (m\r c \n) 



x [^ +) (e m , e n ) J c D ab (R) + K { 2 ) (e m , e n ) D ab (R)J c ] . 

The coefficient r] = 1 for hermitian matrix O and i] = — 1 for anti- hermitian one: 

6 f = r]6 . 



186) 



;i87) 



The regularization functions are given by 



(188) 



and 
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7?Wf e e )-f*HL f duJ 1 



x 



{ 



— oo 
.,2 , Ji / a2 



e -[^+e^]/A 2 



(iu + e n ){iuj' + e m ) A 2 



1 + 



iuj — e„ 
iuj' + e r . 



J 1 dae _[a(tj2+e ™ )+(1 " a)(w ' 2+e ™ )1/A2 |. (189) 



As can be seen in the present scheme, the matrix element of a given current is factorized in a quark 
and collective part. The latter determines the spin-flavor structure. The particular form of the 
matrix element depends on the hermitian properties of the current O and of the symmetry of the 
quark matrix element under m <-> n. 



4 Nucleon properties in the SU(2) NJL model 



Using the formalism developed in the previous section, we are in a position to calculate the nucleon 
observables. We start with the electromagnetic form factors of the nucleon and the related static 
properties like radii and magnetic moments. The E2/M1 ratio for the 'yN — > A transition is also 
considered. Further, the axial properties of the nucleon, namely the weak axial and pseudoscalar 
form factors and the axial and induced pseudoscalar coupling constants, are examined. The results 
for the pion nucleon form factor and the related strong pion nucleon coupling constant are also 
presented. We discuss the PCAC and the Goldberger-Treiman relation and consider the nucleon 
electric polarizability. We make some prediction for the tensor charges which is related to the quark 
transverse structure function. 

4.1 Electromagnetic Properties 



The nucleon electromagnetic Sachs form factors are related to the matrix element of the electromag- 
netic current 

f = ^tyyg^ ( 190 ) 

in the standard way: 

(N'(p')\j°(0)\N(p)) = G E tf)5j, j3 , (191) 

(N'(p')\j k (0)\N(p)) = ^ klm (Aji 3 . h <tGM{q 2 ) , (192) 
Here = 938 MeV is the physical nucleon mass, Q is the quark charge matrix 

and 

q = p — p (194) 

is the momentum transfer. For convenience, we also will use the notation Q 2 = —q 2 . Accordingly, 
the isoscalar and isovector parts of the form factors are defined by 

G E {M) = 2^i(M) + t sG%Im^ . (195) 



In order to examine the divergent structure of the nucleon matrix element of the current it is 
enough to consider the regularized determinant with a small external source £ ( |179| ). In the case of 
the electromagnetic form factors this external source is an electromagnetic field coupled to the 
current j M : 

D{U)-iO - D(U,QA fl )= l0 [ lll (d l ,-iQA ll ) + MU^]. (196) 

The proper-time regularization in the Euclidean space ( |183| ) contains operator D^D, and in order to 
preserve the gauge invariance the Euclidean electromagnetic field A^ must be consider real. Since the 
complex conjugate euclidean Dirac matrices 7* are connected to 7^ by some unitary transformation 
V (^), it is easy to show that 

Tr\ogD(U,QA^) = Tr\ogl(VT 2 )D(U,QA^(VT 2 )- 1 } = Tr\og[D(U,Q'AJ}\ (197) 



Here 



1 r 3 

«="6 + y 



(198) 



and the asterisk stands for the complex (not hermitian!) conjugation. Equality ( |197|) shows that 
the isoscalar electromagnetic form factors originate from the imaginary part of Tr log D(U, QA^), 
whereas the isovector form factors are generated by the real part of Tr log D(U, QA^). Since the 
imaginary part of Trlog.D(?7, QA^) is ultraviolet finite, for the isoscalar form factors we use directly 
the non-regularized expression ( |182| ) , whereas the sea contribution to the isovector form factors must 
be regularized ( |186|) . 



Under the transformation Vt 2 ( |197|) we also have 



(Vr^Vr 2 )' 1 = tJ, (yrVtW 2 )" 1 = -(^T = "(r a f, (199) 

{VT 2 )h{VT 2 y l = h* = h T , (200) 

where the superscript T stands for the transposition in both the matrix indices and the coordinate 
space. The last equation shows that the eigenf unctions of h can be chosen in such a way that 

(yr 2 )$ m (f) = $;(f). (201) 

Hence, the following relations hold: 

[^n{x)lol^ m (x)}(m\T a \n) = 7] [$J ra (f)7o7 At $ n (f)](n|r a |m) , (202) 

and 

[<^)7o7/^ a $ m (£)]Hr>> = -rj l$l(x) lQlfl T a <5> n (x)](n\T a \m) (203) 
where the values for 7] are determined by the hermitian properties of i r )o r )^T a : 

(204) 

The rotational matrix elements, which we need, can be easily computed with the relations 

(TlJ' 3 \J a \T 3 J 3 ) = l -8 T , T3 (T a )j,j 3 , (205) 

{T'^\D ab {R)\T 3 J 3 ) = -i (r a ) nn (T b )j, j3 . (206) 

Now using the symmetry properties of the quark matrix element we can apply the general formulas 
for the nucleon matrix elements derived in the previous section to the case of the electromagnetic 
current ( |190| ) for evaluation of the nucleon electromagnetic form factors (|191|) . 



For the isoscalar electric form factor we obtain a simple result: 



g|=V) = -y I d3 ^1 E - E ^W^SH*)} , (207) 



where the vacuum contribution is subtracted. Since the rotational corrections exactly vanish this 
form factor contains only leading order terms . The baryon number is given by the form factor at 
q 1 0: 

B = G T E =°(0) = /d 3 x{ £ $t (£)<,>„(£) - £ $(°)t(f)$W(f)} = |d 3 x$L(^)^(^) = l-(208) 

e n <e val e (0) <Q 

As it was mentioned in the previous chapter it originates from the imaginary part of the effective 
action and is determined by the valence quarks. 



For the electric isovector form factor we obtain: 



E 



1 



}i(x)r a ^ val (x)](val\r a \n) 



(209) 



where the proper-time regulator IZf is the same as in the case of the moment of inertia (|138|) . Hence, 
the isovector electric form factor has a proper normalization (electric charge): 



G T E =\0) = 1. 



(210) 



It should be mentioned that in contrast to the isoscalar electric form factor the isovector one vanishes 
in leading order (Q° ~ AT?) and the only non-zero contribution comes from the rotational corrections 
(next to leading order terms in Q). 

As in the case of the isovector electric form factor, the isoscalar magnetic one also vanishes in the 
leading order in Q, and we find non- vanishing contributions only from the 1/N C rotational corrections: 



GlTV) = ^^^/d 3 xe- E 



[$J 7l (f)7°7 fc $ n (f)] {n\r a \m} 



(211) 



e m<e va i 



-m L n 



The calculations of the magnetic isovector form factor are more involved and the final expression 
G T M =\q 2 ) = ^Le^ *L J d 3 xe-{[$L(^)7°7 fc r^ m ,(f)] - E^UO H (x) 7 ° 7 V^ n (f) 



E^M 2 (^,en) [^(x) 7 °7V<fn(f)] (n|r» 



- E 



sign(e n J 



l(x)^ k r a ^ai(x)}(val\r c \n) 



(212) 



includes leading order terms (the first two terms in the rhs) as well as rotational corrections. Accord- 
ingly, two different regularization functions appear: 



e n 7 du 



K M = -K l(£ „) = ^ J ^e- 

1/A 2 



(213) 



v a ( x_ 1 f 1 d(3 {l-(3)e m -(3e n 
M2(em ' ^ ~ ^ Jo J^W) (1 " /% + ^ 



-[(l-/3)^+/34]/A 2 



(214) 



Eqs. ( |207| ), ( |209| ), ( |2 1 1| ) and ( |212| ) are our final expressions for the electromagnetic isoscalar and 
isovector form factors in the NJL model. According to (|195|) the proton and neutron form factors are 
expressed in terms of the isoscalar and isovector form factors as follows 
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Figure 9: Electric (left) and magnetic (right) form factors of the nucleon as functions of Q 2 for 
three different values of M: 370 (dashed line), 420 (solid line) and 450 (dash-dotted line) MeV in 
the SU(2)-NJL model (Christov et al, 1995). The experimental data in the case of the electric form 
factor are from (Hohler et al, 1976) for the proton and from (Platchkov et al, 1990, Eden et al, 1994, 
Meyerhoff et al, 1994) for the neutron. For the magnetic form factors they are from (Hohler et al, 
1976, Bruins et al, 1995). The theoretical magnetic form factors are normalized to the experimental 
values of the magnetic moments. 



The nucleon electric and magnetic form factors are displayed in figs.[| a) and b). The theoretical 
curves (Christov et al, 1995b) are given for three different values of the constituent quark mass, 
370, 420 and 450 MeV. The magnetic form factors are normalized to the experimental values of the 
corresponding magnetic moments at Q 2 = 0. 



With one exception, the neutron electric form factor, all other form factors agree with the experimental 
data quite well for the constituent quark mass around 420 MeV. It should be noted, however, that 
the neutron electric form factor, which deviates from the experimental data, is most sensitive to the 
used approximations and numerical errors. According to the formula (|216 ) the form factor has been 
calculated as a difference of the electric isoscalar and electric isovector form factors. Both form factors 
are of order of one, whereas the resulting neutron form factor has experimental values of order 0.04, 
i.e. about 4% of the value of its components. It means that the numerical accuracy together with 
the used large 1/N C approximation behind the model are strongly magnified and could in principle 
be an explanation for the deviation from the experimental data. It should be also stressed that the 
experimental data (see Platchkov et al, 1990), available for the neutron electric form factor, are 
strongly model dependent and a different N - N potential used in the analysis of the data can lead to 
an enhancement of the experimental numbers by more than 50 



/o. 



The mean squared radii and the magnetic moments are obtained from the form factors: 
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(217) 



(218) 



and the results (Christov et al, 1995b) for three values of the quark constituent mass, 370, 420 and 
450 MeV, and pion mass = 140 MeV are presented in table The value ~ 420 MeV is favored 
which agrees with the conclusion drawn from the form factors. With the exception of the neutron 
electric squared radii, to which remarks similar to the case of the neutron electric form factor are 
valid, the contribution of the valence quarks is dominant. However, the contribution of the Dirac 
sea is non-negligible and it varies within the range 15 - 40%. As can be seen, the numerical results 
for the nucleon N-A mass splitting, the mean squared proton, isovector electric radii as well as the 
g-dependence of the proton electric and magnetic and the neutron magnetic form factor differ from 
the experimental data by no more than about ±5%. Finally, for the magnetic moments the results 
are by 25-30% below their experimental values. However, the experimental ratio \x v j \x n is almost 
exactly reproduced which is not the case of other chiral models of nucleon. 



The magnetic isovector form factor is the only one which includes non-zero contributions in both 
leading (iV°) and next to leading order (1/N C ) rotational corrections. However, the numerical cal- 
culations show that the (1/N C ) rotational corrections do not affect the g 2 -dependence (slope) of the 
form factor but rather the value at the origin, G^p 1 (0) = /i T=1 which is the isovector magnetic 
moment. It is not surprising, since (as can be seen from eq. ( |212j )) in both the leading and next to 
leading order terms the shape of the wave functions $ n determines the g 2 -dependence of the form 
factor, whereas the value at q 2 = depends on the particular matrix elements included. In leading 
order the isovector magnetic moment is strongly underestimated (Wakamatsu and Yoshiki, 1991). 
The enhancement for this quantity due to the 1/N C corrections is of order (N c + 2)/N c (Christov et 
al, 1994) and it improves considerably the agreement with experiment. It is also important to note 
that since the 1/N C rotational corrections have the same spin-flavor structure like the leading term, 
they do not violate the consistency condition of Dashen and Manohar (1993) derived in the large- N c 
of QCD. The latter means that other 1/N C corrections (e.g. meson loops) should contribute to both 
the leading and next to leading orders. 



Table 1: Nucleoli observables calculated in the SU(2) NJL model with the physical pion mass. 



Quantity 


Co 

370 MeV 


ustituent Quark M 

420 MeV 


ass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< r 2 >t=o [fm S 


-i 


0.63 


0.05 


0.52 


0.07 


0.48 


0.09 


0.62 


< r 2 > T =i [fm 2 ] 


1.07 


0.33 


0.89 


0.41 


0.84 


0.45 


0.86 


< r 2 > p [fm 2 ] 


0.85 


0.19 


0.70 


0.24 


0.66 


0.27 


0.74 


< r 2 > n [fm 2 




-0.22 


-0.14 


-0.18 


-0.17 


-0.18 


-0.18 


-0.11 


fi T=0 [n.m. 




0.68 


0.09 


0.62 


0.03 


0.59 


0.05 


0.88 


/xt=i [n.m.] 


3.56 


0.77 


3.44 


0.97 


3.16 


0.80 


4.71 


fi p [n.m.] 


2.12 


0.43 


2.03 


0.50 


1.86 


0.43 


2.79 


//„ [n.m.] 


-1.44 


-0.34 


-1.41 


-0.47 


-1.29 


-0.38 


-1.91 




1.47 




1.44 




1.44 




1.46 


< r 2 [fm 2 ] 


1.08 


0.32 


0.66 


0.28 


0.56 


0.25 


0.74 


< r 2 [fm 2 ] 


-1.17 


-0.51 


-0.65 


-0.31 


-0.52 


-0.24 


-0.77 


M A - Mjv [MeV] 


213 




280 




314 




294 


fi'A 


1.26 


0.08 


1.21 


0.11 


1.13 


0.06 


1.26 
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Figure 10: Charge density distribution of the proton (lower) and neutron (upper) for the constituent 
quark mass M = 420 MeV taken from (Christov et al, 1995b). 




M [MeV] 

Figure 11: Electric charge radii of proton and neutron as functions of the constituent quark mass M 
in the SU(2) NJL model. The valence and sea parts are marked by the dashed and dash-dotted lines, 
respectively (Christov et ai, 1995b). 

In fig.|l(] we plot the typical proton and neutron charge distributions from (Christov et ai, 1995b) 
which is similar to those of Wakamatsu (1992) and Gorski et al. (1992) for the constituent quark 
mass M = 420 MeV. For the proton we have a positive definite charge distribution completely 
dominated by the valence contribution whereas in the case of the neutron the Dirac sea is dominant. 
In accordance with the well accepted phenomenological picture, one realizes a positive core coming 
from the valence quarks and a long negative tail due to the polarization of the Dirac sea. Using the 
gradient expansion the latter can be expressed in terms of the dynamical pion field - pion cloud. It 
is reflected on the proton and neutron charge radii shown on fig.|TT[ In the case of proton the sea 
contribution is about 30% whereas for the neutron charge radius the negative sea part is dominating 
and the valence contribution is negligible. 

The theoretical values of the static nucleon properties but with the physical pion mass set to zero can 
be found in (Christov et ai, 1995b) and the corresponding electric form factors in (Gorski et ai, 1992). 
The chiral limit (m^ — > 0) mostly influences the isovector charge radius. In fact, as it is expected 
(Beg and Zapeda 1972, Adkins et ai, 1983), the isovector charge radius diverges in chiral limit and 
the calculations with zero pion mass confirms it. Another quantity strongly influenced by the chiral 
limit is the neutron electric form factor. For the m n — > the discrepancy from the experiment is by 
almost a factor two larger than in the case ^ 0. The other observables differ in the chiral limit 
by about 30%. The reason is that in chiral limit one observes much larger contributions from the sea 
effects, up to 50% of the total value. One concludes that using the physical pion mass one is able to 
get a much better agreement with the experimental data. 



4.2 E2/M1 ratio for the jN — > A transition 



Here, we consider the ratio of electric quadrupole to magnetic dipole amplitude E2/M1 ratio for 
the process 7 + iV — > A. It is a quantity of a current interest, since it is sensitive to a presence of 
charge deformations in the baryon structure and can be measured experimentally. In a very recent 
7r °(+)_photoproduction experiment (Beck et al, 1995) performed at MAMI, Mainz the direct model- 
independent estimate of the ratio E2/M1 = (—2.4 ± 0.2)% indicates a presence of an oblate charge 
deformation in the nucleon or/and delta, and as such it imposes a strong constraint for the effective 
models of baryon structure. Explicitly one has 



E2 _ 1 M E2 (k, A = +1 ; p(J 3 = -§) - A+(J 3 = +§)) 



(219) 



Ml 3 M mfo A = +1 ; p(J 3 = -\) -> A+(J 3 = + 

of the electric quadrupole amplitude M E2 to the magnetic dipole one M M1 . Here k is the momentum 
of a photon of helicity A in the A rest frame: 

Both amplitudes are related to the corresponding matrix element of the isovector current: 

jS(z) = *(z)-f^*{z). (221) 

For the electric quadrupole amplitude M E2 , according to the Siegert's theorem (see for instance 
(Eisenberg and Greiner, 1970)), one can use the zero component jj> =1 as well the space component 
V ■ j't=i- Similarly to (Wirzba and Weise, 1987), we decide to express the amplitude M E2 via the 
N — A transition matrix element of jj<=i (charge density) of the current jj<=i- 

M E2 (k, X;N^A) = Vl^ J d 3 x (A\f 3 (x)\N)Y 2X (x)j 2 (kx), (222) 

The reason is that in the large N c treatment this quantity can be calculated (Christov et al, 1995b) 
directly in terms of quark matrix elements, whereas the matrix element of V • jr=i is suppressed 
additionally by 1/N C and in order to evaluate it one should use supplementary the saddle-point 



equation (|122|) for which the ansatz ([127]) is apparently not a solution 



The amplitude M is directly related to the N— A transition matrix element of the space components 

M M1 (k, A; N - A) = -a| J d 3 x(A\(x x J 3 )a|AWi(^)- (223) 

The matrix elements the current j M are calculated in a way similar to those of the nucleon electro- 
magnetic form factors. 

In the particular case of A = +1, Jf = — | J% — +|, we have 

M E2 = ^ I drr 2 Ukr)p E l(r) , (224) 



and 

M M1 = -3 / drr 2 3l (kr)p%l(r) , (225) 



respectively. The corresponding transition charge density p^A^) nas a m o re complicated structure 

(Watabe et al, 1995a) including a spherical harmonics tensor which acts on the quark wave 
function projector for the charge deformation. 

In the fcr < 1 approximation, one can relate the ratio E2/M1 to the electric quadrupole iVA 
transition moment: 

E2 1 < Q zz > NA fnna\ 

Ml = 2 MNk /ijvA ' (226) 



Table 2: Ratio E2/M1 and some related observables, calculated in the SU(2) NJL model (Watabe et 
al, 1995a) for three different values of the constituent mass M = 400,420 and 450 MeV, compared 
with experimental values. 









Constituent quark mass M 






Quantity 


400 MeV 


420 MeV 


450 MeV 


Exp 






total 


sea 


total 


sea 


total 


sea 




< r 2 >i 


=i M 


0.88 


0.35 


0.84 


0.37 


0.79 


0.41 


0.86 




n.m.] 


2.34 


0.57 


2.28 


0.58 


2.20 


0.60 


3.33 


M A - M N [MeV] 


255 




278 




311 




294 


< Qzz >an [fm 2 ] 


-0.020 


-0.014 


-0.020 


-0.015 


-0.021 


-0.016 


-0.026 


E2/M1 [%} 


-2.19 




-2.28 




-2.42 




-2.4±0.2 



where Q zz >na is the electric quadrupole transition moment and una is the transition magnetic 
moment. 

In table |2|, the results (Watabe et al, 1995a) for the ratio E2/M1 as well as for some related observ- 
ables, namely the isovector charge m.s. radius, the N — A transition magnetic moment /j-na, the iV — A 
mass difference, and the quadrupole electric transition moment (Q zz )na, are given for three different 
values of the constituent quark mass M in comparison with the experiment. For the constituent 
quark mass M around 420 MeV the E2/M1 ratio is between —2.5% and —2.3% quite in agreement 
with the experiment estimate (Beck et al, 1995). From the relation fl226|) we also present an estimate 
for the electric quadrupole transition moment < Q zz >na= —0.026 using the experimental values for 
E2/MI = —2.4 ±0.2 and Una — 3.3 which is not far from the model prediction —0.02. This negative 
non-zero value indicates a charge deformations of oblate type in the nucleon or/and delta structure. 
As can be seen from table |2| the dominant contribution to the < Q zz >na i n the NJL model comes 
from the Dirac sea. It means that the main charge deformation is due to the polarized Dirac sea, 
whereas the valence quarks are almost spherically distributed. Since using the gradient expansion, 
the polarization of the Dirac sea can be expressed in terms of the dynamical pion field - pion cloud, 
one can think of the nucleon or/and of the delta as consisting of an almost spherical valence quark 
core surrounded by a deformed pion cloud. One should note, however, that in the present approach 
the Delta is a stable particle and has no width which is in fact a result of the large N c treatment. 
The implications of this approximation are discussed by Cohen (1995). 

We remind that for the same values of the constituent quark mass M ~ 420 MeV the nucleon prop- 
erties (including also the nucleon form factors) are reproduced fairly well. The only exception is the 
N — A transition magnetic moment which similar to the isovector magnetic moment is underesti- 
mated by 25%. The results for other observables in table |2] show an overall good agreement with the 
experiment. 



4.3 Axial properties 

The axial and pseudoscalar form factors reflect the nucleon structure seen by the weak probes, and as 
such provide an important test for any effective model of the nucleon. Most of the theoretical efforts 
are concentrated on the axial whereas the pseudoscalar one is believed to be generally understood in 
the context of the chiral symmetry. 



The general decomposition of the nucleon matrix element of the axial current 

A^x) = ^(x) l0l , l5 ^(x) (227) 
in terms of the axial Ga and pseudoscalar form factor G p reads 

(N(P,OK(0)\N(p,0) = u(P,0[G A (q 2 h» + ^f\]7 5 f ^,0 • (228) 



where £ stands for both spin and isospin. 



In the large N c limit the time component of the axial current is suppressed by 1/N C compared to the 
space one. It means that in the present large N c treatment the matrix element of the space component 
of the axial current, which we use for the evaluation of the weak form factors, is "enhanced". 

In the calculations, one makes use of the general expressions (|178|) and ( |186| ) for the nucleon matrix 
element derived in the previous section. As for the case of the isovector magnetic current the Dirac sea 
contribution comes from the real part of the effective action and needs regularization. Also i^o^y^T 0, is 
hermitian (77 = 1) and the quark matrix elements have the same symmetry properties under n <-> m 
as those in the isovector magnetic case. Thus, the matrix element of the space components of the 
axial current A 3 includes leading order terms ~ Q° as well as next to leading order ones ~ Q (1/N C ) 



(N(p') \A k 3 (x)\N(p)) = -N C J d 3 xe^j(«I>L(^7V75r3^(x)) 

(*L(^)7°7 fc 75^$n(s))(^ \r c \ val) 



1 cb3 



-E<i( e n)K(^)7V75r 6 $ n (f)) + 



E si S n ( e r 



.n^val €n 6val 



(229) 



where the first and third terms are valence quark contributions in leading and next to leading order 
in angular velocity, respectively. The other terms represent the divergent Dirac sea part. 

Using result Q229|) after some straightforward calculations one can express the axial and pseudoscalar 
form factors in terms of two axial densities A Q (r) and A 2 (r) (Watabe et ai, 1995b) 

G A {q 2 ) = ^r j r 2 dr [jo(qr)A (r) - j 2 (qr)A 2 (r)} . (230) 

and 

\ l '2 .. 1/ 2 1 / ,. 

Gp(q 2 ) = - e(e+ N Mn) J r 2 dr Jo (qr)A (r) - 24-f (1 + J r 2 dr 32 {qr)A 2 {r) , (231) 

where 



A (r) = (N\Al\N) , (232) 
A 2 (r) = (N\±Al-Al i x*x 3 \N}. (233) 



In the above expressions and E are the nucleon mass and energy E = yM^ + cp/4 in the Breit 
frame, respectively. 



Axial form factor and the axial- vector coupling constant. The results for the axial form factor (Watabe 
et ai, 1995b) are shown in fig.|12| for four different values of the constituent mass M, 360, 400, 420 
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Figure 12: Axial form factor in the SU(2) NJL model (Watabe et ai, 1995b) in leading (lower curves) 
and next to leading order (upper curves) in angular velocity in comparison with the experimental 
data (Baker et ai, 1981, Kitagaki et ai, 1983). 



and 440 MeV. The main contribution to the axial form factor in the leading as well as in the next 
to leading order in angular velocity comes from the valence quarks and the Dirac sea contribution is 
almost negligible. As for the other nucleon properties a value for constituent mass around 420 MeV 
is preferred. For this mass value the theoretical curve agrees well with the experimental data without 
any scaling. On fig.[12| we also show the leading order results (Meissner Th. and Goeke, 1992). As 
can be seen the two groups of curves differ in both magnitude and slope. The dipole masses (see 
table |3|) for the axial form factors, calculated in leading order, are slightly larger. 



The axial charge radius calculated as 

= G^)/ r4dr[ ^ (r) + ^ !(r)l + 4^' (2M) 

is presented in Table [3] together with the values from the dipole fit (in brackets) in comparison with 
the estimate from the experimental dipole fit. As can be seen the numbers directly calculated from 
the corresponding densities Aq and A2 are very close to those of the dipole fit which is an indication 
that even at small q 2 the dipole fit is a good approximation to the theoretical curves. 



6 dG A (Q 2 ) 
G A (0) dQ 2 



The axial-vector coupling constant is given by 

g A = G A (0) = 2jd 3 x< N\A 3 3 (x)\N > . (235) 

which means that only the density A$ contributes (|230|) . The strong underestimation of g A in most 
of the chiral models of nucleon is a problem which have attracted a lot of attention and theoretical 
efforts in both the Skyrme-like models (see Holzwarth, 1994 and references therein) and the NJL 
model (Wakamatsu and Watabe, 1993, Blotz et ai, 1993c, Christov et ai, 1994, 1995a,b, Wakamatsu 
1995, Schechter and Weigel, 1995). Wakamatsu and Watabe (1993) made an important observation 
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Figure 13: Axial coupling constant gA evaluated in the leading (0°) and next to leading order + 
in the SU(2) NJL model as a function of the constituent quark mass M (Christov et ai, 1994). 



that since after the canonical quantization, the collective operators do not commute, it may lead to 
non-zero rotational corrections to gA which could resolve the problem of the underestimation of Qa 
in the leading order in Q (Wakamatsu and Yoshiki, 1991, Meissner Th. and Goeke, 1991). However, 
in their calculations, the non-zero result is due to a chosen order of the collective operators being not 
justified by path integrals or many-body techniques. As a consequence their expression violate the 
G-parity symmetry (Schechter and Weigel, 1995) and the Pauli principle (even though, numerically, 
this Pauli violating contribution turns out to be only a tiny fraction of the valence contribution), 
and the rotational (1/N C ) corrections from the Dirac sea are missing. As it was already mentioned in 
chapter 3 such a problem does not exist in the present semiclassical quantization scheme (Christov 
et ai, 1994) and it can be easily checked directly from (|229|) (Christov et ai, 1995a). 



On fig.0 we present the results for the axial-vector coupling constant in leading and next to leading 
order in angular velocity Q as a function of the constituent mass M (Christov et ai, 1994). As 
can be seen, in both the leading and next to leading order the valence part dominates and the 
Dirac sea contribution is almost negligible. The enhancement due to the 1/N C rotational corrections 
improves considerably the agreement with experiment and for the constituent quark mass of about 
420 MeV the experimental value of gA is almost exactly reproduced. It is interesting to notice that 
this enhancement numerically is very close to 



N c +2 
N c 



In the present large N c treatment, only the rotational 1/N C corrections coming from the rotational 
zero modes are taken into account . Apparently, they are only a part of the existing 1/N C corrections 
and the inclusion of all 1/N C corrections could in principle changes the obtained results. In the 

) 



8 2 S 



111 



present approach the only source of additional 1/N C corrections is the term |Trlog 
which represents the meson quantum (loop) effects. One generally expects that the meson loops 
would affect mostly the contribution of the polarized Dirac sea (represents the meson clouds) and not 
much that of the valence quarks. It means that for the physical quantity which are dominated by the 
valence quarks like the axial vector coupling constant one should not expect large changes due to the 



Table 3: Axial properties of the nucleoli, calculated in the NJL model (Watabe et al, 1995b) for four 
different values of the constituent mass M = 360, 400, 420 and 440 MeV, compared with experimental 
values. The obtained dipole mass for the axial form factor and the axial m.s. radius are given in both 
leading and next to leading order in angular velocity. The radius is also calculated from the dipole 
fit (in brackets). The induced pseudoscalar coupling constant gp is also presented. 



Quantity 


Const 

360 


ituent quar 

400 


k mass M [ 

420 


MeV] 

440 


Exp 


Mf [GeV] 


0.91 


1.01 


1.04 


1.072 


1 05+ - 12 

i - uo -0.16 


Mf +QL [GeV] 


0.85 


0.96 


0.99 


1.03 


1 05+ - 12 

i - uo -0.16 


< r 2 > A [fm 2 ] 


0.70 (0.65) 


0.51 (0.51) 


0.46 (0.47) 


0.43 (0.44) 


42+°- i8 


9a 


0.83 


0.80 


0.78 


0.77 


1.26 


9a 


1.31 


1.24 


1.21 


1.18 


1.26 


9p 


6.09 


6.05 


6.01 


6.00 


8.6 ± 1.9 



1/N C corrections from the meson loops. Additionally, the rotational 1/N C corrections obey G-parity 
symmetry and fulfill the consistency condition of Dashen and Manohar (1993) derived in the large- N c 
of QCD. The latter means that the other 1/N C contributions should contribute to both the leading 
and next to leading order terms, and the enhancement due to the rotational corrections will survive 
in a strict next to leading order in 1/N C . 

Praszalowicz et al, (1995) studied the behavior of Qa in the limit of both small and large soliton 
size R. We remind that at MR 3> 1 the NJL soliton shows a structure similar to the topological 
soliton of the Skyrme model. At small R limit the contribution of the Dirac sea vanishes and we 
have a physical picture of three quarks of constituent mass M similar to the one of the constituent 
quark model (CQM). The results of (Praszalowicz et al., 1995) are shown in fig. YA. One immediately 



sees that only in the case with the rotational 1/N C corrections included, the results for g& of 
CQM can be reproduced. It should be stressed that since in the NJL model the nucleon appears 
as a non-trivial bound state of valence quarks and polarized Dirac sea, and the CQM is based on 
pure symmetry assumptions concerning the wave functions, this agreement is by no means obvious. 
Apparently, the semiclassical quantization with expansion up to first order in the rotational velocity 
for the observables and up to the second order for the energy provides a consistent and powerful 
method to describe the rotational degrees of freedom. In the large R limit, the rotational 1/N C 
corrections exactly vanish and we recover the Skyrme case of no rotational 1/N C corrections. 



Pseudoscalar form factor. For the pseudoscalar form factor the chiral symmetry predicts (Amaldi et 
al, 1979) to be entirely dominated by the pion pole and as such to be very sensitive to the pion content 
of the nucleon. The model result for M = 420 MeV is shown in fig.|l5|. The experimental points are 
from the first exclusive experiment (Choi et al, 1993) determining independently both weak form 
factors. The muon capture point (Esaulov et al, 1978) is included as well. We also present the pion 
pole dominance prediction (dashed line). As can be seen, despite the fact that the magnitude of Gp is 
underestimated, its g 2 -behavior is in a good agreement with both the experimental points and the pion 
pole dominance. As it is expected, the pseudoscalar form factor is entirely dominated by the polarized 
Dirac sea (the main contribution comes from A2) and the valence contribution is almost negligible. 
We remind that using the gradient expansion the polarization of the Dirac sea can be expressed in 
terms of the dynamical pion field - the pion cloud. It means that even in the approximation of no 
meson loops the model is able to account for an essential part of the pion dynamics relevant for the 
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Figure 14: Axial vector coupling constant Qa as a function of the soliton size R in the leading 
and next to leading order Q° + Q 1 for a fixed profile function in the SU(2) NJL model (Praszalowicz 
et al, 1995). 



nucleon structure. It should be also noted that the pseudoscalar form factor is very insensitive to the 
particular value of M and that the enhancement due to the rotational corrections in the case of Gp is 
much smaller than for the axial one. The latter feature is related to the fact that the Gp is dominated 
by the Dirac sea. Since the rotational corrections are ~ 1/1 (see eq. ([229| )) and the moment of inertia 
/ is dominated by the valence quarks, one expects that these corrections in the case of Dirac sea 
should be much smaller than in the case of the valence contribution. 



The calculated values of the induced pseudoscalar coupling constant defined at the muon capture 
point 



9p 



2M 



N 



-G P (q 2 = -0.88MJ 



(236) 



are presented in table |3|. As experimental value we take the value at the muon capture point directly 
from the experiment (Esaulov et al., 1978). As can be seen, the model calculations yield numbers 
for gp which are slightly outside the experimental error bars. If we use the estimate of Bernabeu 
gp = 8.2 ± 2.1 (Bernabeu, 1978), extracted from the same experimental data (Esaulov et al., 1978), 
the theoretical predictions are even better. Apparently, despite the agreement in the g-dependence 
on fig.0 this underestimation shows that some physics is still missing in the present model picture. 
Indeed, in the approximation of large N c , used in the calculations, the meson quantum (loop) effects 
are not taken into account. Since gp is dominated by the Dirac sea, this could be the reason (Bernard 
et al, 1994) for the underestimation of gp. 



4.4 Pion nucleon form factor 



We parametrize the pion nucleon vertex in terms of the pion nucleon form factor G^nn and relate it 
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Figure 15: Pseudoscalar form factor of the SU(2) NJL model in leading and next to leading order in 
angular velocity in comparison with the experimental data (Choi et al, 1993, Esaulov et al, 1978). 
The pion pole dominance fit (dash-dotted line) is also presented. 



to the matrix element of the pseudoscalar density: 



{N(p',?)\*i l5 T a *\N(p,£)) 



mil. 



GwNn(c[ 2 ) 



m 2 ^ 1 "* ' m 2 — q 2 



(237) 



where q 2 = p' 2 — p 2 and the denominator m\ — q 2 represents the pion pole. The factor m ^ ensures 
that the definition of G^nn is for a physical pion field: 



(238) 



which is consistent with 



Using the same techniques as before, we evaluate the pion nucleon form factor G^nn and the final 
results reads 



GW(g 2 ) = 2 J^L<h^ {m l _ I d 3 xjl{qr) < iV| 7 ° 75 x 3 r 3 |iY > 



q m 3 

The pseudoscalar density < N\^/ / -f 0/ -f5X 3 T 3 ^/\N > has the following structure 

1 



(239) 
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6 r n 
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[$j Tl (f)7°7 5 i- a — $„(£)] (n\r c \m) 



M2\^rai e n 



- E 



[$t (f) 7 ° 75 r a -$ m ,(f)] (val\r c \n) 



(240) 



Table 4: Monopole cutoff mass and g^NN coupling constant in the SU(2) NJL model 



Quantity 


a 

370 


f [Me< 

420 


A 

450 


Phen 


[MeV] 


835 


1086 


1242 


« 0.9 GeV 


<ffl(0) 


8.1 


7.8 


7.5 


13.1 


*S5T } (o) 


10.2 


10.3 


10.3 


13.1 



including both leading and next to leading terms. 

The pion nucleon coupling constant gvjvjv(O) at g 2 = is given by 

9*nn(0) = GW(0) = 2M N ^ / d 3 x < N\y l5 x 3 T 3 ^\N > , (241) 

mo J 

and it differs slightly (at the level of about 5%) from the physical pion nucleon coupling constant 
g^NN^V) defined via the residue at the pole in (|237|) . 

The calculated pion form factor is presented on fig.|l6| for M = 420 MeV together with separated 
valence and Dirac sea contributions. We plot also the corresponding monopole fit: 

2 A 2 m 2 

GttNn(<1 ) = gnNN ,™° W f" , (242) 

mon Q 

with g n NN = g-nNN^V) being the on-shell pion coupling constant. The monopole fit with a monopole 
cutoff mass of order of 1 GeV is commonly used for phenomenological parametrizations of the pion 
nucleon vertex. The extracted values for the monopole cutoff mass are given in table f|. As can be seen, 
the calculations deviate from the monopole fit. Although the form factor is a smooth monotonically 
decreasing function of the momentum transfer q 2 , the valence and the sea contributions have rather 
different and to some extent opposite behavior. The sea contribution increases rapidly with q 2 — > 
and hence, entirely dominates at small q 2 . Keeping in mind that in the model the polarized Dirac 
sea plays the role of the pion cloud, it can be easily understood as a dominance of the pion-pion 
resonance interaction at small q 2 . On the other hand, at larger q 2 this interaction is suppressed and 
the valence part contribution, almost negligible at small q 2 , increases with q 2 and becomes dominant. 

The obtained values of gvjviv(O) evaluated in leading as well as in the next to leading order in Q are 
presented in table |. As in the case of the axial coupling constant g^, the strong pion nucleon one 
g-xNN is also strongly underestimated in leading order. However, despite the enhancement due to the 
rotational 1/N C corrections, in the next to leading order it is still underestimated by 20 %. Similarly 
to gp the strong constant g n NN is dominated by the Dirac sea, and the meson loop effects (1/N C 
corrections) could be of importance. 



4.5 PCAC and the Goldberger-Treiman relation 

Because of the non-zero current mass m the axial current is only partially conserved ( |T8"D (PCAC). 
According to the PCAC hypothesis its divergence is related to the pion field 
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Figure 16: Pion nucleon form factor calculated in the next to leading order in angular velocity for 
M = 420 MeV in the SU(2) NJL model. The valence (dashed line) and Dirac sea (dash-dotted line) 
contributions as well as the monopole fit are presented. 



which we have already used in the case of the pion nucleon vertex ( [237] ). However, as we discuss 
below, this identity is not necessary valid in the case of nucleon matrix element 

<f < N(p')\A«\N(p) >= m < N(p')\m l5 T a V\N(p) >= m\U < N(p')\>K a \N(p) > , (244) 

evaluated in a particular model. 



In the present model, the nucleon appears as a large N c soliton (localized mean-field solution) and 
the zero modes are taken into account in order to assign proper quantum numbers to the soliton 
solution. In the many-body theory, it is known as projection after variation (Ring and Schuck, 1980). 
In evaluating fl244j) two problems arises: First, the term q°(Ao) is by two orders in 1/N C suppressed 
compared to q l (Ai). Hence, one cannot expect the NJL model to describe q°(A ) properly. Second, 
one also cannot evaluate q L (A i ) systematically. This is due to the fact that one needs the Dirac 
equation ( |9"7| ) in order to replace i^d^n by e n — •y MU^ 5 where U c is the selfconsistent solution 



of large N c saddle-point equation (|122|) . This, however, is inconsistent, since e a — MU C is evaluated 
without rotational 1/N C corrections, whereas the rhs of fl244j ) are evaluated including those corrections. 
Actually the mismatch between the lhs and rhs of fl244j ) tells something about the theoretical errors 
we do by taking from the 1/N C corrections only those which correspond to the rotational zero modes. 
This mismatch may also be seen in connection with a virial theorem used in the chiral sigma model 
by Birse and Banerjee (1985), Birse (1986, 1990) and Fiolhais et al. (1988). 



Let us estimate to what extend the violation of the equality ( p44j) will affect the axial coupling 
constant qa- In the case of a finite pion mass, qa ( |235| ) can be rewritten as 

g A = 2jd 3 x< N\Al(x)\N >= -2 J d 3 xx 3 d» < N\Al(x)\N > . (245) 



Using the PCAC ( |238| ) we can obtain another estimate g\ from the pion field (or from the pseudoscalar 



density) : 

g \ = mlf w I d 3 xx 3 < iV|7r 3 (f)|iV >= m / d 3 x < A^x 3 7 5 T 3 (f)^|iV > (246) 



In the case that the equality (|244|) is valid, we would have g\ = g&. The numerical calculations, 
however, show that g\ is by 15 % lower than g& calculated directly from the matrix element of 
the axial current. The reason for this discrepancy can be easily seen if one calculates explicitly 
< iy|<9 M A 3 (f)|iV > making use of (|7]) with (pjgj ): 

g A -g\ = N C M J d 3 x[< N\^ix 3 l5 r 3 ^\N > S c (r)- < N\^\N > P c (r)} . (247) 

where the scalar S c (r) and the pseudoscalar P c {r) mean-field densities are selfconsistent solutions of 
large N c saddle-point equation (|122|) : 

S c (r) = McosO c (r) and P c (r) = M sin c (r) . (248) 



Using that the pseudoscalar density ( |24(J| ) contains non-zero both leading and next to leading terms 

< A^7 5 x 3 r 3 ^|iV >= P c (r)+ < iV|^75S 3 r 3 ^|A^ > ( ^ } , (249) 

whereas the scalar density is given only by the leading order term S c (r), and also that both S c (r) 
and P c (r) are solutions of (|122|) , we get 



9A -gl = N C M J d 3 x < N\m l5 x 3 r 3 ty\N > {Ql) S c (r) > . (250) 

Obviously, in the leading order (0°) the two estimates, g A and g\, coincide and the PCAC equality 
( 244 ) is valid. However, in the next to leading order they differ and it is probably due to the fact 
that no rotational 1/N C corrections are present in the saddle-point equations. In the present scheme 
the PCAC relation ( |244| ) is fulfilled within 15 % "measured" in terms of the axial coupling constant. 

Alkofer and Weigel (1993) proposed as a solution to this problem to use the rhs of ( [24 7| ) as a new 
"equation of motion" which respects the PCAC. Making use of the fact that the moment of inertia 
7(6) is a functional of the profile function 9 they solved this new equation together with ( ]9"T| ) 
iteratively and found a new selfconsistent solution which satisfies PCAC. In fact, in their calculations 
following the incomplete treatment of Wakamatsu and Watabe (1993), they disregarded the Dirac 
sea contribution. The full calculations, however, including the Dirac sea as well, show no solution, 
which is not surprising, since this new "equation of motion" does not correspond to any stationary 
point of the effective action. 

Apparently, as argued above the present scheme based essentially on the large N c approximation is 
not very appropriate for evaluation of quantities like the divergence of the axial current which requires 
the use of the equations of motion. Hence, on theoretical grounds, the estimate gA calculated directly 
from the matrix element of the space component of the axial current is preferred to g\. Accordingly, 
the numbers of ng.[T2] correspond to gA and not to g\. 

On the other hand, it is not an easy task to generalize the present large N c treatment in order to 
include the zero-mode rotational corrections in the saddle-point equation. It would correspond to 
projection before variation and a deviation from the hedgehog structure already for the saddle-point 
solution would be inevitable. It would also lead to conceptual problems related to the consistency 
of 1/N C expansion which we have used so far. Therefore, the most general solution would be to go 
consistently beyond the large N c approximation including the quantum meson (loop) effects from the 
very beginning. 



Closely related to the PCAC is the Goldberger-Treiman (GT) relation. It makes an important link 
between the strong g^NN an d the axial vector g& coupling constants and the pion decay constant f w . 
In fact, this relation is valid in nature within some percents. Let us assume that the PCAC equality 
( |244| ) is valid and use fl246|) to evaluate gA- Inserting it in fl241| ) one gets 



9a = • (251) 

M N 

which differs from the GT relation in ^jvjv(O) taken at q 2 = and not at the pole q 2 = m 2 . In fact, 
in the model g^ differs from g\ by 15 % which means that with the estimate g& the relation (|251D is 
fulfilled at the 15% level and apparently, due to the approximations used, the present model scheme 
cannot describe the Goldberger-Treiman discrepancy observed in nature. 



4.6 Tensor charges 



The nucleon tensor charges gP (a = 0, 3 corresponds to the singlet and isovector charges, respectively) 
are defined as the nucleon forward (zero momentum transfer) matrix element of the tensor current: 

(N(p)|*(Vy*|JV(p)) = gPu(p)(T^u(p) , (252) 

where = |[7 M ,7^] and t(0) = 1. Although the tensor charge like other nucleon charges is 
a fundamental quantity, which characterizes the nucleon, too little is known about its values and 
its implication to the nucleon structure. The reason is that the tensor charge is difficult to access 
experimentally, since there is no experimental probes coupled directly to the tensor current. However, 
as it is shown by Jaffe and Ji, (1991) the nucleon tensor charge is related to the first moment of the 
transversity quark distribution h\{x): 

i 

Jdx(h[ a \x)-h{ a \x))=gP, (253) 
o 

where a = 0, 3 stands for the singlet and isovector transversity quark distributions, respectively. This 
sum rule gives a hope to gain some experimental information about the tensor charge. It should 
also be noted that the tensor charge depends on the renormalization-scale dependent. However, the 
dependence on the normalization point is rather weak 
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/ , . . 2 \ \ 

,(/)/ 



# « )= \M) # ( "° ) ' (254) 

which means that the initial value of the normalization point is not of big importance. 



In the non-relativistic limit using the identity = £ikjloljl5 one can easily see that the nucleon 
matrix element of the tensor current has a structure similar to those of the axial current: 

gP = "FT E R t(^, €m) < n\r a \m >< m\ lo a a \n >= aN° c , (255) 
9t ) = N c < n l7o<?" r\n > E ~Z — < m U[^ x r] a \n >< n\r a \m >= f3N c + 5N° C 

e n <val 9 ^ c n >val e n e m 

€m<val 

(256) 

In contrast to the case of the axial current, here the singlet tensor charge diverges and needs regu- 
larization: 

2 2 

K(e n , O =~ ^ ~ \ f , (257) 

whereas the isovector tensor charge is finite and is not regularized. In large N c limit g^ ~ N® and 
gj, ~ N c , which is the same as in case of the non-relativistic quark model. In fact, the isovector 
charge, similar to the axial vector coupling constant, contains non-zero contributions from both the 
leading (N c ) and next to leading terms (N®) and the consistency condition of Dashen and Manohar 
(1993) is also fulfilled. 

The model prediction for the nucleon tensor charges for M = 420 MeV (Kim et al, 1995d): 

gP w 1.45, ^ 0) ^ 0.69, (258) 

or 

gP w 1.07, # w -0.38 (259) 



are close to those in bag model (Jaffe and Ji, 1991) and not so far (within the error bars) from the 
QCD sum rule estimates of (He and Ji, 1995). 



4.7 Electric polarizability of the nucleon 



So far we have considered nucleon observables which are related to matrix element of a single quark 
current. In this subsection, we discuss a more complicated case of an physical quantity which is 
related to a matrix element of two quark currents, namely the nucleon electric polarizability. 

Recent measurements of the electric, a, and magnetic, (3, polarizabilities of the nucleon (Federspiel 
et al, 1991, Zieger et al, 1992, Schmiedmayer et al, 1991) narrowed considerably the experimental 
uncertainties in these observables, and were accompanied by a number of theoretical studies. 

The electric polarizability is given by the change of the nucleon energy in an external electric field 
and it can be expressed in terms of the correlator of two nucleon currents (|94]) . Apart from the 
regularization the explicit calculations give for the leading- N c contribution to a 

Q <») = VeWMl") £ Hz^HgW (260) 

where \N) is the collective nucleon state, is an element of a Wigner D-matrix, and \n) and 
e n denote the eigenstates and eigenvalues of the Dirac hamiltonian. This formula can be split into 
the valence and sea parts. The sea contribution is divergent and includes a cut-off function. The 
vacuum contribution has to be subtracted as well. Finally, evaluating the collective matrix element 
one obtains the leading- A^ c valence and sea contributions to the electric polarizability of the nucleon 
in the form 



a 



(o) N c 2 ^ (val | r a z | n)(n | r a z | val) 



val 



^ g2 ^ tlim^l tl^L ; (261) 



n;;:.; = — ^ e 2 V TZj(e m ,e n )(m\T a z\n)(n\r a z\m)- vac, (262) 





'■sea 



where lZf(e m ,e n ) is the proper-time regularization function of the moment of inertia / ( |138|) . The 
subtracted vacuum contribution in ( |262| ) has the same structure as the first term, but the free quark 
states and eigenvalues are used. 

The rotational 1/N C correction a*- 1 ^ is given by a triple sum over quark states (Nikolov et al, 1994). 
The valence part a^j has been calculated exactly, whereas the sea part a^l has been approximated 
by the leading term in the gradient expansion ageti lowest . 

In fig. ^ we present the results (Nikolov et al, 1994) for the average electric polarizability of the 
proton and neutron, a = |(a p + a n ) obtained in the NJL model as a function of the constituent quark 
mass M. The different contributions as well as the total nucleon electric polarizability are presented. 
The long-dashed line corresponds to a' = + a^j + a^ a lowcst + a£j and this is the direct result of 
the NJL model. The solid line shows the total a, roughly corrected for the effects of the A^-A mass 
splitting. For values of M in the physically relevant range 400-450 MeV the changes of a with M 
are small. The present results are in qualitative agreement with the results of other models (Weiner 
and Weise, 1985, Broniowski and Cohen, 1993) with explicit pion degrees of freedom. 
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Figure 17: Electric polarizability of the nucleoli of the SU(2) NJL model (Nikolov et al, 1994) as a 
function of the constituent quark mass M. The a' indicates the direct result of the NJL model, the 
a incorporates corrections due to N — A mass splitting. The polarizability of the nucleon is defined 
as a = (a p + a n )/2. 

The sea contributions clearly dominate over the valence ones both on leading order in N c and for the 
rotational 1/N C corrections. These contributions are very well reproduced by the leading two terms 
in the gradient expansion and the first term is the so called "seagull" contribution to a, discussed in 
many previous papers. 

The inclusion of rotational corrections has important phenomenological consequences. The dominant 
contribution to the electric polarizability is obtained from pion tail effects, and is proportional to 
the axial constant g\. Since the rotational corrections have been shown (Christov et al., 1994) to 
be crucial for reproducing the experimental value of qa in the model, it is important to take into 
account these rotational corrections for the electric polarizability. The calculation shows that indeed 
the rotational corrections give a sizeable contribution of about 50-60% of the leading- N c result. After 
including approximately the corrections due to the N-A mass splitting, the theoretical value is closer 
to experiment than in other studies in soliton models. For the typical choice of the constituent 
quark mass M = 420 MeV we obtain a' ~ 25 x 10~ 4 fm 3 and a ~ 19 x 10~ 4 fm 3 , compared to the 
experimental value a exp = 9.6 ± 1.8 ± 2.2 x 10~ 4 fm 3 (Federspiel et al, 1991, Zieger et al, 1992, 
Schmiedmayer et al, 1991). The model prediction is, however, still too large indicating that the 
inclusion of other 1/N C effects (such as e.g. meson loops) is necessary. In other approaches like e.g. 
the chiral perturbation theory with nucleons (for review see (Bernard et al, 1995)) the pion loops 
provide the main effect. In our model, parts of those meson effects are taken into account already in 
the leading 1/N C order due to the presence of the mean-field. The effects of meson loops in the next to 
leading order in 1/N C remains to be investigated. In particular, one hopes that the meson-loop 1/N C 
effects related to the zero modes would lead to a reduction in the case of the electric polarizability. 



5 SU(3)-NJL model 



In the previous chapters, we have reviewed the Nambu-Jona-Lasinio model with scalar isoscalar 
and pseudoscalar isovector couplings. The extension of this model to the larger symmetry of SU(3) 
introduces not only the kaon and eta mesons to the model but offers in the baryon sector, which is 
the main concern of the present article, the possibility to calculate the properties of strange baryons 
as well as the contributions of strange quarks to the properties of the nucleon and Delta. In fact, at 
present these effects are in the case of the spin structure functions or the strange vector form factors 
under vivid discussion. 

We start with a very short review of the mesonic sector which is already extensively discussed in 
the reviews of Klevansky (1992), Vogl and Weise (1991) and Hatsuda and Kunihiro (1994). After 
performing the bosonization and saddle-point approximation we fix the parameters of the model by 
reproducing pion and kaon properties and have again the constituent quark mass of the non-strange 
quarks as the only free parameter. Actually, the extension from SU(2) to SU(3) is formally not trivial, 
since the current mass of the strange quark, m s , is noticeably larger than that of the up and down 
quarks, m u and m^. This corresponds to a non-negligible explicit breaking of the SU(3)-symmetry 
which plays an important role in the baryon sector. We remind that the semiclassical treatment of 
the NJL model is justified in the large N c limit. In the SU(3) case, apart from the 1/N C expansion 
one has to deal also with the expansion in powers of m s . These expansions in two small parameters 
1/N C and m s are very sensitive to which limit is being taken first: N c — > oo or m s — > 0. This leads to 
certain difficulties in the collective quantization of the soliton rotation. We discuss the perturbative 
approaches (Weigel et al, 1992, Blotz et al, 1992, 1993a), which in view of explicit calculations of 
baryonic observables are most far developed at the present time. Special emphasis then is put on 
the possible representations of the collective hamiltonian. After a discussion of the mass spectra 
of the hyperon multiplets in SU(3) we review baryonic observables like the electromagnetic form 
factors, sigma terms, axial vector couplings, spin properties of the proton and the Gottfried sum rule. 
Predictions will be done for the nucleon strange form factors. 

Besides the perturbative approach one can also consider the hyperons as bound states of kaon mesons 
in the background field of the SU(2) soliton. This approach has been invented by the Callan and 
Klebanov (1985) and has been applied in the chiral sigma model by McGovern and Birse (1990) and 
in the case of the NJL model by Weigel et al. (1994). 



5.1 Extension of the NJL model to SU(3) flavors 



An extension of the SU(2) flavor NJL lagrangian (1) to SU(3) can be done by a direct generalization 
(r a — > A a ) of the SU(2) case with scalar and pseudoscalar couplings. Then the simplest SU(3) form 
of a four-fermion interaction with scalar and pseudoscalar couplings reads 



G 



a=0 



C = V(x)[ip-m]V{x) + — Y i ($>(x)\ a m(x)) + (y(x)\ a i^(x)y 



where the A a are the Gell-Mann matrices with SpA a A 6 = 25 a b and the sum over the A a matrices includes 
A = y/|1 3 as well. This form of the couplings is demanded by requiring a chiral SU(3)r <g> SU(3)l 



symmetry. The mass matrix in ([5]) is now enlarged by the strange quark mass 

/ m u \ 
rh = rrid = rh l + A 3 m 3 + A 8 m 8 




where m = ( m u + m d + m s)/3, m 3 = (m u — m d )/2 and m 8 = (m„ + m rf — 2m s )/(2y / 3). Note that 
in contrast to SU(2), where one could write a minimal lagrangian with the singlet scalar field only 
and a triplet of pion fields, the rank 2 of the SU(3) group compared with rank 1 of SU(2) together 
with chiral symmetry requires the presence of the whole nonet of scalar and pseudoscalar couplings. 
Using integral identities similar to eq. (15) the bosonized form of @ reads 

8 



1 



C! = f (x) [ip - m - V (a a X a + i l5 7r a \ a ) V(x) --Ek« + ^ ~ 



a=0 / 2Lt a=0 



7T„ 



Here we have added "by hand" an extra term C' A [a a ,7i a ]. The necessity of this additional term is 
related to the Ua(1) problem. Besides SU(3)r® SU(3)l<& Uy(l) symmetry, the lagrangian C (||) has 
in addition the Ua(1) symmetry. On the other hand, it is well know that in QCD the U^(l) symmetry 
is broken by the anomaly. As a consequences the rj' meson remains massive even in the chiral limit. 
This means the effective low energy theory should contain some term C A [a a ,ir a ] explicitly breaking 
the Ujl(1) symmetry. In fact, we use the presence of C! A as an argument to consider rf massive and 
not as a Goldstone boson. For the fixing of the parameters of the NJL model only Goldstone bosons 
are needed. In the baryonic sector, we will use the chiral circle (||) and hence the detailed structure 
of C' A is irrelevant for our considerations. 

After integrating over the quark fields the partition function in Euclidean space is given by 

Z = J Va a Vn a e- s ^ 

where 

S(a, tt) = -iV c Tr log Z^, tt) + ^ J d 4 x £ (a 2 a + ^) + j d 4 xC! A \a a , vr a ] 
with D(a, it) = d T + h and 



h = -ij j k d k + 7 m + 7 ^ (a a X a + z7 5 7r a A°) ■ 

a=0 



5.2 Vacuum solutions and fixing of parameters 



In the vacuum we make again a saddle point approximation similar to eq. (28) 



8S(a, tt) 

8<?a 







SS(a, tt) 

fin a 







For three flavors there are the <To, 03 and erg fields which can have non- vanishing vacuum expectation 
values. Therefore, instead of eq. flH| ) one has in general one gap equation for each flavor degree of 
freedom. Neglecting isospin breaking (m = m u = m^) in the following the non-trivial gap equations 
read 



a, - 8N c axh (ax ) 
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G 

—a 2 - 8N c a 2 Ii((T2) 



mo 

G 

m s 

G 



(263) 



where we set o\ = y |o"oc + ^ a sc + m o an d o"2 = y f^Oc — T^sc + "V Deriving these equations we 
neglected the contribution of the C A [a a , 7r a ]. 

The spontaneous breaking of chiral symmetry leads now to the appearance of constituent quark 
masses for the non-strange and strange sector according to 

M u = M d = a u M s = a 2 . 

This resembles two independent Dirac equations with mass gaps M u and M s . In addition, there are 
two quark condensates 

A A 1 1 9 <7 

< uu > = < da > = —— —— — Zj 

V4 2 om 

1 d 

<ss> = -tt^—Z (264) 

which become equal in the chiral limit and actually are of the order of — (250 MeV) 3 . Using the gap 
equations eqs. (|263|) the eqs. (|264|) can be cast into the form 

< uu + dd> = — — (M u — m ) 

< ss > = -— (M s - m s ) . (265) 

Similarly to the definition of the pion decay constant in eq. ([70]), we can define the corresponding 
vacuum matrix element of the axial vector current A a ^ a = 4, 5, 6, 7 with the physical kaon fields 
7r a , a = 4,5, 6, 7. As a result one obtains (Schneider et al, 1995) 



f n = M u g n Z p (q 2 = -m 2 n ) 

M s + M u 2 2 

Ik = ^ 9KZ p {q = -m K ) (266) 



where gx can be determined from the corresponding eq. fl6"l"|). In the chiral limit, from the ( [265| ) 
and the gap equation (|263|) it follows that 

U = f K = M u g^Z p (q 2 = 0) 

and 

m\ _ m s + m 

For finite quark current masses the kaon decay constant deviates from the pion decay constant but 
however it is still underestimated by 10 %. 

The masses of the pion and kaon are determined from the two-point meson functions (50) 

2 m 



GM u Z p (q 2 = -ml) 
m s m \ 1_ 



with A A , = M, - M„. 



Fixing of the parameters. The SU(3) NJL model in the present form contains four unknown param- 
eters: the coupling constant G, the current quark masses mo and m s , and the regularization cut-off 



A. First, instead of G the constituent quark mass M u is used as parameter which is related to G 
by the first gap equation ( |263| ). Second, for a given M u the m , m s and A are fixed to reproduce 
f n = 93 MeV, = 139 MeV and the kaon mass rrix = 496 MeV. The M u is not fixed in the meson 
sector but in the baryon sector, in order to reproduce best the experimental data of the nucleon. This 
results, as in SU(2), in preferred values of M u about 420 MeV. The strange constituent mass M s is 
not an independent quantity but determined by M u according to the second gap equation ( j263j ). 

In contrast to the SU(2) case, where the average current quark mass tuq is rather small and the results 
do not depend on its particular value fixed by the physical pion mass, in the SU(3) case, the explicit 
chiral symmetry breaking due to the strange current mass m s « 175 MeV (Bijnens et al, 1995)Qwill 
be treated perturbatively up to the second order corrections in m s . Apparently we need to reproduce 
the value of the strange quark current mass properly. Actually, if one follows the above scheme and 
uses the proper time regularization (p|) with <p{u) = 6(j^ — u), the absolute values of m u , rrid and m s 
come out about two times larger than expected. To that end we make use essentially of the freedom 
to choose the particular form of the function 4>(u) in the proper-time regulator ([41]) in order to adjust 
the quark condensate and the current quark masses to their phenomenologically accepted values. In 
our practical calculations we use 

(j){u) = c6(u - Af) + (1 - c)9(u - Af) (268) 

in which two additional parameters c and A 2 appear. The c is used to obtain a value for ttiq of 6.1 
MeV. The A 2 is chosen as small as possible in order to facilitate the numerical effort in the calculation 
of the soliton. 



5.3 Restriction to Goldstone modes 



Up to now we have considered the meson sector of SU(3) NJL in terms of independent fields a a and 
ii a with a = 0,1,. ..,8. The parameters of the model are fixed by pion and kaon properties. Like in 
the SU(2) case the model as such does not support solitonic solutions and, again, one has to subject 
the field to the condition of the chiral circle. In analogy to ([76|) this is defined as 

a a (x)X a + iTr a (x)X a = MUf\x) = Me< A<1 

where (x)U^(x) = 1 similar to the case of SU(2) chiral fields. The sum goes over a = 0, 1, . . . , 8 
and M is given by M = M u even in the case of M s ^ M u . This is consistent with the so called trivial 
embedding of SU(2) into SU(3) discussed in the next subsection. 

The condition of the chiral circle corresponds to a clear separation between the Goldstone degrees 
of freedom, i.e. pions, kaons and eta, and non-Goldstone ones like sigmas. Since we are interested 
in describing the ground state properties of the nucleon and the hyperons, the degrees of freedom 
corresponding to heavy mesons are of less importance and hence, will be ignored. With the U^\x) 
of eq.(|l|) in the solitonic sector, the present model is identical to the SU(3) version of the Chiral 
quark soliton model of Diakonov et al. (1988). This model is based on the instanton model of the 
QCD vacuum (Diakonov and Petrov, 1986). We remind that the effective quark lagrangian derived 
by Diakonov and Petrov (1986) from the instanton model of the QCD vacuum has a form of a 't 
Hooft-like interaction which explicitly breaks the Ua(1) symmetry. Its bosonization leads to meson 
fields of the form (JJ). 

1 We note that for the strange current mass we use the value at the same MS scale /j, = 1 GeV (of order of the model 
cutoff) as in the case of the quark condensate. 



5.4 Trivial embedding of SU(2) hedgehog soliton into SU(3) 



By analogy with the SU(2) case we expect that the low-lying baryons can be described as rotational 
excitations of some classical soliton solution. The spectrum of these baryons in the real world imposes 
certain restrictions on the symmetry of the SU(3) soliton. It was suggested by Witten (1983) that 
one should consider the trivial embedding of the SU(2) hedgehog soliton into SU(3) 



Uf\x) = 














where U^>{x) = cosG(r) + irxsmQ(r) coincides with the SU(2) stationary meson field configuration 
U c (x) of eq. (|116|) . The advantage of this ansatz for the saddle point solution is that its rotational 
excitations naturally describe the particles of the SU(3) octet and decuplet. 

If one neglects the quark masses then there is a continuous set of saddle point solutions R^U^> (x)R 
where R is an arbitrary SU(3) matrix. Similarly to the SU(2) case, in the path integral scheme we 
have to take into account the path integral over time dependent fields of the form 

U^\r,x) = R{T)Uf\x)R\r). 

After this path integral is computed one arrives at the baryon rotational wave functions ip(R) de- 
pending on the SU(3) orientation matrix R. As in the SU(2) case the SU(3) flavor generators and the 
spin operator acting in the space of this rotational wave functions are nothing else but the generators 
of the left and right rotations of the matrix R. By analogy with ( |145| ), ( |146| ) we obtain 

exjp(iu a T a )4>(R) = 4> (exp(-iu a \ a /2)R) , (269) 

exp(iuo a J a )ij(R) = V (Rexv(iuj a \ a /2)) . (270) 

Here T a are eight SU(3) flavor generators. As for the right generators J a , the components with 
a = 1, 2, 3 have the meaning of the spin operator. The other components of J a have no direct physical 
meaning. Note that under the rotation associated with Jg the field (r, x) does not change. This 
means that the parametrization of the rotating soliton by the SU(3) matrix R contains "non-physical" 
degrees of freedom. As it will be seen later this leads to a constraint on the allowed eigenvalues of J§. 



5.5 Expansion in angular velocity 



Restricting the path integral over U(x) in ((0) to meson configurations of the form U^(r, x) ([J) we 
obtain by analogy with ( |128| ) 



i=i 



where the V ^ carry the quantum numbers of the hyperons and the Dirac operator in the body-fixed 
frame of the soliton is now given by: 

D{UP (r, x)) = R[d T + h + h 1 + ttt] R ] . 

Here ho is the single particle hamiltonian 

h = -iiol k d k + M ulo Uf^ +7om l ( 271 ) 



which includes the two flavor hamiltonian eq. fl79|) embedded in SU(3). In contrast to SU(2) case in 
(|I|) we have additionally the term 

which appears due to the fact, that the mass matrix rh does not commute with the rotation matrix 
R(t). We also have defined 

A Ma =M S -M U . 

Below we consider hi as a perturbation. Note that there is a freedom of the redistribution of the 
singlet part of the strange quark mass between Hq and h\. Different saddle points corresponding to 
this freedom have been studied by Schneider et al. (1995). In the following we will discuss the case 
where the saddle point is determined by using hamiltonian (|271|) and h\ is treated as a perturbation. 



Expanding the correlation function (|l|) in powers of Q and AM a (m s ) up to the second order in the 
rotational velocity we get for the propagators 



N c l 



TT(0,T/2 I M I 0,-T/2) fi9 , 

exp \-N c e vaX {U)T - \Q l J drQ a Q b - J drQ 8 - J drL ml [A Ms ]j (272) 



as well as from the effective action 

exp(-S(RU^rf)) = exp (-TE sea (U) - i/ a T / drQ a Q b - J drL sea [A Ms 

where I[A^J = L sea [A Ms ] + L va i[A Mg ] = 0(A Mg ) and will be discussed below. The tensor of the 
moments of inertia I^b is diagonal but in contrast to the SU(2) case its components are different 

r J 1; a = b= 1,2,3 
lab = I h, a = b = 4,5,6,7 
[ 0, otherwise. 

The second important difference from the SU(2) case is the linear term in Q 8 in ( p72j ). In Skyrme 
type models, this expression is proportional to the topological winding number of the chiral field 
which in these models coincides with the baryon number. In the present non-topological chiral model 
the term with Q 8 is due to the discrete valence level in the Dirac spectrum. The presence of the Q 8 
term is a reflection of the fact that right rotations associated with the generator J 8 have no effect on 
the chiral field. Therefore, the generator J 8 corresponding the rotational velocity Q 8 is constrained. 
From the term linear in Q s in eq. ( 272|) we see that 



■h 



N, 



c 



2^ ' 

By analogy with the hypercharge 

Y=(2/V3)T S , 

one can introduce the "right" hypercharge 

Y R =(2/V3)J 8 , 

in terms of which the constraint ([]]) can be rewritten in the form 



Y R = -N c /3 



This in turn transforms into a constraint for the possible representations of SU(3). Actually, one 
can show that only representations with zero triality survive. These are the octet and decuplet 
representations with spin 1/2 and 3/2 respectively. 

By analogy with the case of the SU(2) case the effective rotational lagrangian leads to the following 
quantization rules 

f iWh, a = 1,2,3 

J a = \ tV a I 2 , a = 4,5,6,7 

Here, the components J a with a = 1,2,3 correspond to the spin operator. Using these quantization 
rules we obtain a collective hamiltonian from ( |272[ ) and ([TJ) in the form 

H ^-\{{-Ty^ u ^ + k Asm) -h 

Here 

C 2 (5f/(3)) = E^ 2 = ET a 2 (273) 

a=l a=l 

is the quadratic Casimir operator of the SU(3) group and 

C 2 {SU{2)) = j:Jl (274) 

a=l 

is the Casimir operator of the SU(2) subgroup corresponding to spin generators J a (a = 1, 2, 3). 

r(0) 



The upper index (0) in H^J U indicates that we are working in the chiral limit without any mass 

I coll aI1 d H coll- 



corrections. Those will be incorporated later and will yield additional terms and 



The Hamiltonian H^ u commutes with all flavor SU(3) generators T a . Therefore the eigenstates of 
H^ u can be classified according to the irreducible representations of SU(3). Below we use the well 
known (p,q) parametrization the SU(3) representations in which the octet is (1, 1) and the decuplet 
is (3,0). 

Since in our model the spin J a and flavor T a operators are realized as the generators of right and 
left rotations of the same orientation matrix R, the representation (p, q) should be the same for T a 
and for J a . Apart from the representation (p, q) the eigenstates of H^ u are also characterized by the 
quantum numbers corresponding to the following set of commuting operators: 



2 



3 



0=1 

3 



2 

7T 



Y R = ^=J 8 , Y,Jl Js- (275) 



a=l 



The corresponding wave functions are 



^ q T,T 3)i y R ,j,j s) (R) = (-l) yi?/2+J3 ^( ( rT;T3)(-y„.,-,3)(^) (276) 

where the asterisk stands for the complex conjugation and dim(p, q) is the dimension of the represen- 
tation (p, q). 

With these wave functions the Casimir operator has the following spectrum 

C 2 (SU(3))^ T3;YM = hp 2 + q 2 + 3(p + q) + pq)^ T:j){YM) , (277) 



C 2 {SU{2))^% T , i){YM) = J{J + 1)^, T3)(1W , J3) • (278) 

In the chiral limit presently considered, the members of each multiplet are degenerate in energy. Using 
( |277| ), ( |278| ) we find the eigenvalues of the rotational hamiltonian (]T|) which leads to the following 
result for the energy of the center of the multiplets, defined as the average mass of all the members, 

E p , q>J ■ = M d + \(±- £) ^ + l) + ^(p 2 + ^ 2 + 3( P + g)+pg)-^-. 

Note that constraint ([[]) fixes the spin of the particles of the (p, q) multiplet. The analog of the 
nucleon delta splitting (|156|) has to be interpreted in SU(3) as octet and decuplet-splitting, i.e. the 



energy difference between the centers of the octet and decuplet, and follows as 

-^8-10 = -^3,0,3/2 — £-1,1,1/2 = 77^- • 

All 

Formally this coincides with E^k in SU(2), however, the phenomenological value for the difference 
of the average octet and decuplet masses is now El\ = 230MeV. 



5.6 Strange mass terms of the collective lagrangian 



In order to calculate the splitting within the multiplets we have to include the SU(3) symmetry 
breaking terms due to hi in ([!]). Note that we deal with two quantities which may be considered as 
small quantities. The first one is the rotational velocity Q, which is of the order of 1/N C . The second 
one is the mass difference A^ s = M s — M u . Assuming both corrections to be roughly of the same 
order we use systematic expansions up to the second order. 

Expanding the baryon correlator ([!]) and considering sea and valence contributions together, one 
obtains for L[A Ms \ in fl27|) and (@) 

L[A Ms ] = Li[A Ma ] + L 2 [A Ms ,Q) 

+L 3 [A 2 Ms }+0(AZ s n n ) m + n> 3. (279) 

The first term in this expansion 

ii[AM.] = -%(l-BS»Wr))) 

777-0 O ' 

contains the SU(2) sigma term a, which in the present model is given by 

dM cl [m ] 



m - 



drriQ 

Then we have in the next order O (QAm s ) 



mo=0 



L 2 [A Ms ,n} = %2A^ 1 ^Dg ) ( J R(r))a(r) 
V3 i= i 

A-- 7 

4 



-^-2K 2 J2D^(R(r))n a (r), (280) 

V3 a=4 



where Ki and K 2 similar to 7i , I 2 are related to tensor K a b 

Ki, a = 6= 1,2,3 
K llh { K 2 , a = 6 = 4,5,6,7 
0, otherwise. 



The latter can be split in valence 



Kvai _ N c (vol \ X a \n)(n\ AS | val) 
ab ~ 2 % e n -e val 



and sea part 



r^sea _ N c (m | A a | n)(n | A b 7o I m) 
ab ~ 4 ^ e -e 



sm<0 



In contrast to the sea contributions of I\ and I2 that of _K" a b originates from the imaginary part of 
the effective euclidean action. Hence, it is ultraviolet finite and needs no regularization. One should 
note that in Skyrme approaches with vector mesons such contributions from the effective action are 
obtained only by adding the gauged Wess-Zumino term (Park and Weigel, 1991). 



From the quadratic corrections in the symmetry breaking we obtain 
Ls[A 2 Ms ] = (^\ Nl j:Di\R(r))DS\R(r)) 



+ 



+ 



'A 



Ms 



7 



2N 2 J2D®(R(t))DW(R(t)) 



(8), 



a=4 



'A 



M s 



where N% and N 2 are diagonal elements of tensor N a b 



-N (l-D<g{R(T)j) (1-D$>{R(t)) 



ab 




a = b = 1,2,3 
a = b = 4,5,6,7 
a = b = 8. 



It consists of a valence part 



ab 

and a regularized sea contribution 

i\Tsea c - 

Iy ab — ~T 



AT vai _ N c (val I A a 7o \n){n\ A fe 7o I val) 



n^val 



E n — E, 



val 



J2( m I -^ a 7o I n)(n I A 6 7o I m)Tl%{e n ,e r 



n,m 



The regularization function 1Z% is given by 



du 



2^ Jo \fu 



cf>(u) 



(281) 



Before we can study the mass splittings from the collective hamiltonian resulting from (|281 ) , one 
has to remember that in the presence of the linear terms in the rotational velocity Q, the generators 
of the group are modified according to 



Jn 



0»I 1 -2^K l Dg>(R(T)), 
WI 2 -2^K 2 dS(R(t)), 



2V3 



V3 V 
~2~ Ri 



a 
a 
a 



1,2,3 
4,5,6,7 



5.7 Hyperon splittings in linear order - Sum rules 



Here, we consider the hyperon splitting of the octet and decuplet in linear order in A« s and neglect 
the isospin breaking, i.e. set m u = = m^. As shown in the meson sector in linear order of Aa/ s 
the difference M s — M u equals m s . Now the collective Hamiltonian 

ticoii - n coll + n coU 

includes a new term 

o m s / 2m., (K] Kn 



-m s —Y + — (R) . (282) 

The moments of inertia and the matrix elements of the Wigner functions can be found in (Blotz et 
al, 1993a). 

Although H^ n contains three different operators Y,D^(R) and X)f=i Df^(R)J i) the mass splitting 
can be parametrized in terms the two quantities 



3m V h h 

and 

6 h vP = m s ^. 

h 

For the masses of the octet and decuplet baryons, calculated relative to the mass of the £*, one gets 
(Blotz et al, 1993a) 



Am N = A h y p - 6 hyp -£ 8 -io 

Am A = -i A*«p -£ 8 _ 10 

Am s = i A a «p -£ 8 -io 

Am H = | A^ p + <^ --&8-10 

Am A = -| A^ p - <f h w 

Arris* = 

Am H , = | A h yp + 5 hyp 

Am Q = I A^p + 25 hyp 

The reason that we consider the mass spectrum relative to the £* is that the mass of the S* almost 
corresponds to the center of the mass of the decuplet. From the formulas ([!]) above, it follows the 
mass relation of Guadagnini (1984) 

1 

m=* — ms* + rax = -(11mA — 3ms) • 
8 

Here it is interesting to note that this sum rule can be obtained in the pseudoscalar Skyrme model 
only by introducing the hypercharge operator by hand with some unknown coefficients. In the present 
approach, this term arises naturally ( |282| ). 

Furthermore, one obtains also the Gell-Mann - Okubo relations 

3 11. . 

-m A + -m s = -{m N + m H ) 

and 

rriQ — m=. = m=* — ms* = m^* — m A . 

It should be noted that these relations are mass sum rules which rely basically on the fact, that 
the SU(3) flavor symmetry breaking part of the strong interaction can be treated in first order 
perturbation theory. Their validity shows in some way the consistency of the quantization procedure. 
On the other hand the values of the hyperon masses themselves depend on the dynamics of the present 
model and this goes beyond the perturbative sum rule approach. 



5.8 Hyperon splitting in second order - Wave functions 



In the present scheme, the corrections in second order of the strange quark mass have three different 
origins. The first one results from the expansion of the baryon correlator in Am s = M s — M u . The 
second one is due to the quantization formulas ([!]) , which replace the angular velocities in ( [28 1| ) 
by the corresponding generators. Accordingly, the collective hamiltonian contains not only 

«S = ^(l-DffW) 

V3 i\ i=1 V3 h a=4 

but also second order terms in the AM s (m s ) 

H { 1 = 2^^±D£\R)D£\R) + 2^^±DiS(R)D^(R) 



-2N^ £ D$>(R)D$>(R) - 2N 2 ^k £ D g>{R)D§>W 

6 i=l 6 a=4 

-2N Q ^(l-D<8{R)Y . (284) 



To study the third contribution in second order of the symmetry breaking we recall the perturbation 
theory for our system of H^ u + H^ u + H^ u . For the energy of a baryon B in the representation [n], 
which has to be the lowest dimensional one where the quantum numbers of the baryon appears, we 
can write in increasing order of the perturbation: 

4 0) n = (B,n\H^ u \B,n) (285) 
E$. = (B,n\H^ u \B,n) (286) 



E {2) 



m r, I FJM i RnU v I (B,n' | I B,n) | 2 
( B , n | B coll \B,n)+ ^ - Q) . (287) 

n^ri ^B,n ^B,ri 

Here, ^ is the energy of baryon of quantum numbers B in the representation [n] with strangeness 



corrections of the order m taken into account. In the lowest order the energies Eg n (|285 ) are degen- 
erate for each multiplet and correspond to the mass center of the multiplet. In the next order Eg n 
include linear strange mass contributions from the matrix elements of H^ u between the unperturbed 
wave functions | B,n). This removes the degeneracy for the zero-order results Eg n . The second-order 

contributions E B n ( |287| ) are from H^ oll between the wave functions of lowest order as well as from 
the last term in ( 287|) . There the summation is over all representations \n'\ ^ [n]. In fact, this 



term appears since the baryon wave-functions of the collective hamiltonian H^ u + H^ a after diago- 
nalization in the space of all | B, n) are no longer pure octet or decuplet states. Because the operator 
Hcoii commutes with Y, T, T3 and J operators, there are in addition only higher representations like 
[10] of [27] in the case of octet baryons and [27] and [35] in the case of decuplet ones. These higher 



representations enter the baryon wave functions implicit in Q287| ) perturbatively in first order of the 
symmetry breaking and therefore the total contribution to the energy in (|287|) is quadratic in the 
symmetry breaking. 

Actually, if the baryon is a member of the octet, the only non- vanishing matrix elements can be 
parametrized by the quantum numbers of the baryons according to 

I (B, 10 I if« I B, 8) | 2 = 1 (a + 7 |) 2 <Y + \y 2 + T(T + 1) 



and 



(B, 27 I fl« I B,8> |»= -If „ - 7 iW 9 + jj (V(T + 1) - IyA - 7 -Y'< 



100 V 6/ V 2 V v '4/4 



with the abbreviations 



a = 


1A M 

3 m 


K 2 A 

- ~t^m s 




(3 = 


A K2 












7 = 




h) ' 


(288) 



Similarly, if the baryon is in the decuplet, the only non-vanishing matrix elements are given by 

( B ,27|7/S,| B ,10) = i( a + 7 5) 2 ( 1 + 3 r+ l y2 



and 



(B, 35 I I B, 10) = A (« _ i 7 ) 2 (i-\y- ir 2 ) . 



We will consider expectation values of external current operators, for which we need also the wave 
function | B) itself in the next to leading order of the symmetry breaking. Then the wave functions 
can be written in first order perturbation theory as 

\B) -\n,B)+ ^ \ n, B) . 

Furthermore, the wave functions for the octet baryons, which will be the case of special interest in 
the next sections, in linear approximation in Am s read 

\b) = |fi,8)-^=( a + i 7 )(r + ^ 2 + r(r + i)) |s,fo> 

~ (« - ^y) (Ver 2 + 3(1 - 7 -y") - ±t(t + 1)) | b, 27) (289) 

The expectation values of these current operators, corresponding to some observable quantity, can be 
obtained consistently by sandwiching the operators, calculated up to the linear order in the symmetry 
breaking, between the unperturbed wave functions ( ( f276| ) ) and then adding the matrix element of 
the operators, calculated now in zeroth order of the symmetry breaking, between the perturbed wave 
function (( |289| )). In this scheme, one can calculate consistently the same type of corrections from 
the symmetry breaking to the expectation values of these current operators as it was previously in 
the case of the energy. However, because of technical reasons, current expectation values are always 
calculated with linear rotational and linear strangeness corrections. 



5.9 Yabu-Ando diagonalization method 



In contrast to the former perturbative methods of evaluating the collective hamiltonian either between 
the unperturbed eigenfunctions ( [276| ) or the first order perturbated ones fl289| ), in the case of the 



wave-function corrections, Yabu and Ando (1988) developed a method for treating the baryonic wave 
functions to all orders in m s or in our scheme in Am s = M s — M u . Without reviewing technical 
details, which can be found in the original work (Yabu and Ando, 1988), the method is based on 
diagonalizing the collective hamiltonian H^ u + H^ u + H^ u in the basis, spanned by the eigenfunctions 
of H^ u . Therefore, the wave functions for the octet or decuplet baryons after the diagonalization do 
not only have an admixture of the representations mentioned above, but contain representations of 
arbitrary high dimensions. Of course, if the symmetry breaking can be considered as small, these 
higher representations are also suppressed in the wave functions. Nevertheless, this method is able 
to calculate the strange mass corrections from the wave function to all orders. On the other hand, as 
we have seen above, there are contributions of the same order from expanding the baryon correlator 
itself. Therefore it seems to be more consistent to calculate all kinds of strange mass corrections 
up to the second order instead of calculating only one of these corrections up to all orders and the 
others up to the second. The fact that the latter ones are calculated only up to the second order is a 
technical problem to which no solution is seen. 

However, apart from these theoretical considerations a numerical comparison between the Yabu- 
Ando treatment and the perturbative treatment of the wave function correction up to the second 
order has been performed recently (Blotz et al, 1995a). As can be seen clearly in this investigation 
both methods lead to almost identical results for the mass splittings. This means that higher order 
corrections, at least for the wave function, O (m s ) 3 ) do not play a significant role. 



Table 5: Hyperon mass splittings in second order perturbative calculations in the SU(3) NJL model 
(Schneider et al, 1995). 
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594 


Q 


725 


733 



5.10 Results: Mass splittings in the SU(3) NJL model 

Several authors have performed numerical calculations for the hyperon mass splittings. First order 
perturbative results can be found in (Weigel et al, 1992, Blotz et al, 1992, 1993a). These papers 
present also the numbers obtained in the Yabu-Ando approach which has been compared to the the 
second order results by (Blotz et al, 1995a). The inclusion of the second order corrections in the 
strange mass improves the agreement with the experiment and the corresponding second order results 
from (Schneider et al., 1995) are presented in the Table |5|. 

It is also interesting to see to what extent the numbers of table |5] fulfill the the Gell-Mann - Okubo 
mass relations. If the nucleon mass in the NJL model is put to 938 meV one obtains for | + 
|f«s — \{rriN + ma) a value of —9.2 MeV (from table ^) and +6.7 MeV for the experimental masses. 
Both values of these violations are on the scale of 1 % in the baryon masses, nevertheless one should 
note the different sign. For the baryon decuplet ([]]) the three terms in ([!]) are 138, 149 and 150 MeV 
from table [5] compared to 139, 149 and 152 MeV for the experimental masses. Apparently, both the 
experiment and the SU(3) NJL model calculations fulfill the 11-Mann - Okubo mass relations equally 
well. 

The above mass splittings of the baryons belonging to the octet and decuplet are due to the finite 
value of m s . It is also interesting to calculate the mass splittings within the isospin multiplets. It 
consists of the hadronic and the electromagnetic part 

AM B = (AM B ) h + (AM B ) e . (290) 

The origin of this is theoretically clear, since the electromagnetic part is due to the baryon electromag- 
netic self-energy and the hadronic part is due to the current mass difference Sm u d = vn^ — m u . Gasser 
and Leutwyler (1982) quote the following values for the experimental estimates of the electromagnetic 
parts of the octet baryons (in MeV): 

(n-p) e = _o.76 ±0.3, (E~ - = 0.17 ±0.3, (s~ - = 0.86 ± 0.3 . (291) 

These values have to be subtracted from the measured mass differences 

i n -P)exv = L29 > =8.07 ±0.09, (S--E°) =6.4 ±0.6, (292) 

V / exp \ / exp 

and one obtains for the hadronic parts 



) =5.5 ±0.7. 



(293) 




m d - m u [MeV] 



Figure 18: Hadronic part of the octet isospin splittings N(n — p), — H°), £(X~ — £ + ). The 

SU(3)-NJL calculations are presented by solid lines (with m s - corrections) and dashed lines (without 
m s -corrections). The experimental values including error are indicated by dashed-dotted lines. The 
double lines indicate the regions where the SU(3)-NJL calculations with m s -corrections overlap with 
the experimental numbers. 

Unfortunately, these simple estimates do not exist for the baryon decuplets. In chiral quark models, 
the hadronic part of the mass splitting is related to the current quark mass difference Am B( j = m^— m u 
and in Skyrme type models to the corresponding boson fields. Actually, in the SU(2) NJL the hadronic 
part of the neutron-proton mass difference vanishes identically in the leading order in iV c . One way to 
cure this disease is to enlarge the symmetry group. Jain et al. (1989) used a U(2) ® U(2) extension 
of the Skyrme model with pseudoscalar and vector fields. They predicted, however, a n-p hadronic 
mass difference which was 35% too small. Some increase was obtained by the same authors in the 
case of U(3) <8> U(3). Similar results were obtained in the chiral bag model (Park and Rho, 1989), 
where quark degrees of freedom are explicitly taken into account. The symmetry breaking in the 
Skyrme model has recently successfully been reinvestigated by Walliser (1993). 

In contrast to some of these estimates, the SU(3) NJL calculations of Praszalowicz et al. (1993) 
and Blotz et al. (1994) show good agreement with the experimental data not only for the n-p 
system but also for X - — S + and H~ — 5°. In these calculations, the hadronic part of the isospin 
splittings is evaluated perturbatively to orders 0{m 2 s ) and 0(m s Am u d). The results for a constituent 
mass of M — 420MeV they obtained are shown in fig. [H| As can be seen, for a common value of 
Am u d of 4.4MeV the experimental data are reproduced within the error bars. If one ignores strange 
mass corrections, the difference Am M d reduces to about 3.5MeV. Obviously the SU(3) NJL offers a 
satisfactory description for the isospin mass differences in the octet. Predictions for the decuplet can 
be found in (Praszalowicz et al, 1993). 



6 Baryon properties in the SU(3) NJL model 



In the previous chapter, we have presented the collective quantization in SU(3) case and showed that 
the SU(3) model provides a successful description of the mass splittings of the SU(3) baryon octet and 
decuplet. From the mass splittings for the constituent quark mass M = 420 MeV has been chosen 
and this value will be used throughout all calculations, which we will present in this subsection, so 
that there will be no more free adjustable parameter. Since the SU(3)-formalism is a generalization 
of the SU(2) NJL model, though being quite involved and complicated, we will not repeat the details 
of the calculations but concentrate on the results for the static properties as well as for various form 
factors of the SU(3) octet baryons. We will consider the contributions up to the first order in both Q 
and A Ms (m s ). 

The investigations in SU(3) case has two significant meanings: First, stimulated by recent experi- 
ments, it is now an important issue of hadron physics to understand the role of the strangeness in 
the nucleon structure. Second, to study the properties of the other members of the SU(3) multiplets. 

We will start with studying the electromagnetic properties of the SU(3) octet baryons. 



6.1 Electromagnetic properties 

The electromagnetic Sachs form factors Ge{<1 2 ) and Gm{<1 2 ) 

(B'(p>)\f(0)\B( P )} = G E ( q 2 )5j, j3} (294) 

(B'(p')\j k (0)\B(p)) = ^ klm (T l )j^ m G M (q 2 ) , (295) 
are related to the matrix element of the time and space components of the electromagnetic current 

j^x) = #(z) 7 „Q*(a;) . (296) 
Here denotes the experimental nucleon mass and Q is the quark charge matrix 



Q = -| (297) 




The q is the four momentum transfer and we also use Q 2 = —q 2 . Hence, the electromagnetic current 
jn can be decomposed in third and eighth SU(3) octet currents 

Jm = jf + 7^4 8) ( 298 ) 

with 

jf = and jf = ^A 8 *. (299) 

The matrix element of the electromagnetic current are evaluated according to the technique developed 
in chapter 3 (for details, see (Kim et ai, 1995a)). 



Fig. |20] shows the electromagnetic form factor of the nucleon for three different values of the constituent 
quark mass M (Kim et al. 1995a). As can be seen, the model results agree very well with the presented 
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Figure 19: Electric (left) and magnetic (right) form factors of the nucleon as functions of Q 2 evaluated 
in the SU(3) NJL model for three different values of M: M = 370MeV (dashed line), M = 420MeV 
(solid line) and M = 450MeV (dash-dotted line) (Kim et a!., 1995a). For the strange mass m s = 
180 MeV is used. The experimental data are from (Hohler et ai, 1976, Platchkov et ai, 1990, Eden 
et ai, 1994, Meyerhoff et al, 1994, Bruins et ai, 1995). The results for the magnetic form factors 
are presented without any scaling. 



experimental data. In particular, the g-dependence of the form factors are reproduced rather well. 
The deviations from the experimental data, seen in the case of the magnetic form factors, is due to a 
slight underestimation of both proton and neutron magnetic moments. It should be also mentioned 
that the electromagnetic form factors depend weakly on the particular values of the constituent 
quark mass if this is chosen between 400 and 450 MeV. As in the case of SU(2) version, values of 
the constituent mass around M = 420 MeV seem to be preferable. In the case of the magnetic form 
factors on fig. [2(], one should note that the results are presented without any additional normalization 
to the experimental values of the magnetic moments. 

On fig. [H] we plot again the neutron electric form factor with and without the strange mass corrections 
compared to the experimental data as well as to the results of SU(2) case. As can be seen, the m s 
corrections influence the results at low momentum transfer and in particular, the neutron charge 
radius. It should be also noted that in SU(3) case the results for the neutron electric form factors 
are quite different from those of the SU(2) version and are in a better agreement with experiment. 
In fact, this difference can be understood looking at the structure of the quark charge matrix. In 
contrast to SU(2) case where Q — |(T = 0) + y(T = 1), i.e. the isoscalar part is a constant, in SU(3) 
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Figure 20: Neutron electric form factor as a function of Q 2 evaluated in the SU(3) NJL model (Kim 
et al., 1995a) with m s = 180 MeV (solid line) and without m s corrections (dashed line) in comparison 



with the SU(2) result (dash-dotted line). The experimental data are the same as in fig. [LS. 



it includes As, Q = 57/3 = 0) + 4f(T = 1). Thus the isoscalar part of electromagnetic current ( |296| ) 
is not invariant under rotation R. This leads to a reduction of the electric isoscalar form factor in 
SU(3) and results in considerably smaller numbers compared to those of the SU(2) model. 



On fig. |21] and fig. | 
(Kim et al, 1995a). 



we present the electromagnetic form factors for the SU(3) octet hyperons from 



Without m s corrections, for all momentum transfers the form factors fulfill the relations 
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(300) 



which follow from the ?7-spin symmetry. 



In table || we present the magnetic moments and the charge and magnetic squared radii of the SU(3) 
octet baryons. The results agree with the experimental data within 20%. It should be mentioned that 
Bae and McGovern (1996) made a x 2 analysis of the hyperon magnetic moments for effective theories 
and found that the present NJL model provides the best description among the hedgehog approaches 
to these quantities. It should be mentioned that the SU(3) version provides a better description of 
the nucleon magnetic moments compared to the SU(2) one. 
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Figure 21: Electric form factors of charged (left) and neutral (right) octet baryons as a function of 
Q 2 in the SU(3) NJL model for M = 420 MeVwith m s corrections taken into account (Kim et al, 
1995a). 
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Figure 22: The same as in fig. 21 but for magnetic form factors (Kim et al, 1995a). The experimental 
values for the magnetic moments are also shown (full circles). 



Table 6: Electromagnetic properties of the octet baryons in the SU(3) NJL model (Kim et al. 1995e). 



Baryons 


< r > e |tm J 


exp. 




n.m.} 


exp. 


< ^ >m [mi J 


exp. 


V 


n 7Q 
U. to 


U. / 4 


2.39 


1. t\) 


n 7n 
U. (0 


U. /4 


n 


n no 

— u.uy 


n 1 1 
— U.li 




1.76 


1 01 


U. 


n 77 


A 


n ri/i 
— U.U4 






0.77 


— U.D1 


U. ( U 




Y+ 


fl 7Q 




2.42 


9 4fi 


71 

U.I 1 




s° 


0.02 




0.75 




0.70 




E- 


-0.75 






0.92 


-1.16 


0.74 






-0.06 






1.64 


-1.25 


0.75 






-0.72 






0.68 


-0.65 


0.51 





6.2 Strange vector form factors of the nucleon and related observables 



We now proceed to discuss the strange vector form factors which recently have attracted much 
attention. The success of the model in describing the electromagnetic properties of the nucleon 
encourages one to apply it for studying also the strange vector form factors. Experimentally, although 
several experimental proposals have been suggested and a series of experiments is in progress (Musolf 
et al, 1994), there is still no experimental result for the strange vector form factors available. From 
the side of the theory the strange form factors of the nucleon have been studied within different 
model schemes. Using pole fit analyses based on dispersion relations, Jaffe (1989) first calculated the 
strange form factors and gave estimates for the squared strange radius and magnetic moment of the 
nucleon. These quantities have been also studied in the SU(3) Skyrme model without (Park et al, 
1991) and with vector mesons (Park and Weigel, 1992), in a kaon-cloud model including kaon-loops 
(Musolf and Burkardt, 1994) and also based on the vector dominance (Cohen et al, 1993). Forkel 
et al, (1994) have evaluated the strange vector form factors using different models, in particular, 
a vector dominance model with uj — 0-mixing and a kaon cloud. Hammer et al, (1996) updated 
the Jaffe's results using a new dispersion theoretical analysis of the nucleon electromagnetic form 
factor (Mergell et al, 1996). In a very recent work, using some QCD equalities for the baryon current 
elements Leinweber (1995) estimated the strange magnetic moment of the nucleon. In the NJL model, 
including both the m s and the rotational 1/N C corrections the strange vector form factors have been 
studied recently by Kim et al, (1995c). 

The information of the strange vector form factors in the nucleon is contained in the quark matrix 
elements: 

(N(p')\J^\N(p)) = (N(p')\s^s\N(p)). (301) 

The strange quark current j s can be expressed in terms of the baryon current and the hypercharge 
current jY: 

Jl = s7m s = JfWvV ~ T^T/A? = il ~ jl (302) 
Here, the sign convention of Jaffe (1989) is used. 

As in case of the electromagnetic form factors, Sachs isoscalar and isovector form factors G s E (q 2 ) 
and G s M (q 2 ) 

(N'(p')\Js(0)\N(p)) = G s E (q 2 ) 

(N'(p')\jn0)\N(p)) = ^-G s M (q 2 )te tjkq ^s'\a k \ S ) (303) 
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Figure 23: Strange electric form factor G S E evaluated in the SU(3) NJL model (Kim et al, 1995c) as 
functions of Q 2 without m s (left) and with m s corrections (right) for different values of the constituent 
quark mass M: 370 MeV (dash-dotted line), 420 MeV (solid line) and 450 MeV (dashed line). 



can be defined related to the matrix element of the time and space components of the strange current, 
respectively. However, in contrast to the isoscalar electric form factor, which at q 2 = determines 
the baryon number B = 1, the strangeness of the nucleon is zero, G S E (0) = 0. 



Fig. 1241 shows the strange electric form factor G S E with the constituent quark mass M between 370 
and 450 MeV calculated in the NJL model (Kim et al., 1995c). As can be seen, the strange electric 
form factor G S E depends weakly on M but it is rather sensitive to the m s corrections. The inclusion 
of these strange mass corrections leads to a strong enhancement of more than 70 % of the form factor 
at finite momentum transfers. Similarly to Forkel et al, (1994), the strange electric form factor G S E 
in the NJL model (Kim et al, 1995c), has a sign opposite to the one of the pole- fit analyses of Jaffe, 
(1989) and Hammer et al, (1996). 



The strange electric squared radius is given by 

(r 2 )s = ~ 



JG%(Q 2 ) , 

b dQ 2 W=°- 



(304) 



and the results from different model calculations are summarized in table 0. 



In the lack of experimental data, it is useful to compare the results of different model calculations 
presented in table 0. In fact, only two approaches, pole-fit analysis (Jaffe, 1989, Hammer et al, 1996) 
and NJL with m s and rotational corrections (Kim et al., 1995c), predict clear non-zero values for 
the strange radius but of opposite signs which simply reflect the opposite signs of the correspond- 
ing strange electric form factors. It should be stressed that both approaches incorporate a proper 
description of the nucleon electric form factors and the corresponding charge radii. 



In fig. E4] we show the strange electric density from (Kim et al 
charge densities 



1995c). As in the case of the neutron 



(see fig. [I~0| ), this density has a negative tail which is due to the contribution of 
the polarized Dirac sea. In fact, the shape of this tail is determined by the behavior of quark wave 
function at large distances. Because of the embedding ansatz ([]]) it is simply given by the pion 



Table 7: Comparison of the strange magnetic moments and strange electric squared radii calculated 
in different models. In the SU(3) NJL calculations a constituent mass M = 420 MeV is used (Kim 
et al, 1995c). 
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references 


pole-fit 


-0.31 ±0.09 


0.14 ±0.07 


Jaffe 89 


pole-fit 


-0.24 ±0.03 


0.19 ±0.03 


Hammer 95 


kaon cloud 


-(0.31 -> 0.40) 


-(0.027^ 0.032) 


Musolf & Burkardt 94 


VMD 


-(0.24 -> 0.32) 


-0.04 


Forkel 94 


Skyrme 


-0.13 


-0.11 


Park 91 


Skyrme with VM 


-0.13 


0.05 


Park & Weigel 92 


QCD equalities 


-0.75 ±0.3 




Leinweber 95 


NJL (ft + ft 1 ) 


-0.45 


-0.35 


Kim 95 




Figure 24: Strange electric densities split in valence and Dirac sea contributions in the SU(3) NJL 
model (Kim et al, 1995c). The sign convention of Jaffe (1989) is used. 
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Figure 25: Strange magnetic form factor G S M (Kim et at, 1995c) as functions of Q 2 evaluated in 
the SU(3) NJL model without m s (left) as well as with m s corrections (right) for different values of 
the constituent quark mass: M = 370 MeV (dash-dotted line), 420 MeV (solid line) and 450 MeV 
(dashed line). 



mass, exp(— m^r). On the other hand, it is generally expected that the tail of the strange electric 
density should be related to the kaon cloud with asymptotics exp(— m^r). Hence, one can expect 
some overestimation for the quantities related to the strangeness in the NJL approach for which the 
contribution of the Dirac sea is dominant. In particular, this caveat applies to the strange radius. As 



one can see in fig. |23| the present model yields a valence contribution to the strange charge distribution 
localized in the interior of nucleon. In terms of a baryon-meson picture this corresponds to a A or X 
core coupled to the corresponding kaon cloud. 

The strange magnetic form factor (Kim et ai, 1995c) is plotted on fig. As in the case in the case 
of the strange electric form factor the magnetic one shows a weak dependence on the constituent mass 
whereas the inclusion of the m s corrections leads to a common enhancement of more than 50 %. The 
results are similar to those of Jaffe (1989), Hammer et al. (1996) and Forkel et al. (1994). 

The strange magnetic moments for different models are given in table [7]. All presented model cal- 
culations predict a negative sign and similar values for We remind the reader that both the 
pole-fit analyses (Jaffe, 1989, Hammer et al., 1996) and NJL model (Kim et al., 1995c) provide a 
proper description of both proton and neutron magnetic form factors and the corresponding magnetic 
moments. 

There are also some results of strange form factors using the NJL model performed by Weigel et 
al., (1995a). The authors have used non-selfconsistent profiles and have ignored the rotational 1/N C 
corrections. In addition, there seem to be also some purely numerical shortcomings of the calculations 
of Weigel et al, (1995a). 



6.3 Scalar Form Factor and ttN sigma term 



The scalar form factor a(t) is defined as 



a(t) = m (N{p')\uu + dd\N{p)) (305) 

with m = (m u + m^)/2 ~ 6 MeV. The mass m is fixed using the additional free parameter in 
the proper-time regularization integral Here we use the usual notation t for the square of the 

momentum transfer. 



Because of a factor-of-two discrepancy between empirical value for the ttN sigma term (S = 64 ± 8 
MeV) (Koch, 1982) from the analyses of ttN data and naive estimates of the cx-term (a ~ 25 MeV) 
from the baryon mass spectrum (Cheng, 1976), which contradicts to the prediction of the low-energy 
theorem of current algebra, there have been a lot of discussion in the literature (see (Gasser et al, 
1991, Gasser and Leutwyler, 1991) and references therein) devoted to this problem. Donoghue and 
Nappi (1986) suggested that the discrepancy is due to the presence of strange quarks in the nucleon, 
i.e. (iV|ss|./V) ^ 0. However, using both the Skyrme model and bag model they got unlikely large 
(N\ss\N) contribution (almost 30%) to the quark condensate in the nucleon. Gasser et al. (1991) 
reanalysed the a term carefully taking into account all currently available nN scattering data. They 
concluded a stronger t-dependence of the scalar form factor a(t) resulting in 

A CT = a(2ml) - a(0) ~ 15 MeV, 

and a larger value for sigma term cr(0) = 45 ± 8MeV which means that the value at the Cheng- 
Dashen point a(2m 2 T ) = E ~ 60MeV. This value differs from the estimates from the baryon spectrum 
and it is due to both non-zero strange current mass and the strange content (iV|ss|./V) in the quark 
condensate. For y = 2(N\ss\N) / (N\uu + dd\N), the strange contribution to the a term, Gasser et al. 
(1991) suggested a value of about 0.2, which means that the strange contribution to (N\uu+dd+ss\N) 
should be about 10 % and hence much smaller than suggested by Donoghue and Nappi (1986). 



Fig. [26] shows the scalar form factor as a function of the constituent quark mass M (Kim et al, 1996). 
As can be seen, the scalar form factor is very sensitive to the the particular value of M only in the 
time-like region close to the Cheng-Dashen point. For M = 420 MeV both the empirical estimate 
cr(0) = 45 ± 8 MeV and the deduced t dependence A a = 15.2 ± 0.4 MeV as well as the slope at zero 
momentum transfer (strange scalar ms radius (r 2 )^ ~ 1.6fm 2 ) (Gasser et al, 1991) are reproduced. 
The actual values of (Kim et al, 1996) are a = 41 MeV, A a = 18.2MeV and (r 2 ) s N = 1.5 fm 2 . It 
should be noted that the m s corrections do not influence a(t) much. However, the m s corrections 
play an important role of reducing strongly the strangeness contribution (iV|ss|iV). With the m s 
corrections taken into account, Kim et al. (1996) obtain y = 0.27 in case of M = 420MeV. This 
qualitatively agrees with the value y ~ 0.2 of Gasser et al. (1991), whereas without the m s corrections 
y it is strongly overestimated y = 0.48. 



6.5 Axial charges of the nucleon 



The axial charges of the nucleon are defined as forward matrix elements of the axial current 

(N(p)\^ 5 \ a V\N(p)) = <7i o) tZ(p) TM 75«(p), (306) 

where a = 0, 3, 8 and A 3 , A 8 are Gell-Mann matrices, A is in this context 3x3 unit matrix and u(p) 
is a standard Dirac spinor. 

(3) (S) 

The flavor non-singlet axial charges g A and g A can be measured in the experiments on hyperon 
semileptonic decays. In particular, the g A is known with good accuracy from the neutron /3-decay 
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Figure 26: Scalar form factor as a function of Q 2 evaluated in the SU(3) NJL model (Kim et al, 
1996): On the left for different values of the constituent up and down quark mass 370 MeV (dotted 
line), 420 MeV (solid line) and 450 MeV (dashed line) and with m s corrections, and on the right, for 
M = 420 MeV and with (solid line) and without (dashed line) m s corrections. The empirical value 
for cr(0) = 45 ± 8 MeV of (Gasser et al, 1991) is also presented. 



data (Particle Data Group 1994): 



g f =g A = 1.257 ± 0.0028. 



(307) 



In order to extract g A from experiment an additional assumption of flavor SU{3) symmetry is usually 
used for analysing the data on hyperon /3-decays. In particular, Hsueh et al. (1988) give the following 
estimate 



gf = 0.34 ± 0.02. 



(30£ 



and for a recent discussion of the applicability of flavor SU(3) symmetry see (Lichtenstadt and 
Lipkin, 1995). The flavor singlet axial charge g A cannot be obtained directly from hyperon weak 
decay experiments, but it can be related to the first moment of the polarized quark distribution in 
the proton. The latter is defined as 



9i( x ) 
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q=u,d,s 



2 E 4 A 9(x) 

q=u,d,s 



(309) 



In the parton model, the moment Aq = Jq dxAq(x) is related to the light quark (of flavor q) con- 
tribution to the proton spin. The distribution ( p09| ) has been measured in deep inelastic scattering 
experiments. Using ( |309| ) and the relations 



gf = Au-Ad 
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Figure 27: Proton spin fraction carried by light quarks as a function of the soliton size R for a fixed 
profile function of the SU(2) NJL model. 



gf = (Am + Ad — 2As) 

n {0) 
9a 



Au + Ad + As = AE. 



(310) 



the moments of separate polarized quark distributions Au, Ad and As as well as the singlet axial 
charge g^' can be obtained from the experimentally measured moment Jq dxg\ and axial charges 
and g^ . The 1988 EMC result (Ashman et a/.1988) suggested that the flavor singlet axial charge g^ 
is close to zero which indicates that none of the proton spin is carried by its constituent quarks. In 
addition the EMC data disagreed with the assumption As = 0. This has sometimes been called the 
"proton spin crisis" . This EMC result was rather surprising from the point of view of the constituent 
quark model (CQM), which would suggest = 1 and As = 0. Soon after the first data became 
available it was demonstrated by Brodsky et al. (1988) that in the Skyrme model the singlet axial 
charge is suppressed by l/N c in comparison with CQM and this was suggested as a solution of the " 
spin crisis". 

In the NJL model, the g^' has a contribution which originates from the imaginary part of the effective 
action and hence it has no counterpart in the Skyrme model (Blotz et al. 1993b). An analogous 
contribution is also pertinent to the Skyrme model with explicit vector mesons (Park and Weigel, 
1991), though its origin and particular value differ from that of NJL model. It should be noted that 
in the large N c limit this NJL contribution is non- vanishing. However, it is suppressed in large soliton 
size R (or large constituent quark mass) limit g^ ~ 1/(MR) 4 . In fig. |27] the proton spin fraction 
carried by light quarks |AE = in dependence of MR is plotted. At small Mi?, corresponding to 
CQM, it is close to one half whereas it is rapidly vanishing with increasing MR, so that in the limit of 
large MR the Skyrme model result g^ = is reproduced. This figure illustrates how the NJL model 
in some sense interpolates between two models of baryons, namely the constituent quark model and 
the Skyrme model. Indeed, if MR is small, the valence level is situated just below the threshold 
M, and the picture resembles the non-relativistic constituent quark model. In the opposite case of 
MR 3> 1 the NJL picture of nucleon appears to be ideologically close to that of the Skyrme-Witten 



Table 8: Axial constants and the light quark contributions to the proton spin calculated in the SU(3) 
NJL model (Blotz et ai, 1993b, 1996) compared with the experimental values taken from the analysis 
of Ellis and Karliner (1995). 
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9 A 
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9 a 




Au 


Ad 


As 


NJL SU(3) 


1.38 


0.30 


0.37 


0.90 


-0.48 


-0.05 


exp. 


1.26 


0.34 ±0.02 


0.31 ±0.07 


0.83 ±0.03 


-0.43 ±0.03 


-0.10 ±0.03 



topological soliton. Actually, the truth seems to be just between these two cases, because in the 
NJL model for the selfconsistent solution, which provides an overall good description of the baryon 
properties, MR « 1 . 

The results of the selfconsistent calculations of the axial charges and g^' in the three flavor 

NJL model are summarized in table |8| (Blotz et al. 1993b, 1996a). These results are obtained in the 
SU(3) NJL model with both linear in m s and rotational 1/N C corrections taken into account with a 
constituent quark mass M = 420 MeV. 

As can be seen from the table, the values of axial charges calculated in the three flavor NJL model 
are in a good agreement with experimental values from hyperon f3 decays ( Hsueh et al. 1988) and 
recent analysis of experimental data on polarized deep inelastic scattering (Ellis and Karliner, 1995). 

In QCD, the value of is related to the nucleon matrix element of the topological charge, (N\FF\N), 
by the U(1)a anomaly. Recently, using the instanton model of QCD vacuum, a method has been 
developed to represent systematically QCD gluon operators as constituent quark operators in the 
NJL model derived from the instanton model. (Diakonov et ai, 1995b). It was shown that this 
method fully preserves the U(1)a anomaly. Thus, in the context of the instanton vacuum, the result 
gjP = 0.37 of (Blotz et ai, 1993b) can be regarded as a consistent estimate of the nucleon matrix 
element of FF. 



6.6 Gottfried sum rule 



The Gottfried sum rule, which is the unpolarized analogue of the Bjorken sum rule for the polarized 
structure functions, was reexamined in the recent NMC experiment (Arneondo et al., 1994) and 
the result of the measurement is again in contradiction to the naive quark model expectation. The 
Gottfried sum is defined as 

S G = f dx- (Fi(x) - F2{x)) 
Jo X 

where the F^ \x) are the unpolarized structure functions. 

In the language of the parton model and in terms of quark distribution functions: 

S G = fdx £ el (qf + qf - q? - %) (311) 

Jo i=u,d 

= [ dx- (u + u-d-d) (312) 
Jo 3 v ' 

where in the last line the SU(2) flavor symmetry is used. The u, u, d and d are the quark distribution 



functions of proton. Using standard quark sum rules for the valence part of the proton 



J dx(u — u) = 2, J dx(d — d) = 1, 



the Gottfried sum ( |312| ) becomes 

S ° = I + \ I dX ^ ~ ^ 

where the second part in rhs of (^) is related to the flavor asymmetry of the sea. If it is flavor 
symmetric one obtains the Gottfried sum rule 

Sg= 1- 

However, the recent NMC data (Arneondo et ai, 1994) leads to a different value 

S G = 0.221 ± 0.008 
which suggests that the Dirac sea carries some flavor asymmetry. 

Assuming a particular ansatz (Wakamatsu, 1992) for the integrated distributions (second term in (^)) 
Wakamatsu (1992) in SU(2) and Blotz et al. (1996b) in SU(3) version analyzed the Gottfried sum 
rule in terms of the flavor asymmetry of the Dirac sea in the NJL model. The obtained estimations 
are given in table |6|. In the SU(3) case the m s corrections are also included (Blotz et ai, 1995b). 
It should be noted that in the NJL model the deviation from of So = 1/3 occurs even in the chiral 
limit. 



Table 9: NJL estimates of (Wakamatsu 1992) and of (Blotz et ai, 1995b), respectively, for the 
Gottfried sum Sq for SU(2) and SU(3). The values are given for M = 420 MeV in the chiral limit as 
well as with m s corrections. 



M [MeV] 


Q SU(2) 


b G no m s 


S G U ^ with m s 


Exp 


420 


0.2091 


0.2336 


0.2341 


0.221 ± 0.008 



7 Summary 



The present article reviews the theory and applications of the Nambu-Jona-Lasinio model with two 
and three flavors for the description of low-mass baryons. The model corresponds to an effective chiral 
theory where quarks interact with selfconsistently generated bosonic mean fields. After concentrating 
on Goldstone fields (pion, kaon and eta) and after semiclassical quantization of the collective rotational 
degrees of freedom the model coincides with the Chiral quark soliton model of Diakonov, Petrov and 
Pobylitsa obtained from QCD assuming an instanton liquid for the QCD vacuum. After adjusting 
the parameters of the model in the meson sector to the pion decay constant and the pion and kaon 
masses, the constituent quark mass is left as the only free parameter. This parameter is tuned in 
the baryonic sector so that the model using a hedgehog form for the meson fields reproduces best 
the wealth of experimental data with a value of the constituent quark mass M around 420 MeV. 
With this uniquely defined parameter set the calculations are performed for various radii, magnetic 
moments, mass splittings of the baryon octet and decuplet, and axial and spin properties of the 
nucleon. Particular emphasis is put on electromagnetic and axial form factors. Predictions are made 
for strange form factors and tensor charges. Without any readjustment of the parameters, a very good 
agreement between the theory and nearly all low-energy experimental data is obtained. Therefore, 
one can consider the present model as a quite useful tool to describe the ground-state properties of 
the octet and decuplet baryons. 

For the parameter set associated with a constituent quark mass of 420 MeV, the practical outcome 
of the calculations can be summarized as follows. SU(2) case: Within 15 % we reproduce the proton 
charge radius, the axial vector coupling constant, the pion-nucleon sigma term, the delta- nucleon 
splitting, the E2/Ml-ratio, the electric proton form factor and the momentum dependence of the 
the magnetic form factors of proton and neutron. The nucleon magnetic moments deviate by 25 % 
from the experiment. The electric neutron form factor and the corresponding charge radius appear 
to be overestimated by 50 %. Predictions are given for tensor charges. SU(3) case: Within 15 % we 
reproduce the proton charge radius, the axial vector coupling constant, the pion nucleon sigma term, 
the delta-nucleon splitting, the mass splitting within the octet and decuplet, the isospin splittings 
within the octet. The electric proton form factor and the magnetic form factors of proton and 
neutron and the corresponding magnetic moments deviate also less than 15 % from experiment. The 
electric neutron form factor and the neutron charge radius are by 20 % too small compared to the 
experimental data. The flavor singlet and octet axial constants come out within the experimental 
errors. These quantities are connected with the "spin crisis" of the proton. Predictions are given 
for strange electric and magnetic form factors of the nucleon and for electromagnetic form factors of 
octet and decuplet baryons. The corresponding magnetic moments agree well with experiment. 

The present model reproduces most of the relevant low-energy observables of the nucleon and of 
various hyperons. It is important to note that this agreement is rather good and hence the model 
might be used to obtain some insight into the microscopic structure of the baryons. The structure 
of the baryons emerges to be based on chiral symmetry breaking as a dominating mechanism and 
the important resulting forces are those between quarks and Goldstone bosons. Confining forces, 
being absent in the model, seem not to play a relevant role for the baryonic ground states. Also 
vector mesons and scalar mesons are apparently not needed. This simple picture of the interactions 
is associated with a simple picture of the microscopic structure of the baryons. They appear as 
composed of three valence quarks in a bound state together with a continuum of the polarized Dirac 
sea. In fact, the contribution of the polarized Dirac sea to the baryon observables is at about the 
30 % level. It should be stressed that the polarized Dirac sea contribution is absolutely necessary 
since it causes the valence quarks and hence the whole baryon to be bound. The baryon number 
of this system is given by the occupation of the valence level and not by the topological properties 
of the meson fields, as it is in the Skyrme model. For the parameter values used, an approximation 



of the present non-local NJL effective action by a local Skyrme-type action is not justified. Thus, 
the present model picture of the baryon with clearly distinguished valence quarks differs from the 
topological soliton and hence, the present approach does not support topological baryon models like 
the Skyrme one. 

The above considerations are all made within the non-linear version of the Nambu-Jona-Lasinio 
model (on the chiral circle) being in this form equivalent to the Chiral quark soliton model. There 
are several problems still within this model, which need further consideration. One consists in the 
undesired asymptotic form exp(— m^r) of the meson profile in the strange direction (kaon field) as it 
turns out in the present treatment of the strange mass corrections. This affects the applicabilty of the 
formalism to those "strangeness" quantities which strongly depend on the tail of the polarized Dirac 
sea, as e.g. strange radii. Another problem is related to the present large N c treatment in which only 
the l/N c corrections originating from the zero modes are taken into account. This question is closely 
connected with the fact that the absolute masses of the of the baryons come out noticeably too high 
and the meson-loop corrections, being presently very much discussed in other chiral effective theories, 
are needed to provide a remedy for this. In the present model it seems to be quite a formidable task, 
both formally and numerically. 

There is also an interesting field of applications of the model that is related to gluonic degrees of 
freedom in baryons. For example, the evaluation of matrix elements of operators related to higher twist 
contributions to moments of nucleon structure functions has recently attracted much attention. The 
Nambu-Jona-Lasinio model as such, based on a four quark interaction, cannot be used here since it 
does not contain gluon degrees of freedom. However, the non-linear version of NJL model, having been 
bosonized and restricted to Goldstone modes by invoking the chiral circle, corresponds directly to the 
effective quark meson action obtained in the instanton liquid model of QCD. One therefore has a link 
between the description of baryons, as it is presented in this review, and the instanton configuration 
of the gluon field. Making use of that, the present approach has the potentiality systematically 
to represent gluon operators as effective quark operators, so that baryon matrix elements of gluon 
operators can be evaluated using the techniques presented in this review (Diakonov et al, 1995). In 
this way useful information for the understanding of nucleon structure functions may be obtained. 
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